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Abstract 

Magnetic edge states are responsible for various phenomena of magneto-transport. 
Their importance is due to the fact that, unhke the bulk of the eigenstates in a mag- 
netic system, they carry electric current along the boundary of a confined domain. 
Edge states can exist both as interior (quantum dot) and exterior (anti-dot) states. 
In the present report we develop a consistent and practical spectral theory for the 
edge states encountered in magnetic billiards. It provides an objective definition for 
the notion of edge states, is applicable for interior and exterior problems, facilitates 
efficient quantization schemes, and forms a convenient starting point for both the 
semiclassical description and the statistical analysis. After elaborating these topics 
we use the semiclassical spectral theory to uncover nontrivial spectral correlations 
between the interior and the exterior edge states. We show that they are the quan- 
tum manifestation of a classical duality between the trajectories in an interior and 
an exterior magnetic billiard. 

PACS: 05.45+b, 73.20.Dx, 03.65.Sq, 03.65.Ge 

Keywords: magnetic billiards; edge states; quantum chaos; spectral correlations 



Preprint submitted to Elsevier Preprint 



mpipks/0112003 



Contents 



1 Introduction 



2 Motion in the free magnetic plane 



2.1 Classical motion 



2.2 Quantization] 



2.3 The scaling property 



2.4 The free quantum propagator 



2.5 The free Green function 



3 Introducing a boundary 



3.1 Motion in a restricted domain 



3.2 The classical billiard 



3.3 Quantum billiard; 



3.4 Orbital magnetism 



4 Quantization in the interior and the exterior 



4.1 Boundary methods 



4.2 The boundary integral equations 



4.3 The boundary operators 



4.4 Solving the integral equations 



5 Results of the boundary integral method 



5.1 Spectral statistics 



5.2 Wave functions in the interior and in the exterior 



5.3 General boundary conditions 



6 Semiclassical Quantization 



6.1 The semiclassical boundary integral operators 



6.2 From boundary to map operators 



4 
7 
7 

10 
14 
16 
18 
24 
24 
26 
31 
37 
41 
41 
42 
47 
50 
50 
51 
54 
60 
63 
64 
68 



2 



6.3 Trace formula for hyperbolic billiards 



6.4 Trace formula for the integrable cas^ 



6.5 The separable case 



7 A spectral measure for edge states 



7.1 Bulk states and edge states 



7.2 A spectral density based on the boundary conditions 



7.3 Asymptotic properties of edge and bulk states 



7.4 Edge magnetization as a spectral measure 



8 Properties of edge state spectra 



8.1 Universal auto-correlation; 



8.2 The action spectrum 



8.3 Using the edge magnetization 



9 Spectral cross correlations: The interior-exterior duality 



9.1 A semiclassical theory of spectral cross correlations 



9.2 Numerical evidence 



9.3 The pair relation 



10 Conclusions 



Acknowledgements 



A Appendices 



A.l Green function in angular momentum representation 



A. 2 The null field method 



A. 3 Exact quantization of the magnetic disk 



A.4 The stationary phase approximation 



A.5 The product relation of the map operators 



A.6 The straight line with periodic boundary conditions! 



71 
84 
89 
96 
96 
99 
107 
111 
115 
115 
119 
123 
124 
125 
130 
133 
141 

142 
142 
143 
145 
146 
148 
149 
151 



A.7 Scaled spectra 



155 



List of Symbols 



156 



References 



162 



1 Introduction 

Magnetic edge states are formed if a confined two-dimensional electron gas 
is penetrated by a strong magnetic field. Unlike the bulk of the electronic 
eigenstates, which approach the Landau levels as the field is increased, these 
states localize at the edge of the confinement region and carry a finite current 
along the boundary. Due to their quasi one-dimensional extension and the 
ability to mediate transport the edge states play an important role in various 
phenomena of semiconductor physics, most notably the quantum Hall effect 



In the present review we shall not be concerned with the physics of interacting 
electrons in real semiconductor samples. Rather, we study an idealized system: 
A single charged particle moving ballistically in a plane which is subject to 
a homogeneous, perpendicular magnetic field. The confinement is caused by 
impenetrable walls such that the quantum wave function vanishes outside the 
considered region and the corresponding classical particle is reflected specu- 
larly at the boundary. This simple setup permits to study in some detail the 
spectral properties of magnetic edge states and their relation to the corre- 
sponding classical motion, which is typically chaotic. Thus, on the one hand, 
the present study extends the field of quantum chaos [2-7] to magnetic sys- 
tems which could not be accounted for so far. On the other hand, we expect 
that many of the results and insights obtained from the model system will 
carry over to the analysis of the more realistic case of interacting electrons in 
real samples. 

Throughout this report the confining boundary will be a closed line separating 
the plane into two parts - a compact interior and an unbounded exterior. The 
particle can move in either of these domains - forming a quantum dot or 
the respective anti-dot. In the absence of a magnetic field the interior system 
constitutes a billiard problem whose classical and quantum properties are 
a paradigm in the study of chaos and its quantum implications [8,2-7]. It 
exhibits a discrete quantum spectrum in the interior, while from the exterior 
the billiard boundary acts as an obstacle of a scattering problem. It is well 
known that there exists an intimate relation between the interior quantization 
and the exterior scattering system called the interior-exterior duality [9,10]. 
The situation changes if a finite magnetic field is present since now the exterior 



[!]• 



4 



classical motion is also bounded: The classical particle is either trapped on a 
cyclotron orbit or it performs a skipping motion around the billiard boundary. 
Consequently, the quantum spectrum is purely discrete in the exterior as well. 
It is natural to ask whether any correlations are to be expected between the 
interior and the exterior spectra, and the investigation of this issue is one of 
the motivations for the present work. 

We shall show that a duality in the underlying classical dynamics of the skip- 
ping trajectories leads to nontrivial cross-correlations between the interior and 
the exterior spectra. In order to observe this relation it is crucially important 
to have a proper quantitative definition for the notion of edge states at hand. 
Although the classical trajectories exhibit a clear partitioning into the skip- 
ping type and the cyclotron orbits, such a sharp division is no longer valid 
in the quantum treatment, and one finds many eigenstates which interpolate 
between states in the bulk and the proper edge states. In the present work we 
offer a very natural way of treating this gradual transition between the edge 
and the bulk states. It yields an objective and physically meaningful definition 
for edge states which permits a semiclassical description. 

Apart from presenting our results on the properties and the dual nature of 
the edge state spectra, the present report is aimed at providing a consistent 
and self-contained formulation and exposition of the following subjects: 

(a) the exact quantization of interior and exterior magnetic billiards based 
on boundary integral equations, 

(b) the semiclassical quantization of interior and exterior magnetic billiards 
in terms of the classical dynamics, 

(c) a spectral measure for edge states and its semiclassical form, and 

(d) the relation between the interior and the exterior edge state spectra. 

Structure of the article 

In the next two chapters, we give a survey of the classical and quantum dynam- 
ics in the free magnetic plane and in magnetic billiards, respectively. Although 
many of the statements in Chapter 2 are elementary, we shall present them 
in some detail for the sake of completeness and to introduce a consistent set 
of notations. These chapters include also the discussion of concepts, such as 
the scaling properties or the semiclassical approximation, to which we refer to 
frequently in the remainder of the report. In the first part of Chapter 3 the 
classical interior-exterior duality is explained. Turning to the quantum prob- 
lem, we introduce general boundary conditions and discuss the asymptotic 
properties of magnetic spectra. The introductory chapters conclude with the 
definition of a scaled edge magnetization. 

In Chapter 4 we solve the quantization problem in the interior and exterior of 
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arbitrary magnetic billiards by means of a boundary integral method. We ex- 
plain why spurious solutions arise initially and how they can be systematically 
avoided. The application of this method in numerical simulations, its accuracy 
and its performance is demonstrated in Chapter 5. We focus mainly on two 
issues: the computation of wave functions in the extreme semiclassical regime 
and the extraction of large sequences of eigenvalues. The former serves to vi- 
sualize the properties of edge and bulk states and the latter enables the study 
of spectral statistics and their relation to the underlying classical motion. 

Chapter 6 is devoted to the derivation of the semiclassical trace formula for 
hyperbolic and integrable magnetic billiards by means of a surface-of-section 
method. We start from the boundary integral operators and formulate the 
semiclassical quantization condition in terms of map operators which are semi- 
classically unitary and which refer to either the interior or the exterior. These 
operators are related in a way which reflects the underlying classical interior- 
exterior duality. The integrable disk billiard is then quantized for a second 
time making use of its separability. In conjunction with the former results, 
it allows the trace formula to be extended to general boundary conditions. 
This chapter is rather technical but it lays the foundation for the subsequent 
analysis. 

The spectral density of edge states is introduced in Chapter 7. It gives the 
concept of edge states a quantitative meaning and is appropriate, both in 
the deep quantum and in the semiclassical regime. As a matter of fact, we 
propose two different methods to define the edge spectral densities and discuss 
their relative advantages and connections. The new measures allow a spectral 
analysis to be performed also in the exterior. The consistency with random 
matrix theory is checked in Chapter 8 and the quantum edge state densities 
are compared to the predictions of the semiclassical trace formula. 

In Chapter 9 wc finally identify non-trivial cross-correlations between interior 
and exterior edge state spectra. We show that they are based on a classical 
duality of the periodic orbits. In order to observe the correlations the spectral 
density of edge states or an equivalent measure, such as the edge magnetiza- 
tion, is of crucial importance. We conclude this report with a summary and a 
list of open problems for further research. 

Most of the material which we considered of technical nature is deferred to 
the appendices. However, the reader may (justifiably) become impatient with 
some of the derivations deemed by us to be needed for the coherent exposition. 
The busy reader is encouraged to skip directly to Chapters 4, 6, and 7 and 
to go back to the earlier parts whenever needed. Note that a list of the most 
frequently used symbols can be found in the appendix. 
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2 Motion in the free magnetic plane 



We start by collecting a number of elementary statements on the classical and 
quantum motion in the magnetic plane. This allows to introduce the notation 
used throughout the report, and to set the stage for the discussion in the 
following chapters. In particular, the treatment of the quantum time evolution 
operator in Section 2.4 yields the opportunity to discuss the semiclassical 
approximation. In Section 2.5 the Green function of a particle in the free 
magnetic plane is derived in both its semiclassical and its exact form. 

2.1 Classical motion 

Consider the motion of a non-relativistic, spinless, charged particle in the 
two-dimensional Euclidean plane,Q which is subject to a magnetic field. Its 
Lagrangian has the form [18] 

£ = ^v2 + gvA(r), (2.1) 

where rrie and q denote mass and charge, respectively. The vectors r = [x, y)'^ 
and V = r give the position and velocity of the particle. Both of them deter- 
mine the canonical momentum 

dC 

p = — = meV + gA(r) . (2.2) 

The classical time evolution is given by the Lagrangian equation of motion 

p = gV(vA(r)) . (2.3) 

Here, the magnetic field is described by the (time- independent) two-dimensional 
vector potential A(r). It follows immediately that the equation of motion for 
the velocity v depends only on the rotation 5 = V x A of the vector potential. 
It reads 

mer = gi?V(rxv) (2.4) 



The motion on magnetic surfaces of finite curvature received some attention in re- 
cent years both in the classical [11-14] and the quantum treatment [11,15,12,16,17]. 
One motivation for introducing a non- vanishing curvature is the possibility to study 
the quantum spectrum of the free magnetic motion on a compact domain (a modular 
domain in the case of constant negative curvature). This has considerable mathe- 
matical advantages since the spectrum remains discrete in the limit of vanishing 
field. 
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which is Newton's equation of motion under the action of the (magnetic) 
Lorentz force. The latter acts perpendicularly to the velocity and is propor- 
tional to the magnetic field B (the magnetic induction). 

Throughout this report we are interested in the case of a homogeneous mag- 
netic field B (with qB > 0). Equation (2.4) is then easily integrated, yielding 
the cyclotron motion 



r{t) = r(0) + - 



v(0) (2.5) 



sin(cc;ct) 1 — cos(a;ct) 

1 + cos(a;ct) sin(a;ct) 
r(0)-p(0)+p(t) (2.5a) 



with r(0) and v(0) the initial position and velocity, respectively, and cUc—qB / rrie 
the cyclotron frequency. The particle moves clockwise on a circle with constant 
angular velocity Uc- Below, the velocity will be needed as a function of the 
initial and the final position, r(0) and r{t). Apart from the points in time 
which are multiples of the cyclotron period 2t: /uj^ it is given by 



y— sin(|a;c^) cos(|ci;ct)^ 



sin(|a;c^) 
The radius vector 



(r(t)-r(O)). (2.6) 



Pit) 1 



-Vy{t) 

Vx{t) 



(2.7) 



points from the (instantaneous) center of motion to the particle position. 

Clearly, the position of the center c(t) = r{t) — p{t) is a constant of the 
motion. To verify this in a more formal way one may consider the classical 
Hamiltonian 

H = pr-£=^(p-gA(r))^ (2.8) 

as a function of the canonically conjugate variables r and p. A short calculation 
shows that the Poisson bracket indeed vanishes, 

A(,_p).iL^ = {H,c} = 0. (2.9) 

Similarly, the (kinetic) energy E :— H(r, p) = is constant, as well as the 
cyclotron radius \p\ and the kinetic angular momentum with respect to the 
center of motion p x v, which are functions thereof. In contrast, the canonical 
momentum p itself is not a constant of the motion. In general, it does not have 
a kinetic meaning since it depends on the vector potential, cf (2.2), which is 
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not uniquely specified by the magnetic field. Rather, the gradient of any scalar 
field X(i")) ie, any "gauge field", may be added to the vector potential without 
affecting the classical equation of motion (2.4). 

We note that the general vector potential for homogenous magnetic fields may 
be written in the form 

A(r) = |(";^)+Vx(r). (2.10) 

The choice of X is a matter of convenience. An important case is the symmetric 
gauge, X = 0, which distinguishes merely a point in the plane (the origin). 
Choosing X = ~^^y-i oii the other hand, yields the Landau gauge which 
distinguishes a direction (the orientation of the y-axis). These two gauges 
are particularly important because they turn components of the canonical 
momentum into constants of the motion. In the Landau case is given by 
the (constant) y-component of the center of motion, 

A = ALan = b(^^^ ^ p^ = -TUeUJc Cy , (2.11) 

while the symmetric gauge fixes the (canonical) angular momentum with re- 
spect to the origin, L = r x p, 

A = A,,^^|(-j) ^ L:=rxp = ^(|cr-|pr) . (2.12) 

It is determined by the distance |c| of the center of motion from the origin, 
and the cyclotron radius p = \p\. 

Below, it will be important at several points to state equations in a manifestly 
gauge invariant fashion. This is done by keeping X unspecified and verifying 
that the resulting expressions do not depend on its choice. As the only restric- 
tion, X will be assumed to be a harmonic function, ie, V^X = 0) throughout. 
This rules out conveniently the occurrence of singularities in X but keeps the 
essential gauge freedom. Moreover, it ensures that the vector potential (2.10) 
is a transverse field, ie, divergence free, VA = 0, which facilitates a number 
of mathematical transformations. 

Turning to the quantum mechanical description, the quantum time evolution 
will be treated in terms of the path integral formulation in Section 2.4. Before 

that we discuss the stationary solutions of the Schrodinger equation (in a 
specific gauge, to prove the rule stated above). This permits to obtain the 
spectrum and the scaling properties of the Hamiltonian in a straightforward 
manner. 
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2.2 Quantization 



In the quantum description the canonical variables r and p become observables 
expressed as operators in £2(I^^)- They turn the Hamiltonian (2.8) into an 
operator whose spectrum determines the energies E of the stationary states. 
In position representation, p = — i/iV, the stationary Schrodinger equation 
reads 

'-\nV-qKfi){v)=E^{r). (2.13) 



Imp 



In addition, the solution '?/'(r) must be normalizable to qualify as a stationary 
quantum state. The energy eigenstates in the magnetic plane were obtained 
not before 1930, when Landau published his article on orbital diamagnetism 
[19]. Although he used the gauge (2.11), the symmetric vector potential (2.12) 
will prove more convenient in the following. First, we introduce a (quantum) 
length scale 



b--= -T^] (2.14) 




and call it the magnetic length, although it differs from Landau's definition[^ 
by a factor of y/2. It allows to transform position and momentum operators 
into dimensionless quantities, denoted by a tilde, 

_ r _ 6 

r := - and p := — p . (2.15) 
n 

In the symmetric gauge the Hamiltonian (2.8) now assumes a particularly 
simple form, 

H = ^ ^ (p2 + f ^) - h{i X p) = Hose -ujL. (2.16) 

It is given by the energy of a two-dimensional harmonic oscillator Hose minus 
its angular momentum L = r x p, in quanta of the same size. The oscillator 
eigen-frequency differs from the cyclotron frequency by a factor of 2. It is given 
by 

Imp 2 



^ Landau's definition of the magnetic length = ^/\/2 is appropriate for the 
Landau gauge (2.11). The length h (which is the suitable scale of the symmetric 
gauge) proves more convenient since it avoids the appearance of the factor 2 and 
\/2 at various places. It gives the radius of a disk the area ft^vr of which assumes the 
role of Planck's quantum, cf (3.11a). (The flux through the disk equates the "flux 
quantum" $o = h/q = Bb'^n). 
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and known from the precession of magnetic moments as the Larmor frequency. 
In order to construct the complete set of energy eigenstates on the magnetic 
plane consider the annihilation operators of the left- and right-circular quanta, 

^(l) = T + ^^P^ T ^^y)) ' (2-^^) 

with [ol, Sl] = [or, Sf^] = 1 as the only non-vanishing commutators. It is well 
known [20] that the simultaneous eigenstates of the left- and right-circular 
number operators {a\^ai^) and (0-1^0^) form a complete basis set of C2{M.'^). 
An oscillator eigenstate corresponding to n left-circular and m right-circular 
quanta is given by 

Km) = ^^(a[r(4nO,0) , (2.19) 
ynl ml 

with n,m & Nq. Here, |0,0) denotes the harmonic oscillator ground state, a 
Gaussian in position representation, (r|0,0) = exp(— ir^/6^)/\/6%. Like all 
the states (2.19) it is square-integrable and normalized. 

Inverting equations (2.18) the Hamiltonian of a particle in the magnetic plane 
may be expressed in terms of the circular operators. It assumes a form 

H = Hose — uo L = hhj (apj^o-R -|- SlSl + 1) — Tujo (dj^aR — o-LdL) 

= hw, [aiai, + ]^ , (2.20) 

which depends only on the number operator of the left-circular quanta. It 
follows that the states (2.19) form a complete set of eigenstates of the magnetic 
plane. Their energies are determined by the number n of left-circular quanta, 
called the Landau level, 

E = nw^(^ + ]^ . (2.21) 

This proves that the spectrum of H is discrete and equidistant .0 The fact that 
the energy does not depend on m shows that each Landau level is infinitely 
degenerate (with a countable infinity). This degeneracy is due to the energy 
independence of the position of the center of motion. To show that the latter 
is indeed determined by the right-circular quanta alone we note the operators 
corresponding to the classical radius vector (2.7) and the center of motion 
c = r — p, respectively. 



and c = - = - 
t J 6 2 



Sr + Or 

i(dR - 



(2.22) 



^ For mathematical literature on the spectral properties of magnetic Schrodinger 
operators see [21,22]. 
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Here, (2.2) was used to express the velocity in terms of momentum and posi- 
tion. Clearly, c commutes with the Hamiltonian like in the classical case. The 
components Px and p^, on the other hand, are not constants of the motion, 
although the cyclotron radius |p| is again fixed and determined solely by the 
energy. This can be seen from the squared moduli of the vectors, 

|p|2 = a[«L + ^ and \c\^ = a^i^ + ^ , (2.23) 

which contain only the number operators of left- and right-circular quanta. 
Consequently, the states (2.19) with fixed n and m arc eigenstates of these 
operators. They are characterized by definite expectation values for the cy- 
clotron radius and for the distance from the origin to the center of motion. 
Moreover, these stationary states are eigenvectors of the (canonical) angular 
momentum given by the difference |cp — |pp = in analogy to the classical 
result (2.12). The general eigenstate of |pp (with eigenvalue n + \) is given by 
a superposition of states (2.19) with different quantum numbers m. We will 
call any such stationary state a Landau state in the Landau level n. 



Coherent states 

Since the states (2. 19) are eigenstates of the radial components of the operators 
p and c their azimuthal components are maximally uncertain. It is known from 
the two-dimensional harmonic oscillator that the common eigenvectors of ol 
and Or have the property to minimize the uncertainty product [20]. These 
coherent states are given by the superposition 



exp 



> — |n, m) 



(2.24) 



with ttL, Or G C the associated eigenvalues. If considered in the magnetic 
plane, the expectation values of p and c are determined directly by these 
eigenvalues. 



Q!l;Q!r 



a;L;Q;R 



Q!l;q!r 



Q;L;a;R 




Im(Q;R) / 



(2.25) 



as can be found immediately from equation (2.22). The corresponding uncer- 
tainties Apa; = Apy = Acx = Acy = 6/2 are minimal, indeed. Furthermore, 
the wave functions (2.24) remain of the coherent type as they evolve in time. 
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Prom (2.20) one observes that the state at time t, 



-iHt/h 



-iujct/2 



-iojct 



(2.26) 



is merely characterized by a different phase of ckl. It is a localized wave packet 
rotating with cyclotron frequency ujc around the constant center of motion c. 
As such it embodies the closest quantum analogy [23] to the classical motion 
discussed in Section 2.1. 



Gauge invariance 

So far, the quantum problem was discussed for the symmetric gauge (2.12) 
only. We will now admit an arbitrary gauge again and consider the conse- 
quences of a finite choice of X- Although the canonical momentum is gauge 
dependent, its representation as a differential operator, p = — i/iV, contains 
no dependence on the vector potential. This can be understood by the obser- 
vation that the velocity operator 

undergoes a unitary transformation as one changes the gauge: 

— f - iW - qA(r)) = e'^^^"-)/'' — f - iW - qA,yM) e-'^^^"-)/'^ (2.28) 
me ^ ^ rue^ ' ' 

Consequently, in order to preserve the gauge independence of the velocity 
expectation value also the wave functions must be transformed unitarily as 
the gauge is changed. This is found immediately by applying (2.28) twice to 
the time dependent Schrodinger equation at arbitrary gauge, 

^ e'^^W/^ -^(-iW - qKyraf c-^^^W/'^ |^^) . (2.29) 
2jme 



Comparing the wave function with the one of the symmetric gauge, 

i/i9t|V'o) = ^(-i^^ - ?Asym)' iV'o) , (2.30) 
we see that they are related by a local, unitary transformation 

= e^^^W/'^lV^o) = e'^(^)|z^o) (2.31) 

which is determined by the gauge field \ (in dimensionless units '■ = 

2x(r)/(-B6^)). It follows that the velocity expectation value is gauge invari- 
ant. The same holds for all observables which commute with r, due to the local 
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nature of the transformation (2.31). As an immediate consequence, the prob- 
abihty density [■(/'^(r) and the probabihty flux, j(r) are also gauge-invariant. 
The latter may be identified from the continuity equation Vj = —dt\ip\^., which 
follows from (2.29), as 

j := Re(V'*vV') = — Im(V'*VV') - ■ (2.32) 

me me 

Like all observablcs which include the gradient in position representation it 
contains the vector potential explicitely to account for the gauge-dependent 
phase of the wave function. 



2.3 The scaling property 



The magnetic Schrodinger operator conventionally contains the four param- 
eters h, rrie, q, B, along with the energy E as the spectral variable. Due to 
the homogeneity of the vector potential (2.12) it is possible to reduce those 
to the two principal length scales which we encountered in the previous sec- 
tions. Those are the cyclotron radius p (2.7) and the magnetic length b (2.14), 
respectively, given by 

9 2meE , ,9 2/i , , 

P'--=^ 6^:=-. (2.33) 

The cyclotron radius is a quantity of classical mechanics. The magnetic length, 
in contrast, has a pure quantum meaning. As discussed above, it determines 
the mean extension of a minimum uncertainty state, and vanishes as — > 0. 

In the preceding section the dimensionless variables r = r/b and p = bp/h were 
introduced. In fact, the homogeneity of the potential (2.12), in conjunction 
with the requirement [x,px] = [y,Py] = i, leads necessarily to the magnetic 
length as the appropriate scale. The only freedom is a numerical factor in the 
definition of b. We took it such that the induced time scale i = ujt is given by 
the (classical) Larmor frequency u (2.17). It is appropriate to measure time 
in terms of the Larmor period T = 27r/a;, rather than the cyclotron period 
3cyc = \T, because the former is the fundamental time scale of the quantum 
problem: It takes two cyclotron periods, as one observes from equation (2.26) 
(and more generally from the propagator (2.47)), before a wave packet returns 
to its initial state with correct parity. 

The respective dimensionless Lagrangian, furnished with a tilde like all scaled 
units, reads 

~ C 1 ~ 1 ~ 

— = -v2 + fxv + v Vf X = o + vA(f ) . (2.34) 
niu 2 2 
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It contains no parameters any more, but for the definition of the scaled gauge 
field, 



X(f) := -^Xibi) (2.35) 

(which is not necessarily homogeneous of order two). This implies the defi- 
nition of the general scaled vector potential A(f) = 2A{br)/{Bb). The scaled 
Hamiltonian, given by 

shows that the proper, scaled energy reads E = E/{huj) = 2p^/6^. We will 
state the energy in terms of the spacing between Landau levels. 



E E p 



2 



(2.37) 



and call v — E/2 the scaled energy, nonetheless. This way we conform with 
the popular convention that the Landau levels start at one half, rather than 
at one. 

Below, it will be important to distinguish between the two independent short- 
wave limits of magnetic dynamics. Prom expression (2.37) one observes that 
the spectral variable v can be increased by either increasing p at constant 
magnetic length 6, or by decreasing h at fixed cyclotron radius p. The former 
direction is realized by raising the conventional energy at constant magnetic 
field. It is the standard high-energy limit. Here, the curvature of the classical 
trajectory tends to zero, which shows that in this limit the dynamical effect 
of the magnetic field vanishes. On the other hand, one may increase both the 
conventional energy and the field at a fixed ratio of E/B"^, thereby keeping 
the cyclotron radius fixed. This way the underlying classical phase space is 
kept invariant, while the magnetic length tends to zero. It is a realization of 
the semiclassical limit since 6^ plays the role of h as the semiclassically small 
parameter. 

In order to be able to consider both limits most equations will not be written 
in scaled variables, since they might depend on the choice of the independent 
variable. Rather, the formulas will be stated in terms of combinations like r/b 
so that they can be immediately replaced by scaled variables. This includes 
the scaled gradient, Vf = 6Vr, written as 

:= 6Vr , (2.38) 

which is an admittedly unusual but consistent notation. The spectral variable 
is always stated as u. 
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2-4 The free quantum propagator 



We return to the Lagrangian formulation of mechanics in order to calculate the 
time evolution operator \J{t; 0) := exp[—iHt/h] for arbitrary gauge. According 
to Feynman its position representation (for f > 0) is given by the path integral 
[24-26] 

q(t)=r 

U(t, r; 0, ro) = / P[q] exp (^W^[q]) ■ (2-39) 

q(0)=ro 

Here, the functional W attributes a classical action 

W[q\:^ /£(q(0,q(i'))dt' (2.40) 

Jo 

to all paths q(t') going from Tq to r in the given time t. (All equations are 
stated for a time independent Lagrangian, and the zero indicating the initial 
time will be omitted in the following.) 

The formulation in terms of a path integral permits the calculation of the time 
evolution operator in a straightforward manner. Its most important advantage 
is that the semiclassical approximant of the propagator can be obtained in a 
transparent way. The situation is called semiclassical if h is small compared to 
the actions (2.40). In this case the dominant contributions to the path integral 
are represented by those paths for which the phase in (2.39) is stationary. 
They are solutions of the variational problem [q] = with fixed initial and 
final position and time. According to Hamilton's principle these are classical 
trajectories. The integral is then evaluated by expanding the variations of 
(2.40) to second order. Provided the trajectories are isolated one obtains the 
asymptotic expression of the propagator to leading order in h [27]. 



det (- ^'^''■l*! 



drdrn 



exp (^W^[qci] - if^qei) (1 + 0(h)) 

(2.41) 



It is a sum over all classical trajectories qd going from Tq to r, in the given 
time t. The only quantum ingredient is the finite size of h, which sets the scale 
of the associated classical action in the phase factor. The additional phase 
shift is determined by the number i/q^j of negative eigenvalues of the matrix 
— 9^1^[qci]/(9r9ro)) [27]. The latter has a dynamical meaning [25,26], it is 
the inverse of the Jacobi field of qd, which describes the linearized deviation of 
classical trajectories with different initial momenta. The points on qd where 
classical trajectories coalesce are called focal or conjugate. They determine 
i/q^j geometrically by virtue of the Morse theorem [28]: The value of u^^^ is 
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equal to the number of conjugate points the particle encounters on its journey 
(counted with their multiplicities [28]) and is called the Morse index. 



We are now in a position to derive the time evolution operator in the free 
magnetic plane. The quadratic dependence of the Lagrangian (2.1) on position 
and velocity renders the expression (2.41) for the time evolution operator exact 
rather than asymptotic. First, the (scaled) action of a trajectory is needed as 
a function of the initial and the final position, ro and r, and the time of flight 
t. Prom the classical solution in Section 2.1 one obtains 

^iy(i,r;ro) =iy(i,f;fo) = /J(^v(?)' + v(?)A(f (?))) df 

= I [\\i')di' + [ A(f')df'+ /A(f')df' 

2 Jo JdQ J 

ro^f 

= ^(f - fo)' cot(t) - f X fo + X(i-) - Xih) . (2.42) 
Here, the action integral was split into three parts: 

[y{i',r-,ro)di'^^^^l (2.43) 

Jo sm (t) / 

/A(fOdf' = -<i^i+<i-£2)!cot(t) (2.44) 

JdQ sm (t) 2 2 

I A(f') df' = -f X fo + X(f) - X(i-o) (2.45) 

In the first, the modulus of the velocity is constant. Its value (2.43) follows from 
(2.6). The second part was made a closed line integral, encircling a domain Q, 
which is confined by the trajectory and the straight line from f back to fo- By 
Stokes' theorem one obtains (2.44), with the minus sign due to the clockwise, 
ie, negative sense of integration. The remaining part (2.45) is a line integral 
along the straight path from Tq to f . Unlike the other contributions, it depends 
on f and fo individually and carries the gauge dependence. 

In principle, more than one classical trajectory could connect the two points f 
and fo in a given time. However, since the determinant of the matrix in (2.6) 
is non-zero for i ^ mr, n = 1,2,..., the initial velocity is uniquely specified 
for those times. At integer multiples of the cyclotron period, in contrast, any 
trajectory returns to its starting point. Excluding these instances for the time 
being, the time evolution operator is determined by only one trajectory. For 
the matrix of second derivatives one obtains 

det = -^T-r . (2.46) 

\drdroJ sm^{t) ^ ' 

The determinant of its inverse has doubly degenerate zeros dXl—WK. Hence, 
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the Morse index reads z/q^j = 2[t/7r] (with [•] the integer part), and one arrives 
immediately at the time evolution operator in the free magnetic plane 



U(t,r;ro) 



1 



27ri62 sin(a;i) 
X exp 



exp 



262 



cot{ujt) — i 



.r X ro 



62 



(2.47) 



As noted above, this expression is identical to the exact path integral [29,24,30,31]. 
It is valid except for the times equal to integer multiples of the cyclotron pe- 
riod. At these instances the propagator is just a unit operator. 



1 (-)" 

hm U(t, r; Fq) = hm — — ^— exp 
U!t-*nn £-»0 27rio^ sm(£j 



262 



■ cotie 



.r X To /r\ /ro 
'r-ro 



(2.48) 



with a sign which is positive only after even multiples of the cyclotron period. 
This means that any wave function which is propagated by multiples of the 
Larmor period T — 27r/u} — 2Tf.yc returns precisely to its initial state. Equation 
(2.48) follows from a special representation of the two-dimensional 5-function 
which is given in the appendix, see (A. 33). Note that the propagator (2.47) was 
derived for positive times t > only. It is valid for all times, nonetheless, since 
it clearly obeys the unitarity relation U(— r; Fq) = [U(i, Fq; r)]*. Furthermore, 
it is given for arbitrary vector potentials. The dependence on X shows how 
the propagator transforms as the gauge is changed. It is consistent with the 
gauge dependence of the wave functions (2.31) discussed in Section 2.2. 



2. 5 The free Green function 



We are now in a position to calculate the Green function of the free magnetic 
plane. It will be an important ingredient in the theory of the exact and semi- 
classical quantization of magnetic billiards. We define the Green function to 
be the Fourier transform of the free propagator 

G(E, r; Fq) := lim / U(t, f; Fo)e^(^+'")*/'^dt . (2.49) 

2ime £-10 Jo 

As such, it is a resolvent of the Hamiltonian obeying the inhomogeneous 
Schrodinger equation 

(H - E)G{E, f; Fq) = -7^5(r - Fq) . (2.50) 
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For later reference we note that there exists a second, independent solution 
of (2.50) which differs appreciably from G. We shall call it the unphysical or 
irregular Green function G^'"-*. 

One procedure to obtain the Green function is based on the observation that 
the differential equation (2.50) separates in polar coordinates if the symmetric 
gauge is used. This way one is led to an angular momentum decomposition of 
G, which is of little use for our purposes. It was derived (with some errors) 
in [32,33] and is summarized in Appendix A. 2. Here, we perform the Fourier 
integral (2.49) directly. It yields the Green function in a clear-cut fashion, in 
Cartesian representation and arbitrary gauge. Substituting scaled variables 
the integral (2.49) reads 



G^(r;ro) := G(2^i/,r;ro 
-1 dt 



f 

Jo 



Att Jo sin(t) 



exp 



cot(i) - f X fo + X - Xo + 2i^t 



(2.51) 



with the abbreviations X '■= X(r),Xo '■— X(i'o)- (The energy u is assumed to 
have an infinitesimally small positive imaginary part.) 

Like in the case of the propagator, stating the Green function as an integral 
has the advantage that its semiclassical approximation can be obtained in a 
straightforward way. This is shown in the following. The exact integration will 
be carried out afterwards. 



2. 5. 1 The semiclassical Green function 

The semiclassical approximation to the Green function G^f^^ is obtained by per- 
forming the Fourier transform in the stationary phase approximation, which 
is summarized in Appendix A. 4. It yields an asymptotic expansion to leading 
order in the semiclassically large parameter Requiring the integrand of 

the Fourier integral (2.51) to have a stationary phase leads to a condition 

which selects the times of flight of classical trajectories connecting the initial 
position ro with the final point r at fixed energy i/. It can be satisfied only if 
the distance between the two points is smaller than the cyclotron diameter. 
If this is the case, the time derivative of the phase in (2.51) vanishes at an 
infinite number of (discrete) times, 

ig^^ = arcsin(C) -|- nn , 

^l"^ —77 — arcsin(C) -|- nn , with n — 0,1, . . . (2.53) 
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Here, 



C := ^ (2.64) 

measures the distance between the initial and the final point relative to the 
classical cyclotron diameter. The two times of flight fg^-* and i^^ belong to the 
two distinct trajectories which connect the initial and the final point directly. 
They are "short" and "long" arcs, respectively, ie, span an angle smaller and 
larger than ir (cf Fig. 6.1). For n > the trajectories perform in addition n 
complete cyclotron orbits. After the Fourier transform the trajectories entering 
the semiclassical Green function exhibit an action S — W + 2ui which is a 
function of energy u — E/2 rather than time. As specified by (2.53) the actions 
read 

^^1^^ = 27TU (a^s) + n) + X - Xo • (2-55) 
Here, we introduced the notation 

as(r; Fq) := ^ ^arcsin(C) + C^l - - 

aL(r; ro) ^ (^r - arcsin(C) - (^1 - C - (2-56) 

for the geometric part of the action. Note that as and aL depend on the 
initial and the final point individually, due to the term r x Fq, which means 
that they are not translationally invariant. However, one observes the relation 
as(r; ro) + aL(ro; r) = 1. It follows that the (scaled) action of a closed cyclotron 
orbit - a short arc followed by a long one - is given by 2^1/. 

To compute the stationary phase approximation (A. 29) we also need the sec- 
ond derivative of the phase in (2.51). It is given by (f — fg)^ cos(t) / sin'^(f) and 
at times (2.53) assumes the values ^Auy/l — C^/C (where the positive sign 
stands for trajectories of the short type). It follows that in the semiclassical 
approximation an infinite number of trajectories contributes to the Fourier 
integral. 

Gi^^Hr; ro) = ^ E (-)" (^7^) ' { (2^iKas + n) + ix - % + i^) 

+ exp (|27rii/(aL + n)+ix- ixo - i^) | 

(2.57) 

The sum over the repetitive cyclotron orbits n converges since u was assumed 
to have a small positive imaginary part. It adds a factor (1 + e^'^''')^^ which is 
singular at the energies of the Landau levels. The semiclassical Green function 
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is therefore given by a sum of two contributions, belonging to the short and 
the long arc trajectory — the principal classical trajectories connecting Tq and 
r: 



Gl^^nr;ro) 



1 



2(1 + e 



(2-)^ (CVT^f ^ ^ 

(2.58) 



This form will be used in Chapter 6 for periodic orbit theory. Alternatively, 
one can combine the short and long arc contributions pulling out that part of 
the phase which was time independent in (2.51). This leads to the expression 



G(^^)(r;ro)=exp 



-1 I X + Xo 



pO(sc) 



(r - rpj 

62 



(2.59) 



with 



■1 (27r)5 




(2.60) 

It shows that the semiclassical Green function is given by a phase factor which 
contains the gauge dependence and a real function G^^^'^-' which depends only 
on the distance between the initial and the final point. The exact Green func- 
tion has the same property, as manifest in (2.51). 

Note that the expressions (2.58) and (2.59) are defined only for separations 
smaller than the cyclotron diameter |r — ro| < 2p. For larger distances, the 
semiclassical Green function vanishes by definition, since the stationary phase 
condition (2.52) has no solution. As the distance between the initial and the 
final points approaches the cyclotron diameter, the short and long arcs coalesce 
and are therefore no longer isolated. In this case the approximation (A. 29) 
fails, which is indicated by the diverging pref actor of G^^'^^ as C ^ 1. If a 
semiclassical expression is needed for the domain |r — ro| ^ 2p, eg to describe 
tunneling effects, uniform approximations [34] must be employed as discussed 
in Appendix A. 6. 3. 



2.5.2 The exact Green function 

When evaluating the exact Green function we may separate the part of the 
phase in (2.51) which is not explicitely time dependent, like in the semiclassical 
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case. 



C(r; ro) = exp 



r X Fq 



X + Xo 



(or,. 



62 



(2.61) 



Now, the integral can be performed exactly by contour integration [35] 



G'(z) = — 



1 r"^ dt 



47r Jo sin(t) 
1 

47r 



^ exp 



i ( -cot(t) + 2z/t 



(2.62) 



Here, W^^ is the (real valued) irregular Whittaker function [36; eq (13.1.34)]. 
This expression was also obtained [37] using the separability of (2.50) in the 
symmetric gauge. 

Both, the function (2.62) and its semiclassical approximant (2.60) exhibit sim- 
ple poles as the energy u approaches the Landau levels. It is often convenient 
to remove these poles by considering the regularized version of G^, 



Gl{z) := lim cos(7r/i) G^z) 



/J,—*!/ 



(2.63) 



We finally state the regularized Green function in terms of the irregular conflu- 
ent hypergeometric function U [36] which is more common than the Whittaker 
function: 



IT 



An r(i/+|) 



■-/^m-u,i-z) 



(2.64) 



2. 5. 3 Properties of the free Green function 

Figure 2.1 displays the gauge- independent, regularized part of the exact and 
semiclassical Green function. As one expects, the exact Green function decays 
exponentially once the points are separated by a distance, |r — ro| > 2p, (ie 
z > 4z/) which cannot be traversed classically.0 At small distances, r ^ fq, 
it displays a logarithmic singularity, similar the (complex valued) field-free 
Green function [38]. We find 



^0 _ cosjiTu) 



\og{z) + + u)- 2^(1)) - ^J^^^Spa. + O(zlogz) (2.65) 



The above mentioned irregular solution of (4.2) grows exponentially beyond the 
classically allowed region. Its derivation is sketched in Appendix A.l. 



22 




01 2345678 



Fig. 2.1. Regularized gauge-independent part of the free Green function. Top: Ex- 
act (solid line) and semiclassical (dashed line) functions at = 10.1. Bottom: 
Error of the semiclassical approximation. Even at this moderate value of strong 
deviations occur only at the classical turning point -^z = ~ 6.36 and at small 
distances. (The deviations arise since the semiclassical approximation does not ac- 
count for the logarithmic singularity at z = and the tunneling into distances larger 
than the cyclotron diameter.) 

as z 0, with "^{z) the Digamma function [36]. Our method to evaluate the 
free Green function numerically with high precision and efficiency is discussed 
in [35]. 

The gauge invariant part of the Green function has the remarkable feature that 
its derivatives can be expressed by the function itself, at a different energy. 
For the regularized version one finds 

^^G° (;^) = - (i - u) (g° + GU) - '-Gl (2.66) 

=(| - ^)(| - ^) (G° + 2G^, + GU) 

+ z{\-u){Gl + GU) + '^Gl. (2.67) 

These formulas were obtained by employing the differential properties of the 
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confluent hypergeometric function [36]. Their asymptotic behavior reads 



d g ^ cosM 



1 - ^ 1/ ( log(^) + - i/) - 2*(1) - l)] + 0{z^ log z) , 

(2.68) 



z^-^Gl{z)^ ^ + 0(2;log^), asz 



0. 



(2.69) 



It can be deduced from the logarithmic representation of U in terms of the 
regular Kummer function [36; eq. (13.6.1)] and will be needed below. 



3 Introducing a boundeiry 



The motion in the magnetic plane turns into a non-trivial problem once the 
particle is restricted to a bounded domain. 



3.1 Motion in a restricted domain 



Let us assume that the particle is confined to move in a compact and sim- 
ply connected domain P C with smooth boundary T = dT>. The classical 
equation of motion (2.4) apphes in the interior of the domain V. Here, the 
particle moves on arcs of constant curvature, which may at some point im- 
pinge on the boundary. At these instances the trajectories must obey the law 
of specular reflection to qualify as a classical solution. This follows directly 
from Hamilton's principle, as will be shown in Sect. 6.3.1. Clearly, any tra- 
jectory which was reflected once must run into the boundary again. It follows 
that the phase space is in general split up into two disjunct parts. One part 
consists of skipping orbits. Their classical motion is no longer described by 
a continuous Hamiltonian flow (but by a discrete map) and may range from 
regular (integrable) to completely chaotic (hyperbolic). We will briefly review 
this classical billiard problem below, in Sect. 3.2. The remaining part of phase 
space describes the trivial motion on closed cyclotron orbits. It has a flnite vol- 
ume whenever the cyclotron radius is small enough to enable a disk of radius 
p to fit into the domain. We will call the magnetic field strong, accordingly, if 
the cyclotron radius is comparable to or smaller than the size of the billiard - 
a criterion which is purely classical. 

In the corresponding quantum problem the eigenfunctions are required to sat- 

o 

isfy the Schrodinger equation in the open domain I?, together with a boundary 
condition on the border line T (as discussed in Sect. 3.3). One observes that, 
at strong fields, the spectrum reflects the partitioning of the classical phase 
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space. There are eigenstates which hardly touch the boundary and have en- 
ergies very close to the Landau levels. They are called bulk states because 
in the limit of strong fields they constitute the major part of the spectrum. 
We will see that these states are based on that part of phase space which 
is given by the unperturbed cyclotron motion. At the same time, one finds 
eigenstates which arc localized at the boundary. These edge states correspond 
to the skipping trajectories and are expected to reflect the underlying billiard 
motion. Albeit being an effect of the boundary they may be quite significant. 
For instance, they typically exhibit a directed probability flux causing a large 
magnetic moment. This way they balance the magnetic moments of the bulk 
leading to a vanishing mean magnetization, as discussed in Section 3.4. 

The separation into edge and bulk states is intuitively clear and often used. 
Early studies were concerned with the surface electron states inside met- 
als [39,40], and after the discovery of the Quantum Hall Effect [41,42] the 
notion of edge states was employed to explain this phenomenon [1,43-47]. (In 
the latter problem the Hamiltonian must include an additional impurity po- 
tential.) However, the above characterization of edge states is not precise and 
we are not aware of a general quantitative deflnition in the literature. In due 
course, we will introduce a spectral measure, which permits to quantify the 
edge character of a state [48]. Having a meaningful spectral density of edge 
states at our disposal, it will be worthwhile to consider the quantum problem 
also in the exterior. 



Motion in the exterior 

The exterior billiard problem is obtained by restricting the particle to the 
domain \ P - henceforth called the exterior domain. From the classical 
point of view there is little difference between the interior and the exterior 
dynamics. A particle impinging on the boundary from outside is reflected 
specularly and performs a skipping motion around the billiard. Like in the 
interior the skipping trajectories cover a finite volume in phase space and are 
described by a discrete billiard bounce map. Complete cyclotron orbits, on 
the other hand, now exist for any p. The corresponding phase space volume 
is unbounded because the cyclotron center may be located at an arbitrarily 
large distance from the billiard. 

The fact that a "free particle" cannot escape to infinity but is trapped on a 
cyclotron orbit is reflected by the exterior quantum spectrum. It is discrete, 
in marked contrast to the fleld-free scattering situation. The exterior quan- 
tum problem requires the stationary wave function to satisfy the Schrodinger 
equation in \ T>, again with a boundary condition on F. In addition, the 
normalization condition implies that the wave functions must vanish at infin- 
ity. In the absence of a boundary the spectrum would be given by a discrete 
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set of Landau energies, each infinitely degenerate, as shown in the precceding 
chapter. The presence of a bilhard hfts this degeneracy turning each Landau 
level into a spectral accumulation point. This means that there arc infinitely 
many discrete eigenenergies in the vicinity of each Landau energy. 

We shall address the general quantum problem in Section 3.3. There, the main 
concern will be on the boundary conditions and the average spectral behavior, 
whereas the actual quantization is performed in Chapter 4. To prepare for the 
semiclassical quantization in Chapter 6 let us first take a closer look at the 
classical problem. 



3.2 The classical billiard 



Classical magnetic billiards were first examined by Robnik and Berry [49] and 
are still the subject of active research [50-57,14]. In this section we collect 
basic results, limiting the discussion to those aspects which will be needed 
later on. 

The classical dynamics is completely specified by the size of the cyclotron 
radius p and by the shape of the billiard. Throughout this work, the bilhard 
boundary F is assumed to be smooth, so that its normals fi exist everywhere. 
We define them to point outwards (ie, into Keeping their orientation 

fixed will allow to distinguish the interior from the exterior problem. The 
boundary is parameterized by the arc length s, 

r : s e [0;J^] ^ r(s) e R% (3.1) 

such that the derivative yields the normalized tangent 

We define the local curvature 

(r(s) - r(so)) n(s) 
Kiso) := 2 hm ^ \ ' (3-3) 

(r(.)-r(.o))' 

to be positive for convex domains. The area of the domain is denoted by 
and ^ represents its circumference. 



3.2.1 The billiard bounce map 

As mentioned above, the particle's skipping motion may be described by the 
mapping of a Poincare surface of section onto itself. Like in the case of field- 
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free billiards [58,59,8,60] it is natural to use the Birkhoff coordinates {s,ps) to 
define the surface of section. They are given by the position on the boundary 
s (the curvilinear abscissa) and the (normalized) tangential component of the 
reflected velocity Ps — Vo(s) t(s) at the point of reflection. The variables s and 
Ps are canonically conjugate in the sense described below. It is worth noting, 
therefore, that ps is defined as a component of the velocity vector, rather than 
the (gauge-dependent) canonical momentum. 

A point {s,Ps) in the Birkhoff phase space describes the position of incidence, 
and the direction of the velocity after refiection (once it is agreed on whether 
to consider the interior or exterior problem). Tracking the classical trajectory 
until its first intersection with the boundary specifies the next point of re- 
fiection s' uniquely, and p'^ follows from the law of specular refiection. Since 
any reflected trajectory is included this way the complete billiard dynamics is 
described by the bounce map 

B:{s,p)^{s',p') (3.4) 

which maps the Poincare surface of section [0; =Sf] x (— 1; 1) onto itself. In order 
to see that the map generates a discrete Hamiltonian evolution, one may look 
for a generating function '^^{s,s'), which yields the (canonically) conjugate 
coordinates by differentiation, 

d^(s,s') , , d^(s,s') ,^ , 

The relation (3.5) is the discrete analogue to the case of continuous Hamil- 
tonian dynamics, where the canonical momenta are similarly given by the 
derivative of the action. If the mixed second derivative of ^ has a deflnite sign 
the equations (3.5) may be globally inverted [60], yielding the bounce map 
(3.4). 

The billiard dynamics may now be studied conveniently by investigating the 
properties of the map. In Fig. 3.1 we show surface of section plots of an 
interior ellipse at different values of the the cyclotron radius. One observes the 
standard picture of mixed chaotic dynamics [61-63]. The trajectories either 
lie on invariant curves (characterizing regular motion) or cover a whole area 
in the surface of section (chaotic motion). Stable periodic orbits, in particular, 
are characterized by surrounding invariant lines ("stability islands"). 



3.2.2 Integrable and hyperbolic billiards 

In the field free case the ellipse is known to be the only smooth and simply 
connected billiard with two integrals of motion (including the circle as a spe- 
cial case). At finite magnetic fields, the ellipse turns chaotic, as we have just 
seen, except for the circle billiard. The latter exhibits the canonical angular 
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Fig. 3.1. Birkhoff phase space portraits of the interior ehipse (strong eccentricity 0.8, 
area jz/ = tt), for different values of the cyclotron radius p = 0.40,0.44,0.50,0.54 
(left column, top to bottom) and p = 0.6, 1.0, 2.0, 10.0 (right column, top to bottom). 
The motion turns (more) regular as the limit of a strong field, p — >■ 0, and a vanishing 
field, p ^ oo, is approached, [figure quality reduced] 
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Fig. 3.2. Parts of dual trajectories in the interior and exterior of a stadium-billiard 
at /9 = 0.5 (sequence of 75 reflections). The billiard shape is defined in Fig. 8.5. 

momentum (2.12) as the second integral of the motion (provided the circle 
is centered at the origin of the symmetric gauge). This suggests that circular 
shapes, ie, the disk and the annular billiard, are the only boundaries which 
yield integrable motion in the magnetic field. 

The other extreme type of motion is called hyperbolic, or displaying hard chaos. 
It is present if the stable part of phase space has zero measure rendering al- 
most all trajectories unstable. Hyperbolic billiards are popular, although they 
form a small class. Early examples of field-free billiards displaying hard chaos 
were given by Sinai [58] and Bunimovich [59]. Conditions for the instability 
of orbits in magnetic billiards are discussed in [53-55] . In his recent work [14] 
Gutkin applied a general hyperbolicity criterion [64] to construct classes of 
hyperbolic magnetic billiards. The critical parameter in these sets is given by 
the sum of the reciprocal cyclotron radius and the (local) curvature of the 
boundary. Hard chaos is guaranteed in these cases only for cyclotron radii 
above a certain minimal value. Most of the billiards studied numerically in 
this report are hyperbolic at zero field, but assume a mixed chaotic phase 
space at any finite cyclotron radius. An example of a billiard shape which 
generates truly hyperbolic motion even at fairly strong fields is given in the 
right part of Fig. 5.1. 

Since the above statements apply equally to the interior and the exterior 
dynamics there was no need to distinguish between them. We now turn to the 
question of how the classical interior and exterior problems are related. 
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Fig. 3.3. Breakdown of the duality in segments of partially corresponding trajectories 
(stadium of Fig. 3.2 at p = 0.8.) Only the two left (top) arcs in the interior (exterior) 
meet with a dual partner. The breakdown occurs because a cyclotron orbit, which 
is obtained by continuing the arcs (dashed line) intersects the boundary more than 
twice. 

3.2.3 The classical interior- exterior duality 

When comparing the interior and the exterior motion the size of the cyclotron 
radius p plays a crucial role. An important situation is encountered if the 
cyclotron radius and the billiard shape are such that any circle with radius 
p intersects the boundary at most twice. For convex domains, a sufficient 
condition is the cyclotron radius being greater than the maximum radius of 
curvature, or less than the minimum radius of curvature. However, convexity 
is not necessary for the above condition — which we shall assume to hold for 
the moment. 

Now consider a segment of an interior trajectory going from r(s) to r(s'). The 
same two points are connected by a valid exterior trajectory which travels 
backwards in time. Necessarily, the two arcs form a complete circle of radius 
p. (They do not intersect with the boundary, except at the points r(s) and 
r(s'), because the above criterion was assumed to hold.) The interior trajectory 
is reflected specularly and finally runs into the boundary at r(s"). Clearly, the 
time-reversed exterior trajectory obeys the same law of specular refiection, 
leading to the same boundary point r(s"). It follows that the interior dynamics 
and the time-reversed exterior one are described by the same Poicare surface 
of section. Every interior trajectory is linked with a dual exterior trajectory, 
which travels backwards in time. We call this property the classical duality of 
interior and exterior motion. Pairs of dual trajectories are displayed in Figure 



30 




Fig. 3.4. Fraction of the complete interior phase space belonging to arcs which 
violate the duality criterion, as a function of the cyclotron radius p. (Calculated 
for the stadium billiard in Fig. 3.2; the squares indicate the error of the statistical 
sampling.) 

3.2 and 9.8. 

As an immediate consequence of the classical duality one finds for any given 
interior periodic orbit a dual periodic orbit in the exterior, and vice versa. 
Being periodic, both may now be thought of as running forward in time, but 
then with opposite orders in the sequence of reflection points. Clearly, these 
dual partners are intimately related. We will see that they have the same 
stability properties and that the sum of their actions is an integer multiple of 
the action of a full cyclotron orbit (with the integer given by the number of 
reflections). Examples of dual periodic orbits are given in Figure 9.8. 

Figure 3.3 shows that the duality breaks down once the duality condition that 
"any circle of radius p intersects the boundary at most twice" is no longer 
fulfilled. Typically, only a small fraction of the phase space corresponds to arcs 
which violate the duality condition. Fig 3.4 gives an impression of the fraction 
of phase space belonging to arcs whose extension intersects the boundary more 
than twice. 

3.3 Quantum billiards 

An early study of a magnetic quantum billiard was carried out by Nakamura 
and Thomas [65] (see [66] for a correction) . Later works are concerned with the 
spectral implications of the absence of time-reversal invariance [67-69]. Special 
geometries, such as the disk [70,71] or, more recently, the square [72,73], re- 
ceived attention as well. All these studies were limited to the first few hundred 
eigenvalues, and only to the interior problem. 
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3.3.1 General boundary conditions 



The mentioned works use Dirichlet boundary conditions, ic, demand tiie wave 
function to vanish on the boundary. It is the natural choice from a physical 
point of view, which takes the boundary as due to an infinite potential wall. 
However, it will prove worthwhile to consider slightly more general, "mixed" 
boundary conditions which include the Dirichlet choice as a special case. They 
are defined by the equation 



The lower sign stands for exterior problem and the symbols dn/b '■= bh{r)'Vr 
and An = fi(r)A denote the scaled normal derivative and the normal compo- 
nent of the scaled vector potential, respectively. 

The "mixing" parameter A interpolates between the two extremes, Dirichlet, 
A = 0, and Neumann boundary conditions, A~^ = 0. In principle, A may 
be a function of the position on the boundary, but will be assumed con- 
stant throughout. At non-vanishing A our boundary conditions (3.6) are the 
gauge-invariant generalization of the mixed boundary conditions known for 
the Helmholtz problem [74-76]. They imply that the normal component of 
the current density j„ = Im^ip* dn/b'ip) — A.n\'ip\'^ vanishes for any A. (Take the 
imaginary part after multiplying (3.6) with ijj*.) The resulting conservation 
of the probability density explains why the condition (3.6) keeps the problem 
self-adjoint for any A. The explicit appearance of the vector potential in (3.6) 
is needed to ensure the gaugc-invariance of the boundary conditions. The fact 
that the definition docs not depend on the gauge freedom x is easily seen ob- 
serving the gauge dependence of a general wave function (2.31). Finally, note 
that A has the dimension of a length, cancelling the dimensionality introduced 
by the normal derivative. The magnitude of the latter depends on the modulus 
k = y/2meE/h of the wave vector. To account for this trivial energy depen- 
dence of the eigenstates on the boundary condition it will be convenient (later 
in the semiclassical treatment) to use the dimensionless mixing parameter 



We did not state the definition (3.6) of the boundary condition in terms of 
A because its dependence on the spectral variable u would destroy the self- 
adjointness of the problem rendering different eigenstates non-orthogonal. 

A quite different type of boundary conditions for magnetic billiards was pro- 
posed recently by Akkermans et al. [77]. It was designed specifically to be 
sensitive on the "chirality" of the wave functions. For the special situation of 
a separable problem (disk billiard) they allow to split the interior eigenspace 
into two subspaces with definite chirality. We will see that this is quite close 




(3.6) 




(3.7) 
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to the desired separation into bulk and edge states. However, it cannot be 
generalized to billiards with arbitrary shapes, and the resulting spectrum has 
no relation to the standard Dirichlet problem. Below, we take a different ap- 
proach to separate edge and bulk, by adjusting the spectral measure according 
to our needs, rather than modifying the spectrum. 



3.3.2 The quantum spectrum 

Unlike their field-free relatives, magnetic quantum billiards offer two indepen- 
dent external parameters - the cyclotron radius and the magnetic length. As 
discussed in Section 2.3, one must specify which one is to be fixed in order 
to define a quantum spectrum. In the main part of this report the formulas 
for spectral densities are constructed at a fixed magnetic length h. This is 
done to avoid clumsy notation (and to minimize the danger of confusion). A 
summary of formulas for spectra defined in the semiclassical direction is given 
in Appendix A. 7. Still, some of the numerical investigations presented below 
are carried out on spectra defined in the semiclassical direction, which will be 
clearly indicated. 

The simplest function to characterize a spectrum is the spectral staircase (or 
number counting function) which gives the number of spectral points below 
the specified energy. For a set of eigenvalues {f„} it is formally defined as a 
sum 

oo 

n{v):^Y.®(ij-v^) (3.8) 

n=l 

over Heaviside step functions ©. Note that N(i/) is a well-defined function only 
for the interior problem, due to the infinite number of exterior bulk states close 
to each Landau level. The spectral density is conveniently defined as the energy 
derivative of the counting function, 

d °° 

d{v) — N(i.) = E ^(^ - ^n) , (3.9) 

and should be understood in the sense of distributions. Formally, such a sum of 
Dirac 5-functions could be defined for the exterior problem as well. However, 
this density would be meaningful at most in a local sense since the convolution 
with a test function would diverge at all the Landau energies. Therefore, the 
following discussion of the smooth, asymptotic properties of magnetic spectra 
must be restricted to the interior problem. 
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Fig. 3.5. At strong magnetic fields, p <C the major part of the available phase 
space consists of complete cyclotron orbits. The skipping orbits give rise to a net cur- 
rent along the boundary. It has a counter-clockwise sense of orientation, in contrast 
to the cyclotron orbits. 

3.3.3 Asymptotic counting functions 

The spectral staircase is described asymptotically by the mean number count- 
ing function N(z/), which is uniquely defined [78]. For Dirichlet boundary con- 
ditions it is given by the asymptotic expression [79] 



The expression includes only geometric quantities and the conventional wave 
vector \/2meE /h = 2a/z7/6, which are all independent of the magnetic field. 
The field independence of the leading order term follows immediately from 
Weyl's law, as discussed below. However, it is not obvious that the next two 
orders are identical to the field free case as well. This was proved only recently 
in [79], and for circular billiards in [80]. 

Note the hierarchy of the geometric quantities appearing in (3.10). The lead- 
ing and the second term are proportional to the area and the circumference, 
respectively. The constant is determined|^ by the mean curvature /pK(s)ds = 
2tt. Moreover, the higher order terms are typically proportional to higher mo- 
ments of the curvature [81]. This hierarchy reflects the systematic method to 
derive the boundary corrections to asymptotic quantities (see eg. [75]): The 
boundary is locally approximated first by a by a straight line, then a circular 
arc, and so on. 



^ The constant term in (3.10) is modified if there are corners in the boundary [79]. 




(3.10) 
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Weyl's law revisited 



Let us consider Weyl's law more explicitely. It states that the number of quan- 
tum states below a given energy is determined, to leading order, by the volume 
of phase space, within the energy shell, divided by (a power of) Planck's quan- 
tum 

'^"''''■'"''^ (2??^//®*'^""*'*''''*''' (3.11) 

This is the first term in the asymptotic expansion (3.10). Changing the in- 
tegration of the canonical momentum to the velocity vector in the first line 
renders the phase space integral independent of the magnetic field (since the 
Jacobian is constant [82]). This shows immediately that the leading order term 
of the counting function (like any quantity which may be written as a phase 
space integral of position and velocity) cannot depend on the field strength. 
In (3.11a), however, we transformed the variables of integration to the radius 
vector p', cf eq (2.7), and the cyclotron center c = r — p', which do depend on 
the magnetic field. As a result, the role of Planck's quantum is now played by 
the area b'^n. This second form of the phase space integral has the advantage 
that it permits to separate the volumes of skipping and cyclotron motion. The 
center c is a constant of the motion for all cyclotron orbits. Hence, integrating 
only the cyclotron part of the centers one obtains the area s^cyc{p) of the set of 
points in V with a distance from the boundary greater than p. Consequently, 
the number of quantum states which correspond to cyclotron motion is given, 
to leading order, by the integral 

Ott f'P 

Ncyc(p^ b') = ^ <c(p')p'V . (3.12) 
We note from (3.10) that the total number of states reads to leading order, 

N»(P^^.=) = ^^. (3.13) 

Hence, the number of states associated with the skipping part of phase space 
can be written as an integral 

Nskip(p', &') = Ntot - Neyc = [ =<kip(p')p'V , (3.14) 

involving the area £^kip(p) := — •^cyc{p)- By definition, this area is given by 
those points in the interior domain which are closer to the boundary than the 
cyclotron radius, cf Fig. 3.6. It determines the mean density of those states, 
which correspond to the skipping part of phase space. 
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Fig. 3.6. The dark shaded region indicates the area determining the phase space 
volume of interior (left) and exterior (right) skipping orbits. It is given by those 
points of the interior and exterior domain, respectively, which have a distance less 
than one cyclotron radius p to the boundary. 

dskip{i^) = :3-Nskip(z^o ,b ) = 73 (3-15) 

dz/ b^TT 

This is a remarkably simple and intuitive formula. It should be made clear, 
however, that we do not yet have a criterion at our disposal, which provides 
a clear distinction of edge and bulk states. Clearly, a reasonable definition 
should pass the requirement of being consistent with (3.15). 

Furthermore, a proper "density of edge states" will have to be well-defined also 
in the exterior. Let us therefore comment on the expected mean number of 
exterior states which correspond to skipping motion. By symmetry, it should 
be determined by the area ^^Tp those points in the exterior domain which 
are closer to the boundary than p. This can be confirmed for the circular 
geometry, where the integral over the skipping part of phase space in (3.11a) 
can be performed explicitely. For a disk of radius R one obtains 

-z/2 liJulxR 

Tf™* J 3 27ro 2 (o ^c\ 

Nskip = { ^ (3-16) 
\t7^v \i^fvh>R 

for the interior and the exterior problem, respectively. Note that the interior 
number is determined by the area =2/ of the domain once the cyclotron radius 
p = ^ b exceeds the radius R of the disk, preventing any cyclotron orbits in 
the interior. 

At very strong fields, b p <^ =Sf, in contrast, it is the circumference term 
which dominates. Since in this case we may neglect the mean curvature, the 
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average number of skipping states is approximately given by 



Nskip „ 



4 ^ 



3 



(3.18) 



3 27r6 



This expression coincides with the phase space estimate for a straight hne 
with periodic boundary conditions (see Appendix A. 6). 

Let us turn to another quantity which serves to characterize interior mag- 
netic bilhards — the orbital magnetism which measures the response of the 
spectrum to changes in the magnetic field. Its asymptotic properties may be 
related to a phase space integral as well. 



3.J^ Orbital magnetism 

Employing the notion of orbital magnetism we slightly abuse a thermody- 
namic concept for our one-particle problem. Nonetheless, it is worthwhile to 
ask for the magnetic response of the billiard dynamics in the sense of statisti- 
cal mechanics. We consider only micro-canonical ensembles (since we are not 
concerned with effects of finite temperature) which means that averages are 
performed on the energy shell in phase space, ie, among all orbits of a given 
cyclotron radius. 

Let us first consider the classical motion along a single periodic^ trajectory. 
Being charged the particle constitutes an electric current which in turn induces 
a magnetic moment. Will it serve to strengthen or to weaken the applied 
magnetic field? Clearly, the latter is expected in the case of a cyclotron orbit. 
Here, the (scaled) magnetic moment turns negative. 



which shows that the cyclotron part of phase space is diamagnetic. The skip- 
ping orbits, on the other hand, will in general give rise to both signs. At 
strong fields (if the cyclotron radius is shorter than the minimum diameter 
of the billiard) skipping trajectories carry a net current along the boundary. 
It is orientated clockwise, ie, opposite to the cyclotron orbits (see Fig. 3.5). 
A detailed analysis [83] shows that, in any case, a subtle cancellation mech- 
anism between cyclotron and skipping orbits is at work, which guarantees 
that classically there is no net orbital magnetization. This is the van Leeuwen 
theorem [83,82]. 



" We may confine the discussion to periodic orbits because the set of periodic orbits 
is known to be dense in phase space. 




(3.19) 
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The statement is proved immediately by evoking the thermodynamic definition 
of the magnetization as the derivative of a thermodynamic potential (the free 
energy or the grand canonical potential) with respect to the magnetic field. 
The potentials are determined by the partition sum, which is a phase space 
integral in the classical case. As such it cannot depend on the magnetic field 
for the reasons given in the preceding section [82]. 

Before we turn to the precise quantum definition it should be emphasized that 
orbital magnetism in its proper sense is an effect of many particles at finite 
temperature. Assuming the temperature to be much larger than the spacing 
between Landau levels, T ^ hwc/kB, Landau showed [19] that a degenerate 
Fermi gas exhibits a smalQ net diamagnetic response. This Landau diamag- 
netism is an effect of the bulk. Asymptotic corrections due to the existence of 
a boundary are discussed in [85-87,79,88,80]. Recently, the effect met some re- 
newed interest since the geometry of mesoscopic devices may greatly enhance 
orbital magnetism. Semiclassical treatments in terms of periodic orbit theory 
may be found in [89-91,72,92]. In these works the magnetic field was assumed 
to be very weak such that the bending of the trajectories could be neglected. 
An exception is the study of the quantum and semiclassical magnetization of 
the magnetic disk in [93]. A comprehensive review on the subject of orbital 
magnetism is given in [94]. 

In the following we shall use the concept of orbital magnetization merely as 
a means of characterizing magnetic billiards. We shall argue that it is advan- 
tageous to adopt a modified definition of orbital magnetization. In order to 
motivate this we start with the conventional one. 



Conventional magnetization 

Given the spectrum {En} at finite magnetic field B one may conventionally 
define the magnetization as 



M^onAE, B):=-J2 m{E'; B) dE' . (3.20) 



This is the one-particle and zero-temperature limit of the standard thermo- 
dynamic definition. By means of equation (3.20) the function m{E, B) is in- 
troduced which we call the magnetization density, 

m(E, B) ■■=^{E, B)=-Y.^5{E- E„) . (3.21) 

The relation of m(S, B) to the electrodynamic interpretation of the magneti- 
zation is seen once we note the derivative of the Hamilton operator (2.8) with 



^ The effect is one third of the Pauh spin paramagnetism [84]. 



38 



respect to the magnetic field, 



^ = -| (r X v)sy„. . (3.22) 

It is the operator of the magnetic moment, where (■)sym indicates the sym- 
metrized form. It follows that the energy derivatives dEn/dB in eq (3.21) are 
given by the corresponding expectation values of the magnetic moment, ie, 
the magnetization density (3.21) reads 

m{E,B) = J2l{^^\{r X ^r),,,^\^^) 6{E - E^) . (3.23) 



The fact that the mean magnetization (density) vanishes follows immediately 
from the field-independence of N (3.10), as noted above. At strong fields the 

negative moments of (many) bulk states are balanced, consequently, by the 
large, positive magnetic moments of relatively few edge states. This is seen 
much more clearly once we modify the definition of the magnetization such 
that is complies with the scaling properties of the system. 



Bulk and edge magnetization 

We proceed to define a scaled magnetization which has considerable advan- 
tages compared to the conventional one. According to (2.37) the spectrum 
{un} depends parametrically on the magnetic length, z/„ = Vn{h'^). It is natu- 
ral to define the scaled magnetization density such that it yields the density 
of the scaled magnetic moment (3.25), in analogy to (3.23). Hence, one is led 
to the definition 

fil(l/,62) ■.= Y.\{'^n\{iy-'^\yMK^-^n) ■ (3.24) 
n ^ 

From the explicit form of the scaled Hamiltonian one can easily show that 

^(V^„|(f X v)syn.|V'n) = - l^n ■ (3.25) 

Thus, the scaled magnetization density can be written in terms of derivatives 
of the number counting function. 
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The scaled magnetization follows by integrating the density. 



M{u, b') b') du' = (b'^ - (3.27) 



= Mcdgc + -Mbulk (3.27a) 

As indicated in the second line the scaled magnetization splits up naturally 
into two parts which we call, respectively, the edge magnetization. 



A<edge(^, b') := E = - I ^'^N(i.', b') du' , (3.28) 

and the hulk magnetization, 

>tbuik(^, b") ■■=-E^n = - r i^'^^ii^': b') du' . (3.29) 

This naming is appropriate since any Landau state (2.19) exhibits a scaled 
magnetic moment (n, m||(f x ^)syra\'nim) = —{n + |) = like the classi- 
cal cyclotron orbit (3.19). Each eigenstate contributes to both magnetization 
densities, 

medge(z^, b'^) ^ ^b"^^ - ^n) (3.30) 



and 



nibulk('^, b'^) = -^Vn^i^ - l^n) ■ (3.31) 



The energies of bulk states he close to the Landau levels and the nearer they get 
to the level the less they depend on b^ (since the Landau energy is independent 
of 6^). Hence, they give rise to a negligible edge contribution. Edge states, in 
contrast, contribute to the edge magnetization much stronger than to the bulk. 
This follows from the mean values of the magnetization. For the smooth edge 
magnetization density one finds, cf (3.10), 

_ / ,2'9N ^ 1^1 

niedge(-,M = -^'^ = ^--2^-^- (^-''^ 

Remarkably, the bulk mean value assumes a form, 

mbuik(^, b ) = b ) = --^u +2 2^^^, (3.33) 

which cancels the mean edge magnetization identically. Hence, the mean (to- 
tal) magnetization, — Aledge+A^buik vanishes like in the conventional case. 
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This holds strictly for any field, independently of whether or not there is a 
classical separation into skipping and cyclotron orbits. 

The edge magnetization (3.28) defined in this section embodies a first quan- 
tity which allows to distinguish edge states quantitatively. It gives the excess 
magnetization of the states which arises due to the existence of a boundary, 
ie, as compared to the expected diamagnetism of a state in the infinite plane. 



4 Quantization in the interior and the exterior: The boundary in- 
tegral method 

In the present chapter, we show how to solve the quantization problem for in- 
terior and exterior magnetic billiards by means of a boundary integral method. 
It provides the spectra and wave functions of arbitrarily shaped billiard do- 
mains, and includes the general boundary conditions discussed in Section 3.3.1. 
Moreover, the boundary integral formalism constitutes the basis for the semi- 
classical theory discussed in Chapter 6. 

4-1 Boundary methods 

As compared to the field free case, it is surprisingly difficult to obtain the 
quantum spectra of magnetic billiards. So far, numerical studies were restricted 
to the interior problem and performed almost exclusively by diagonalizing the 
Hamiltonian [65,67-69,95]. This requires the choice and truncation of a basis, 
which is problematic for general billiards, where no natural magnetic basis set 
exists. Consequently, results were limited to the first few hundred eigenvalues. 

In the case of field free billiards quantum spectra are usually obtained by 
transforming the eigenvalue problem into an integral equation of lower dimen- 
sion. The corresponding integral operator is defined in terms of the free Green 
function, and depends only on the boundary [96-101]. This method is known 
to be more efficient than diagonalization by an order of magnitude [102,103]. 
We proceed to extend these ideas to magnetic billiards. A step in this di- 
rection was taken by Tiago et al. [33], who essentially propose a null-field 
method[3 [lOA] for (interior) magnetic billiards. It involves the irregular Green 
function (A. 14) in angular momentum decomposition. A drawback of the ap- 
proach is that this function must be known for large angular momenta, which 
turns out to be numerically impractical. Moreover, the method does not apply 
for the exterior problem. 

^ The authors of [33] inaccurately call their scheme a "boundary integral method" . 
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In the following we derive the boundary integral method for magnetic billiards. 
Like in the field free case, it involves the regular Green function in position 
space representation. We present the method for the interior and the exterior 
problem, and general boundary conditions. 



4-2 The boundary integral equations 



4-2.1 Single and double layer equations 

The stationary eigenfunction of a magnetic billiard at energy v is defined by 
the differential equation 

(i( - iVr/b - A(r))' - ^{r) = , (4.1) 

and a specification of the wave function on the billiard boundary F. The 
free Green function, Gj,, was shown to satisfy the inhomogeneous Schrodinger 
equation 

(i(-iV./,-A(r))'-2^)G.(r;ro) = -|5(^) . (4.2) 



Our goal is to cast the quantization problem into an integral equation defined 
on the billiard boundary. To that end, we take the complex conjugate of (4.1) 
and multiply it (from the left) with G^. Similarly, equation (4.2) is multiplied 
with ip* and subtracted from the former expression. One obtains an equation 

r v,%G, - G, v.^V'* - 2i ^r/b{ArG.) = r ^{^^) > (4-3) 

which has a form suitable for the Green and Gauss integral theorems. It holds 
everywhere in the plane, except for the boundary F, where the boundary 
condition (3.6) introduces a discontinuity in the derivative of ip. 

We start by considering the interior problem and sketch the treatment of the 
exterior case afterwards. Choosing the initial point of the Green function away 
from the boundary, fq G ]R^\F, the integral of (4.3) over the (interior) domain 
V may be transformed to a line integral. 



^ [iP* {dn/bG^ - i in G^) - G^ {dn/b^* +iAnip 



dF 

T 

V'*(ro) if To e I) 
if To e \ P 



2 . ^ (4-4) 



It is defined on the boundary F (with the normal components of the vector 
potential and the gradient denoted as A„ = n(r)A and dn/b ■— 6n(r)Vr, 
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respectively). Note that the vector potential part of the integrand was split 
which is necessary for a gauge invariant formulation of the integral equations. 



The single layer equations 

We choose Fq e F and define Tq := tq ± efio, for small £ > 0. By adding the 
two equations in (4.4), one obtains 



/dr 
, r{dn/bGl-iAr,Gl)-Gl{dn/br + '^Ar,r))\-^ = W{^o) ■ (4-5) 



Here, we used the abbreviation = \Gu{y]Yq) + |Gi,(r;ro). Equation (4.5) 
holds for all (sufficiently small) e > 0, hence the limit e exists. Moreover, 
observing the asymptotic properties of the Green function (cf Sect. 2.5.3), it 
can be shown, that the integration and the limit G^ G^, dn/hG^, dn/bGu 
may be interchanged. Inserting the boundary condition (3.6) we obtain, after 
renaming the limiting function u — dn/bi^* + i^nV'*) '^Q '■— u{yq), 



L 



A 



h 



dr A 



G,--{dn/bG,-iAnG,) u—^-{-\uq), (4.6) 



h b 



an integral equation defined on the boundary T [35]. 

In order to derive the corresponding equation for the exterior problem, con- 
sider a large disk /Cp D P of radius p, and integrate (4.3) over /Cp H T>. Once 
ro lies in the vicinity of F, the contribution of dICp to the boundary integral 
vanishes asp ^ oo, due to the exponential decay of the regular Green function 
Gu- Similar to eq (4.5) one obtains an equation which permits the limit £ — > 
to be taken before performing the integration. The resulting boundary integral 
equation differs from (4.6) only by a sign. In the following, we shall treat both 
cases simultaneously, with the convention that the upper sign stands for the 
interior problem, and the lower sign for the exterior one. 



A 



dr A 



G^::f-{dn/bG^-iAnG^) u — = ji-^uo) . (4.7) 



b b 



In analogy to the Helmholtz problem [96], we will refer to these equations as 
the single layer equations for the interior and the exterior domain. 
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The double layer equations 



A second kind of boundary integral equations can be derived by applying the 
differential operator + iA^g) := n(ro)(Vro/fe + iA(ro) on equation (4.5), 



This equation is true for all e > 0, which means that the limit e ^ exists. As 
for the second integral, we may again permute the limit and the integration 
which yields a proper integral. Consequently, the limit of the first integral is 
finite, too. However, in the first integral we are not allowed to exchange the 
integration with taking the limit because the limiting integrand (4.26) has a 
l/(r — ro) ^-singularity which is not integrable (see below). 

Integral operators of this kind are named hypersingular [105]. Similar to a 
Cauchy principal value integral, they are defined by taking a special limit. 
However, compared to the principal value the singularity is stronger by one 
order in the present case. Below, in Section 4.3, we define which limit is to 
be taken. It is denoted by ^ and should be read "finite part of the integral" . 
With this concept and equation (3.6), we obtain the double layer equations. 



which are again integral equations defined on the boundary F. 
The spectral determinants 

It is useful to introduce a set of integral operators (whose labels D and N 
indicate correspondence to pure Dirichlet or Neumann conditions): 



/dl 
, -0* (<9no/b + i ^no) {dn/bGl - i A„ G^) — 

/- dl 



±i(9„„/, + ii„jV'*(rJ) ■ (4.8) 




T^uo, (4.9) 




U 



(4.10) 



(4.12) 



(4.11) 



(4.13) 
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They act in the space of square-integrable periodic functions, u e J^2(r), with 
the period given by the circumference 



Nontrivial solutions of the single layer equations (4.7) and double layer equa- 
tions (4.9) exist for energies where the corresponding Predholm determinants 
vanish, 



det 
det 



QdiTAQ^i±^id 



= (single layer) 
= (double layer). 



(4.14) 
(4.15) 



Hence, these are secular equations although the explicit dependence on the 
spectral variable is not shown in our abbreviated notation. However, each of 
the determinants (4.14) and (4.15) may have roots, which do not correspond 
to solutions of the original eigenvalue problem given by (4.1) and (3.6). For 
finite £, the equations (4.5) and (4.8) are still equivalent to the latter. They 
acquire additional spurious solutions only as they are transformed to boundary 
integral equations by the limit £ ^ 0. 



4-2.2 Spurious solutions and the combined operator 

The physical origin of the redundant zeros is apparent in our gauge invariant 
formulation: They are proper solutions for the domain complementary to the 
one considered. This is obvious for the single layer equation with Dirichlet 
boundary conditions (A = 0), where the spectral determinant docs not de- 
pend on the orientation of the normals. The same is true for the double layer 
equation with Neumann boundary conditions (A"^ = 0). 

In general, the character of the spurious solutions may be summarized as 
follows: Independently of the boundary conditions, the single layer equation 
includes the Dirichlet solutions of that domain which is complementary to 
the one considered. Likewise, the double layer equation is polluted by the 
Neumann solutions of the complementary domain, irrespective of the bound- 
ary conditions employed. This statement is easily proved by observing that 
the single-layer-Neumann operator and the double-layer-Dirichlet operator are 
adjoint to each other, = (Qdi)^ while the operators and Q[5 arc self- 
adjoint (see below). Now assume that m is a complementary Dirichlet solution. 
In Dirac notation, 

Q°|«) = A Q^^\u)^l\u)^0 (4.16) 
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Applying the dual of u to the single layer operator yields 

(«|Q°TA{(«|Q,^TiM} = 0, (4.17) 

which implies that the Predholm determinant of the single layer operator van- 
ishes. Similarly, if is a complementary Neumann solution, 

±Q^|m) + ||m)=0 a Qdil^^) = (4.18) 
±(M|Q° + i(M| =0 A (m|Q^i = 

then its dual satisfies the double layer equation, again for any A, 

±{ ± («|QS + |(«|} T X{u\Ql = • (4.19) 

Since the spurious solutions are never of the same type, it is possible to dispose 
of them by requiring that both, the single and the double layer equations, 

should be satisfied by the same solution u. Therefore, one obtains a necessary 
and sufficient condition for the definition of the spectrum by considering a 
combined operator 

Qc^ := (qS T AQS ± ^id) + iae T XQ^i + ^id) , (4.20) 

with an arbitrary constant etc- It has a zero eigenvalue only if both, single and 
double layer operators do. In practice, the spectrum is obtained by finding the 
roots of the spectral function 



C{iy)^det . (4.21) 



It is worthwhile noting that (for the interior problem) spurious solutions will 
not appear if one uses the irregular Green function. The reason is that the 
gauge-independent part of this function is complex, which destroys the mu- 
tual adjointness of the operators. This is why the irregular Green function 
had to be chosen for the null- field method [33]. For the boundary integral 
method, the option to use this exponentially divergent solution of (4.2) is ex- 
cluded, since the corresponding operator would get arbitrarily ill-conditioned 
once the diameter of the domain V exceeds the cyclotron diameter. The exte- 
rior problem cannot even formally be solved using G^'''"'-* (due to an essential 
singularity at the origin). 

A last remark is concerned with the important case of Dirichlet boundary 
conditions. Here, one could as well derive a pair of boundary integral equations 
that are not gauge- invariant. (Just set ip* = in (4.4) and consider u = 
dn/b'4'*-) Of course, these equations would yield all the proper gauge-invariant 
eigen-energies of the problem. However, the energies of the additional spurious 
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solutions would depend on the chosen gauge, and a characterization of the 
latter in terms of solutions of a complementary problem would not be possible. 

The fact that the spurious solutions can be removed by considering a combined 
integral operator is of great practical importance for numerical calculations 
[35,106]. An individual spurious solution of the single or the double layer 
operator may be identified as well after evaluating the corresponding wave 
functions by observing in which domain it vanishes. 



4.2.3 Wave functions 

The eigenfunctions at points off the boundary, ^(rg ^ F), are determined 
by the null vectors u corresponding to the roots of the spectral determinant. 
Prom equation (4.4) we obtain immediately an integral representation of the 
(un-normalized) wave function. 



V'(ro) = ± 



r b 



u 



(4.22) 



for To ^ r. According to eq. (4.4) the integral vanishes identically either in 
the interior or in the exterior. This is indeed confirmed by our numerical 
calculations which are reported in the next chapter. 

In order to calculate the current density (2.32) one needs the gauge invariant 
gradient of the wave function. An integral formula is obtained from equation 
(4.4), after applying the differential operator V^o/b + iAo, 



V,.o/6V(ro) - iA(ro)V'(ro) = ± 



dr 

rT 



± + iAo)(9n/6G^ - iinG^) 



~ iyro/h^v + iAoG,,) 



u 



(4.23) 



The densities of other observables can be obtained by similar boundary inte- 
grals. 



4.3 The boundary operators 



In the following, we give explicit expressions for the boundary integrals. This 
allows to define the "finite part integral" appearing in the double layer equa- 
tion (4.9). 
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The integral operators (4.10) - (4.13), 



(QM) (ro) = J^dr q(r; ro)u{v) , (4.24) 

are defined by their integral kernels q(r; Tq). The form of the Green function 
(2.61) leads to the expressions 

q°(r;ro)=E(r;ro)G°(^) (4.25) 
qN(r; ro) = E(r; ro) | - i ^''^"^^^ G° (.) + .^G° (.) } (4.26) 



qS(r; ro) = E(r; ro) { - i^l^H^ G°(.) - 2^^^^^ z^GUz) ~ 



(4.27) 

q?,(r; r„) = E(r; ro) { (- ' ^""^ "f ^""^ - i^) G;(.) 

((r-ro)n)((r-ro)no) , pO/ x\ ^on 

with fi = n(r), yiq = n(ro), ^ := (r — ro)^/^^, and the abbreviation 
E(r;ro) :=exp 



-i ( ^ - X(r) + X(ro) 



(4.29) 



for the gauge dependent part. Note that the gauge freedom X has cancelled 
in the prefactors and appears in the phase only. It can be absorbed by the 
substitution ■u(r) — > exp(+iX(r))M(r), proving the manifest gauge invariance 
of the boundary integral equations (4.7), (4.9). Note also that expressions 
(4.26) and (4.27) are related by a permutation of r and tq with subsequent 
complex conjugation (since G° is real), hence the operators are the adjoints 
of each other. The self-adjoint nature of (4.25) and (4.28) follows hkewise. 

The derivatives appearing in (4.26) - (4.28) may be stated in terms of the 
gauge independent part of the Green function, G°, itself, at different energies u. 
They are given in Section 2.5.3 together with their asymptotic properties. 
displays a logarithmic singularity as r — > ro, while the differential expressions 
are bounded. In that limit, most of the quotients vanish for a smooth boundary, 
others tend to the curvature (3.3) at the boundary point ro- As a consequence, 
all the terms in (4.25) - (4.28) are integrable — but for the one containing 
the (nno)/(r — tq) ^-singularity. The latter gives rise to the need for a finite 
part integral. 
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The hypersingular integral operator 



For finite A the double-layer equation contains a hypersingular integral defined 
as 



f dT 

QdlM = -^i^no/b + iA„o)(^n/bG^ - iAnG„)u 

r dr 

:=lim / -rridno/b + iA„J(9„/f,G^ - iv4„G^)M 



(4.30) 



We want to replace the integrand by its limiting form. To this end the bound- 
ary is split into the part 7c£, which lies within a (ce) -vicinity around tq (with 
arbitrary constant c), and the remaining part Fee, 



lim 



dF 
dF 



+ / u(^^o/b + iAno){dn/bG^ - iA„G^)(M - Mo) 

'7ce 

f dF 

+^0 / -7Ti9no/b + iAno){dn/bGl-iAnGl) 
J "fee 



(4.31) 



with uq := u{ro). For sufficiently small e the boundary piece may be 
replaced by its tangent^ and the Green function by its asymptotic expression, 
cf Sect. 2.5.3. This way the third integral in (4.31) may be evaluated to its 
contributing order, 



dF 



(dno/b + iAno){dn/bGl - iA„G^ 



'Ice 

1 /rono 

+ 0(eHoge) 



/ no X ro 

— ^ 



1 r"^ - 



1 1 



27r J-^ (s2 + £2)2 



+ 0(£'log£) 



+ £2 



1 1 



TT ce C 



2\2 



ds 



0(5^ loge) 



VO{e logs) 

71 ce 



(4.32) 



Here, the explicit form of the integrand was obtained from (4.28) by the re- 
placement fq — > Tq . The last approximation in (4.32) holds because c may 
be chosen arbitrarily large. In a similar fashion it can be shown that the sec- 
ond integral in (4.31) is of order 0(£:). In the first integral we may replace 
(again for large c) the integrand by its limit, because e is small compared to 



We emphasize that we assume the bomidary to be smooth throughout. 
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min(|r — ro\) = c e. Therefore, the hmit in (4.30) may be expressed as 



where we replaced ce by e. This equation defines the finite part integral. It 
completes the derivation of the boundary integral equations. 



4-4 Solving the integral equations 

As discussed above, the integral equations (4.7) and (4.9) determine the spec- 
tra and wave functions of arbitrary interior and exterior magnetic billiards. In 
the stated form the equations are not yet suitable for numerical evaluation, 
though, since the integral kernels display (integrable) singularities. 

Fortunately, it is possible to treat the singular behavior analytically which ren- 
ders a highly accurate and efficient numerical scheme. In brief, the boundary 
integral equations are regularized using the known asymptotic behavior of the 
Green function and its derivatives, cf Sect. 2.5.3. Representing the periodic 
boundary functions in a Fourier basis then leads to an exponential localization 
of the integral kernels. This permits a well controlled truncation of the cor- 
responding matrix. The roots of the (Fredholm) determinant are accurately 
obtained by singular value decomposition. We refer the reader to our recent 
publication [35] for the technical details and a convergence analysis. p"^ 



5 Results of the boundary integral method 

The numerical implementation of the boundary integral method provides thou- 
sandths of eigenfunctions at high accuracy with moderate computational ef- 
fort. This includes the bulk states as long as the small energy difference to the 
Landau level can be represented numerically. In the following we demonstrate 
the performance of the boundary integral method by exhibiting numerical re- 
sults on magnetic billiards which have been inaccessible by other methods [35]. 



Note that the equations in [35] are stated in complex conjugated form since the 
focus is there on the wave functions rather than the Green function. 





1 



(4.33) 
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Fig. 5.1. Definition of domain boundaries considered in Chapter 5. The magnetic dy- 
namics in the asymmetric stadium (left) exhibits an anti-unitary symmetry, but no 
unitary one. In contrast, the skittle shaped domain (right) is free of any symmetry. 
It generates hyperbolic classical motion for p > 2. 

5.1 Spectral statistics 

Wc start by applying some of the standard tools of spectral statistics to large 
data sets of interior spectral points. The spectra are expected to reproduce the 
features of random matrix theory (RMT) if the underlying classical motion 
displays hard chaos [107]. In this section we define the spectra in the semi- 
classical direction 6 — > 0, keeping the cyclotron radius p constant. This way 
we can ensure that the classical dynamics are completely chaotic throughout 
the spectral intervals considered. 

We consider the two domains described in Fig. 5.1. One is an asymmetric 
version of the Bunimovich stadium billiard (ri = 0.75, r2 = 0.25, = 2.10957, 
=5f = 5.39724). In the magnetic field its dynamics is free of unitary symmetries 
but exhibits an anti-unitary one (time reversal and reflection at ?/ = 0). The 
skittle shape, in contrast, (made up of the arcs of four symmetrically touching 
circles, ri = 1.0, = 0.5, ^ = 4.33969, JSf = 9.42478) does not display 
any symmetry. It generates hyperbolic classical motion even for fairly strong 
magnetic fields [14]. (The asymmetric stadium is not strictly hyperbolic, but 
any possibly regular part in phase space is much smaller than the uncertainty 
product {b'^nY throughout the considered spectral interval.) 

We calculated 12300 and 7300 consecutive interior Dirichlet eigenvalues at p = 

1.2 for the asymmetric stadium and the skittle shaped domain, respectively. 
Using the boundary integral method it is possible to converge states even with 
much greater quantum numbers [35]. The time consuming task is really to find 
all energies, including the near-degenerate ones, in a given interval. 
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Fig. 5.2. Fluctuating part of the spectral staircase in the asymmetric stadium at 
p = 1.2. The displayed range contains the first 12000 points of the interior spec- 
trum, with the heavy line a running average over 500 neighboring points. A missing 
spectral point would show up as a distinct step by one. 




Fig. 5.3. Nearest neighbor spacing distributions of (a) the asymmetric stadium and 
(b) the skittle shaped domain (right), at p = 1.2. The histograms should be com- 
pared to GOE and GUE predictions of random matrix theory, respectively (dashed 
lines.) The corresponding cumulative probabilities are given by the monotonic lines. 
Here, the differences between data and RMT are of the order of the error of Wigner's 
surmise. 
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The integrity of the obtained spectrum may be checked by calculating the 
fluctuating part N^scl'^) = N'^'')(z/) — 'N^^\iy) of the spectral counting function. 
This quantity must average to zero indicating whether spectral points were 
missed. It is defined in terms of the mean staircase (given in equation (A. 53) 
for fixed p). Figure 5.2 displays N^^). for the asymmetric stadium. The strongly 
fluctuating function indeed vanishes on average which indicates the complete- 
ness of the spectrum. This can be seen from the heavy line which gives a 
running average over 500 neighboring points. The oscillations of the running 
average can be related semiclassically to the existence of bouncing ball modes, 
which are discussed below. A very similar result like Fig. 5.2 is obtained for 
the skittle shaped domain (not shown). 

The large spectral intervals at hand permit the direct calculation of some 
of the popular statistical functions used to characterize spectra. Due to the 
underlying classical chaos and the symmetry properties mentioned above one 
expects the statistics of the Gaussian Orthogonal Ensemble (GOE) for the 
asymmetric stadium, and of the Gaussian Unitary Ensemble (GUE) for the 
skittle. Figure 5.3 shows the distributions of nearest neighbor spacings P{s) 
of the unfoldedp^ spectra. Indeed, one finds excellent agreement with random 
matrix theory. The differences between the numerical and the RMT cumulative 
functions /(s) = Jq P{s')ds' stay below 2% (ie, below the error of Wigner's 
surmise [107]). 

In order to characterize the spectrum more sensitively one often considers 
the form factor K{t), ie, the (spectrally averaged) Fourier transform of the 
two-point autocorrelation function of the spectral density [108,109]. Figure 
5.4 gives the spectral form factors of the asymmetric stadium and the skittle 
spectra. The thin and heavy lines correspond to different degrees of averaging, 
using a spectral window of width 3 and 30, respectively. The RMT predictions 
are shown as dashed lines, and one observes again very good agreement. Since 
most of the other popular spectral functions, such as Dyson's A3 statistic, are 
transformations of the form factor we do not present them here. 

We emphasize that the good agreement with RMT is not only a consequence 
of the large spectral intervals the statistics are based on. It is as important 
to have the spectra defined at fixed classical dynamics. Had we calculated the 
spectra at fixed field, they would have been based on a classical phase space 
that transforms from a near-integrable, time-invariance-broken structure to a 
hyperbolic time-invariant one as p increases with energy. This transformation 
of spectral statistics from GOE to GUE as the field is increased was studied 
in [67-69]. 



The spectra are transformed to unit density; see also the discussion in Sect. 8.1. 
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Fig. 5.4. Spectral form factor of (a) the asymmetric stadium and (b) the skittle 
shaped domain based on 12300 and 7300 spectral points, respectively. The heavy 
lines belong to the same data after stronger spectral averaging. One observes good 
agreement with the random matrix prediction of the Gaussian Orthogonal and the 
Gaussian Unitary Ensemble, respectively, indicated by the dashed lines. 

5.2 Wave functions in the interior and in the exterior 



The skittle 

To get an overview of the various types of eigenstates one may encounter 
in magnetic billiards we proceed to present a selection of stationary wave 
functions. We focus on the semiclassical regime of large scaled energies at 
cyclotron radii small enough to observe strong effects of the magnetic field. 
We start with the skittle shaped domain choosing again p = 1.2, such that 
the corresponding classical skipping motion is chaotic in the interior as well 
as in the exterior. 

Figure 5.5(a) shows the density plot of a typical interior wave function around 
the one-thousandth eigenstate. As expected for a classically chaotic system it 
spreads throughout the whole domain and has the features of a random wave. 
Occasionally, one also encounters so called bouncing-ball modes. These states 
are localized on a manifold of marginally stable periodic orbits (which have 
zero measure in phase space). In the skittle a prominent manifold consists of 
orbits with period 2 bouncing in the larger circular part of the billiard. The 
wave function of a corresponding bouncing-ball mode is given in Fig. 5.5(b). 

We turn to the eigenstates of the exterior billiard. The wave function of a 
typical example may be found in Fig. 5.5(c). It belongs to an energy close to 
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Fig. 5.5. Interior and exterior wave functions of the skittle shape around the 
one-thousandth interior eigenstate, at fixed cyclotron radius p = 1.2. The plot- 
ted shade is proportional to \^p\ and the thin line indicates the boundary F. (a) A 
typical interior wave function, v ~ 32.9880. (b) A bouncing ball mode, v ~ 33.1203. 
(c) A typical exterior edge state, u ~ 32.8474. (d) A typical exterior bulk state, 
u 2± 32.50025. [figure quality reduced] 
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Fig. 5.6. (continued) (e),(f) The same exterior edge and bulk states as in (c) and 
(d), respectively, on a larger scale. (g),(h) The current densities of the edge and 
bulk states (e) and (f), respectively. Here, the length of the arrows is proportional 
to the magnitude of the current density (measured at the at the tails of the arrows) . 
[figure quality reduced] 
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the one of Fig. 5.5(a) and is displayed on the same scale. Again, one observes 
the typical features of a chaotic wave function. When viewed on a larger scale, 
cf Fig. 5.5(e), we find that this state seems bound to the billiard and vanishes 
rapidly after a distance smaller than a cyclotron diameter. In addition, circular 
structures are faintly visible in the probability distribution with a radius given 
by the classical cyclotron radius p = 1.2. This state is clearly an edge state 
corresponding classically to a skipping motion around the billiard. This is also 
evident from Fig. 5.5(g) which displays the distribution of the current density 
of the state. 

Figure 5.5(d) shows a quite different exterior state. Its energy is still in the 
same range as that of Fig. 5.5(c) but now close to a Landau level. One observes 
that, unlike the edge state, the wave function shows no appreciable amplitude 
close to the boundary. Moreover, it displays a rather regular structure consist- 
ing of rings of maximal probability density which encircle the billiard. This 
is seen clearly on a larger scale, cf Fig. 5.5(f). The band running around the 
billiard has a width of the cyclotron diameter, 2p, and in general consists of 
N -\- 1 rings if the energy is close to v = N + ^ (here N = 32). This band 
moves outwards and gets more circular as one goes to bulk energies which are 
increasingly close to the Landau level (this way the bulk states in the sequence 
sweep over the whole plane). Clearly, we are dealing with a bulk state. Its wave 
function corresponds to a superposition of unperturbed cyclotron orbits which 
are placed around the billiard. This view is supported again by the distribu- 
tion of the current density, cf Fig. 5.5(h). Near to the boundary it displays 
an opposite orientation compared to at a distance of 2p. This renders the the 
net current around the billiard exponentially small — unlike the edge state 
Fig. 5.5(h) which displays a large finite current. 

In order to see this separation into edge and bulk states more clearly we turn 
to even more semiclassical energies and a symmetric shape of the boundary. 

The ellipse 

Next, we choose an elliptic boundary (of eccentricity 0.8 and area tt) and 
even more semiclassical energies. The cyclotron radius is taken to be p = 0.6 
which is small enough for complete cyclotron orbits to fit into the interior 
domain of the billiard. The classical dynamics of the skipping motion is mixed 
chaotic in this case [49], see Fig. 3.1 for a phase space portrait. Going to 
the extreme semiclassical limit - the ten-thousandth interior eigenstate - we 
expect the wave functions to mimic the structures of the underlying classical 
phase space. Indeed, Figure 5.7(a) displays a wave function which is localized 
along a stable interior periodic orbit. This orbit has period 36 traveling six 
times around the billiard, with six reflections each time, before it repeats 
itself. In contrast, the state in Fig. 5.7(c) is localized on a (mixed) chaotic 
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Fig. 5.7. Probability densities (left) and current density distributions (right) in an 
elliptic domain at p = 0.6, around the ten-thousandth interior eigenstate. (a,b) An 
edge state localized along a stable periodic orbit, v ~ 60.0602. (c,d) A state which 
covers the whole domain ~ 60.1664. Note that the gross direction of the current 
changed compared to (b). (e,f) An interior bulk state, u ~ 60.50031. The wave 
function is exponentially small close to the boundary and represents a superposition 
of cyclotron motion, [figure quality reduced] 
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Fig. 5.8. Exterior probability densities (left) and current density distributions (right) 
at p = 0.6 and at similar energies as in Fig. 5.7. For each state only the first quarter 
of the picture is shown. (a,b) A typical edge state, u 2± 60.2087. (c,d) An edge 
state localized along a stable periodic orbit, u ~ 60.2220. (e,f) A typical bulk state, 
1/ ~ 60.50016. [Possible Moire patterns in (e) are an effect of low resolution printing.] 
[figure quality reduced] 
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part of phase space. Wc note that the gross circulation of the current density, 
Fig. 5.7(d), is opposite to that of the first state, Fig. 5.7(b). Since p is small 
enough for closed cyclotron orbits to fit into the ellipse we find bulk states also 
in the interior, see Fig. 5.7(e) for an example. Again, it almost vanishes at the 
billiard boundary and may be viewed semiclassically as due to a superposition 
of closed cyclotron orbits. This view is supported by the distribution of the 
current density displayed in Fig. 5.7(f). 

Similar states are also found in the exterior, as displayed in Figure 5.8. To 
show more details we give only the righthand-top quarter of the figure (the 
others follow by symmetry). The first edge state. Fig. 5.8(a), corresponds to 
a classical motion with creeps along the boundary. It is the analogue of a 
whispering-gallery mode. Figure 5.8(c) displays an exterior edge state which 
extends much further into the plane. Like Fig. 5.7(a) it is clearly localized 
on a stable skipping periodic orbit. The bulk state Fig. 5.8(e) is close to 
1/ = 60 + 1. It consists of 61 concentric rings of increased amplitude and shows 
no appreciable net current around the billiard. 

We emphasize that all the interior and exterior wave functions shown above 
are calculated throughout the entire displayed area. They turn out to be nu- 
merically zero in the complementary domains (ie, either the exterior or the 
interior) as expected from theory. 



5.3 General boundary conditions 
The Neumann ground state 

So far, we only considered Dirichlet boundary conditions. They are the natural 
choice from a physical point of view if one considers the billiard boundary as 
due to an infinite wall potential. On the other hand, the Neumann boundary 
conditions, A^^ = 0, are frequently employed in spectral theory [22]. They 
have the advantage that the ground state energy lies below the first Landau 
level, which facilitates its mathematical analysis. Here, we are able to observe 
the manifestation of a recent theorem of spectral theory [110]. It states that 
the Neumann ground state of a magnetic billiard is exponentially localized 
around the boundary point of maximum curvature. In order to deal with a 
unique boundary point of maximum curvature we choose the union of a half- 
circle and a half-ellipse (with half-axes a — 2,b — 0.5) as billiard boundary. 
Choosing a magnetic length of 6 = 0.05 (which corresponds to a very strong 
field), we find the ground state energy v = 0.2763. The Figures 5.9 and 5.10 
display the ground state wave function in the billiard and on the boundary, 
respectively. Indeed, one observes an exponential localization over six orders 
of magnitude. 
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Fig. 5.9. Contour plot of the ground state wave function (absolute value on a linear 
scale) for Neumann boundary conditions and strong field, h = 0.05, u 2± 0.2763. 
The wave function is localized at the boundary point of maximum curvature, as 
predicted by a recent theorem [110]. Here the billiard domain is given by the union 
of a half-circle and a half-ellipse (dashed line). 




Fig. 5.10. Boundary values of the ground state wave function of Fig. 5.9. As predicted 
in [110] it is localized exponentially at the point s = of maximum curvature. 
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Fig. 5.11. The parametric dependence of the exterior spectrum on the boundary con- 
dition (for the asymmetric stadium, ie =Sf = 5.39724, at fixed b = 0.25). The param- 
eter A interpolates between Neumann (arctan A = — f ) and Dirichlet (arctan A = 0) 
boundary conditions. The right graph shows details around the fourth Landau level, 
[figure quality reduced] 

Parametric dependence on the mixing parameter 

As a last point, we show the parametric dependence of a spectrum on the 
type of boundary conditions. Figure 5.11 presents the exterior spectrum of the 
asymmetric stadium as a function of the scaled mixing parameter A G (— cxd, 0], 
cf (3.7). It is chosen negative to ensure that the transformation from Neumann 
(A = — oo) to Dirichlet (A = 0) boundary conditions is continuous. For positive 
A this would not be the case, which is a restriction similar to the one for the 
field free case [76]. ( We use the the arctan function in Fig. 5.11 to transform 
the infinite range of A into a bounded interval.) 

One observes that all the energy levels increase monotonically as A is increased. 
The energies clustering around the Landau levels v — N -\-\,N ^f^Q belong to 
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bulk states. One observes that they are hfted from the Landau levels to higher 
energies at Dirichlct boundary conditions, whereas in the Neumann case they 
are always shifted to smaller energies. A semiclassical theory which describes 
the exponential approach of the bulk states to the Landau levels and their 
behavior as a function of A is given in Appendix A. 6. 3 and Section 7.3.1. We 
shall come back to Fig. 5.11 not only there, but also in Chapter 7, when we 
define the edge state density. 



6 Semiclassical Quantization 

In Chapter 4 the boundary integral equations were found to yield an efficient 
method for obtaining the exact quantum spectrum of magnetic billiards. It 
will be shown in the sequel that the same equations are as important for the 
semiclassical quantization: They serve as the starting point for the derivation 
of the semiclassical trace formulas. 

Periodic orbit formulas for magnetic billiards 

The celebrated trace formulas of Gutzwiller [111,27] and Berry & Tabor [112,113] 
allow the semiclassical quantization of systems in terms of their classical mo- 
tion. They were derived assuming a continuous Hamiltonian fiow. The corre- 
sponding formulas for field free billiards are known to exhibit additional phase 
factors which account for the billiard boundary conditions. In order to show 
that the same holds for magnetic billiards we shall explain how the correspond- 
ing trace formulas are obtained from the exact boundary integral formalism. 
To our knowledge no such derivation has been published for magnetic billiards 
to date. The natural approach is to follow the hues of Bahan and Bloch's 
treatment of field-free biUiards [114,115], in analogy to the surface-of-section 
method [116] and the scattering approach [117] for non-magnetic systems. 
Those attempts failed so far for magnetic billiards due to the appearance of 
an abundance of unphysical "ghost" orbits which could not be handled. To re- 
solve this problem we take advantage of the analysis performed in Chapter 4. 
There it was found that the boundary integral equations include spurious solu- 
tions which belong to a particular complementary problem. We will show that 
the semiclassical spectral determinant can be factorized, accordingly, into an 
interior and exterior part. Each of them leads to a trace formula incorporating 
only the physical periodic orbits in the appropriate domain. 

Like in the field- free case [114-120] the semiclassical quantization will be based 
on the double layer boundary integral equation. Apart from the spurious so- 
lutions, the main complication arising at finite magnetic field is the inherently 
asymmetric form of the respective integral kernel. Unlike the case of field-free 
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billiards [60], the latter is not simply related to the semiclassically unitary 
map operator derived from the generating function. 

In Section 6.1 we deduce the semiclassical approximants to the boundary 
integral operators of Chapter 4. After that, in Sect. 6.2, special map operators 
are introduced in order to transform the spectral function of the double-layer 
boundary integral equation. As a result, the number counting function is given 
in terms of the traces of powers of the map operators. The traces are evaluated 
semiclassically in Section 6.3 assuming hyperbolic skipping motion. We show 
why only classically allowed periodic orbits contribute and how their stability 
properties enter. The section concludes with the trace formula for the density 
of states and the magnetization density. In Section 6.4 the traces are evaluated 
assuming integrable dynamics. As a result we obtain the explicit periodic orbit 
formula for the spectral density of states in the magnetic disk billiard. Section 
6.5 gives the corresponding WKB solution. 



6.1 The semiclassical boundary integral operators 



In Section 4.3 the boundary integral operators were defined in terms of the 
free Green function and derivatives thereof. To obtain the semiclassical ap- 
proximations of the operators one simply replaces the Green function by its 
approximant. The latter is the leading order asymptotic expression in the 
semiclassically small parameter which was derived in Sect. 2.5.1. To re- 
main at a consistent level of approximation, the derivatives appearing in the 
single-layer Neumann and the double-layer operators (4.11) - (4.13) are to be 
evaluated to the same leading order. This means in practice, that only the 
phase of the Green function (2.58) must be differentiated. Accordingly, in the 
remainder of this report all equalities involving semiclassical quantities are 
understood to be semiclassical in the sense that they hold to leading order in 



In order to obtain expressions which have a semiclassically intuitive and useful 
form it will be important to use the representation (2.58) which contains the 
actions of the short and long arcs separately. We found the geometric parts of 
the corresponding scaled actions (2.56) to be given by 



/ . /^|r-ro|\ , |r-ro| 



as(r; fq) = - arcsin — + 



and 

aL(r;ro) = 1 - as(ro;r) . (6.2) 
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Fig. 6.1. The angles a, (3^, and /3 are defined with respect the vector r— fq connecting 
the initial and the final point. They measure its length and the relative direction of 
the boundary normals, respectively. These quantities do not depend on the type of 
the arc (left: long, right: short), unlike the relative direction of the initial and the 
final velocities and v. The latter may be expressed in terms of a, and P, cf 
(6.7). (The dotted line indicates the billiard boundary.) 

As a first step, we note their gradients with respect to the initial and the final 
points. 



Vroa^S)(r; Tq) 



1 

np 



ro 



2p 



Vr a/ s ■ 



[r; ro) = — 

vrp 



ro 



A 



ro 



2p 




(6.3) 



(6.4) 



Here, the upper and lower signs of the first summands stand for the short arc 
and long arc contribution, respectively. It will be useful to state the distance 
between the initial and the final point in terms of the positive angle 



a[r;ro) 



arcsm 



rol 



2p 



(6.5) 



In addition, the direction of the normal vectors at the initial and the final 
points are measured by their (signed) angles with respect to the distance 
vector connecting the two points. 



P{r; ro) := <(n; r - ro) /3°(r; ro) := <(no; r - ro) (6.6) 

Now consider the classical arcs connecting the initial and the final points. They 
define the direction of the classical velocities at the points of refiection and 
incidence (see Figure 6.1 for a sketch of the situation). The normal components 
are given by 

vs fi = cos{P + a) Vg no = cos(/3° — a) 

Vl n = — cos(/9 — a) Vl fio = — cos(/?° + a) (6.7) 
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for short and long arcs, respectively. Here, the velocity at the initial point 
of the arc is denoted with a zero superscript, and the hats indicate that the 
velocity vectors are normalized. 



The semiclassical Dirichlet operators 

We proceed to calculate the semiclassical approximation to the kernel of the 
double layer Dirichlet operator (4.12) by inserting (2.58). For the short arc 
term one has to evaluate the gauge invariant derivative 



A ro • r ( ~ ^o) iio (r-ro)xno\ 
dno/b[27fwas - Xo\ + i^no = ^W^ , r- cos(a) 

= -2iv^ cos(/5° -a) = -2iv^ (v^ fio) , (6.8) 

where we used eqs (6.3) and (6.7). Apart from the sign, it is given by the normal 
component of the classical velocity after reflection since 2^/l' is the magnitude 
of the scaled velocity. Note that (vg fio) is a non-symmetric function of the 
initial and the final point and depends on the energy through a. For the long 
arc term one obtains the analogous expression 

9„o/427riz/aL - Xo] + i^no = +2iv^ cos(/3° + a) = -1\^fv (vl ^o) • (6.9) 

It follows that the semiclassical approximation of the double-layer Dirichlet 
kernel (4.27) can be stated in a particularly simple form, 

„D{sc)/ N — o r;(sc) I - A r;{sc) ^_ pi{x-xo) 

<4dl ir, roj .- Cno/folor^ lor^ ~ 2 (1 + e^™) (27ri) 5 

i (sin a COS a) 2 (sin a cos a) 2 J 

It will be an important ingredient in the derivation of the trace formulas.P^For 
completeness we note that the semiclassical single-layer Dirichlet kernel is 
simply given by the semiclassical Green function itself, q2*^'^'^^(r; rg) := G^f'^^ 
as an immediate consequence of (4.10). 



The semiclassical operators derived in Sect. 6.1 allow the computation of spectra 
within a "semiquantum" approximation by means of the boundary integral method 
of Chapter 4. One merely replaces the exact kernels (4.25) - (4.28) by their ap- 
proximants (2.59), (6.10), (6.13), (6.15) and calculates the respective determinants 
numerically without further approximation. As an advantage of this scheme it ap- 
plies irrespective of the type of classical motion (including mixed chaotic dynamics). 
However, it may not be regarded as a proper semiclassical quantization since the 
degree of approximation is not consistent throughout the calculation. Also the Fred- 
holm determinant must be evaluated to leading order in i/, see Sections 6.3 and 6.4. 
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The semiclassical Neumann operators 



The kernels of the single- and double layer Neumann operators, eqs (4.11) and 
(4.13), involve gauge invariant gradients with respect to the final point of the 
Green function. One finds 

9„/6[27rii/as + x] - i^n = +2iv^ cos(/3 + a) = +2iv^(vs n) , (6.11) 
dn/b[2niu&L + x] - iin = -2iv^ cos(/3 - a) = +2iv^ (vl n) , (6.12) 

similar to eqs (6.8) and (6.9). This way the semiclassical single-layer Neumann 
kernel (4.26) assumes the form 



qf""\r; Tq) := dn/bG^J'^ - iinGif^^ = ^ , . , -^-^ 
Hsl V > n/b u n u 2(1 +e2'^^'^) (27ri)5 



gi(x-Xo) 



X ( , e^™ + e-i , e'---^- ] (6.13) 

[ (sin a cos a) 2 (sin a cos a) 2 J 

It is worth noting how the mutual adjointness of the operators (4.26) and (4.27) 
shows up in the semiclassical case. By permuting r and tq the prefactors of 
the short and long arc terms change their roles, 

(vs n) = vs(r; Fq) n(r) = (fq; r) n(ro) = (v° Aq)^ , (6.14) 



and likewise (vl n) = (vg fio)^. As for the phases, it is the factor (1 + c^^^'^)^^ 
whose conjugation provides the term needed in conjunction with the 

relation (6.2) to prove the mutual adjointness. 

The kernel of the semiclassical double-layer Neumann operator follows from 
applying the gauge invariant derivative (6.11) to the single-layer Dirichlet 
expression (6.10), cf (4.13). One obtains 

qS^^^V; ro) := {d^^/b + iA^,){dn/bG^:'^ - iinG^^^)) 



2(1 + e2--) (2vri)l 



g-(x-xo) 



^ 1 -(vgfio)(vsA) ^ -Kfto)(vLA) I _ 

t (sin q; COS a) 2 (sin a cos a) 2 J 

Like the exact kernel (4.28) this semiclassical version is self-adjoint. This fol- 
lows again from the observation that the two summands simply change roles 
when the adjoint operator is formed. 
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6.2 From boundary to map operators 



Let us now consider the semiclassical double-layer equation for Dirichlet bound- 
ary conditions in more detail. As known from Chapter 4 the corresponding 
Fredholm determinant (4.15) is a spectral function. Its roots yield the Dirich- 
let spectrum of the domain considered, conjoint with the Neumann spectrum 
of the complementary domain. 

The semiclassical spectral function is obtained by substituting (6.10) into 
(4.15) with A = 0: 



^M(z/) :=det 



(6.16) 



Here, the upper or lower sign indicates whether the originally considered do- 
main is of the interior or the exterior type (like in Chapter 4) . 



The map operators 

By defining the operator 

P :=2(l + e2'^'")Qdi^''^ (6.17) 
wc factorize the determinant 6.16 into two parts, 

^(-)(zy) = det 

This reflects the partitioning of the underlying classical phase space into cy- 
clotron orbits, which are detached from the boundary, and skipping trajec- 
tories, see Sect. 3.1. The first determinant in (6.18) does not depend on the 
boundary. Its operator is diagonal and singular at the energies I'jsf = N + ^, 
G No, of the Landau levels. Apparently, it represents the semiclassical con- 
tribution of the bulk states to the spectrum and its divergence at the Landau 
energies is due to the infinite number of degenerate bulk states found in the 
exterior. (The exponentially small lifting of the degeneracy observed in the 
exact spectrum is not seen here since the semiclassical Green function (2.59) 
does not describe tunneling effects.) 

As will become clear in the following, the second factor in (6.18) yields the con- 
tribution of the skipping trajectories to the spectral function. It is described 
by the map operator P defined in (6.17). We will see that it can be related 
to the classical billiard map (3.4) describing the motion of skipping trajecto- 
ries. The map operator consists of a short arc and a long arc term and it is 



2(1 + e2^''^) 



det 



1 + e^^''' ± P 



(6.18) 
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advantageous to split it accordingly, 



P = Ps + Pl , (6.19) 
with the corresponding integral kernels given by 

Ps(s, So) r e^™ e'*-'*° (6.20) 

(27ri)2 (sin(Q;) cos(a))2 



and 



Pl(s, so) r e'^^ e"^"""^ e'*"'^" , (6.21) 

(27ri)2 (sin(Q) cos(q))2 

see (6.10), with r = r(s),ro = r(so)- 

We note that the operator P differs from the standard map operator defined 
in terms of the generating function ^ of the classical map [60] , 



1 / 



(27ri)2 V^^^^o, 



(6.22) 



which is the analogue of Bogomolny's transfer operator T [116]. The main 
difference is not the fact that P consists of two distinct parts. Rather, it is 
the inherent asymmetry in the coordinates s, sq which introduces difficulties 
not encountered in the field-free treatment: The pref actors cannot in general 
be stated as mixed derivatives of the relevant phase. Moreover, the parts Ps/l 
are not semiclassically unitary (but satisfy equation (6.25)). 

From now on we focus on the non-singular factor in (6.16). It gives the spectral 
function 

eig(i.)=det[l + e2-'^±P] (6.23) 

of the states which correspond to skipping motion. In the standard procedure 
to obtain a trace formula one would now compute the imaginary part of the 
logarithm of (6.23). Making use of the identity log dot = trlog [121] one would 
like to evaluate the trace of powers of the operator (6.19) in stationary phase 
approximation. However, unlike the case of field-free billiards the correspond- 
ing saddle point condition selects classical periodic orbits in the interior and 
in the exterior. This is not surprising given the fact that the double-layer 
equation includes solutions of the complementary domain. To make matters 
worse, an abundance of saddle-point configurations arises which do not have 
a physical meaning at all. In order to avoid these severe difficulties it is vital 
to write the spectral function (6.23) as a product such that each factor yields 
the spectrum in either the interior or the exterior domain. 
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Factorizing the spectral function 



To facilitate the factorization of the determinant (6.23) we spht the short and 
long arc operators once more, Ps = P™* — Pg^* and Pl = p^* — pj!^*. Ultimately, 
the parts labeled by "int"and "ext" should exclusively account for the motion 
in the interior and in the exterior, respectively. To that end, the splitting is 
defined by the signs of the prefactors of the integral kernels which are functions 
of the initial and the final points. 



e(-noVs)Ps(-5, So) 
-e(novg)ps(s. So) 



pL'^*(^,^o) 

pr(^,^o) 



0(-noVL)PL(s, So) 
-e(novJ^)pL(s, So) 



(6.24) 



(The minus sign in front of the exterior kernels is introduced for convenience.) 
According to these definitions the "interior" part of the operators vanishes 
whenever the initial and the final points have positions such that the corre- 
sponding classical arc points into the exterior domain, and vice versa. This 
crucial property is embodied in the operator equations 



pint pext _j_ pint pcxt 



,27rii/ 



id 



(6.25) 



and 



pint pext _j_ pint pcxt 
S S L L 







(6.26) 



derived in Appendix A. 5. With their help it follows immediately that the 
determinant in eq (6.23) factorizes into an interior and an exterior part. 



(sc) / 
skipv 



det 
det 
det 



1 + e"^""'" ± (P^'^* - P^^* + P[^* - Pf 



1 ± (P'i' + 



1 ± P 



int 



det 



oxt 



(6.27) 



Here we merely replaced the term e^'^"^ by the operators (6.25) and included 
(6.26). In the last equality we introduced the interior and exterior map oper- 



ators P 



mt . _ 



PT + Pr and P*^ 



pcxt ^ pcxt^ respectively.[3 



The factorization of the spectral function (6.27) into an interior and exterior 
part is in accordance with the observation that the double layer equation 
provides the spectra of both the interior and the exterior problems, see the 
discussion in Sect. 4.2.2. In the following section we shall show that using 



13 "Y^e remind the reader that the signs it in (6.27) choose the domain of interest and 
stem from the boundary integral equations. The labels "int/ext" refer only to the 
type of the classical arcs which are included in the definition of the corresponding 
operators in (6.24). 
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the upper sign, ie, starting originally with the interior problem, one gets the 
trace formula for the interior Dirichlet spectrum from the first factor and the 
exterior Neumann spectrum from the second one. In the same way, if we use 
the lower sign we get the trace formula for the exterior Dirichlet spectrum 
from the second factor and the interior Neumann spectrum from the first. 

We finally note that the operators P"'* and P'^^' are semiclassically unitary, 
as shown in Appendix A. 5. Moreover, the interior and exterior map operators 
obey the relation 

pintpext _ _e2^i^id ^ (6.28) 

which follows from (6.25) and (6.26). It has an intuitive form: Propagating a 

boundary state first in the exterior and then in the interior one arrives again 
at the same state, augmented by the global phase 27rz/ of a complete cyclotron 
orbit (plus the Maslov correction tt). It is the semiclassical manifestation of 
the classical interior-exterior duality. 



6.3 Trace formula for hyperbolic billiards 



To obtain the fiuctuating part of the number counting function of the skipping 
spectrum one has to take the imaginary part of the logarithm of the spectral 
function (6.27). We start by computing the contribution of the first factor in 
(6.27), 



1 



Nosi^^^^^n^) = - - Imlogdet 



TT 



1± P 



int 



^ilmVi^^tr 

^ n=l 



(P^°*)"] (6.29) 



To obtain the periodic orbit formula we can now follow the lines of the deriva- 
tion of the trace formula for field free bilhards [114,1 16,115,99,118,119,122,123, 120] . 
The trace in (6.29) amounts to an n-dimensional integral of the form 



tr 



(27ri)«/2 



dsi . . . ds^ 



exp(i^(X(si+i)-X(s,)) 



^(v°no),e(-(v°no) 

(^sm(aj) co£ 
-(vOno),0 



exp ( 2iTwas{sj+i; sj] 



+ 



(si 



sm[aj) cos[aj 



- exp I 27rii/aL(sj+i; Sj) - i- 



)) 



(6.30) 



Here, the abbreviation (v°no)j := v''(r(sj+i); r(sj)) n(r(sj)) is used, together 
with (6.32) and the convention sq = Sn- Note that the gauge dependent factor 
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Fig. 6.2. The angles entering the jth saddle point condition. (The dashed line indi- 
cates the boundary.) 

(involving the Xi^j)) vanishes identically as a consequence of the cyclic per- 
mutability of the integration variables. This renders the trace (6.30) a gauge 
invariant quantity. It is now evaluated to leading semiclassical order using the 
stationary phase approximation (A. 30). 



6.3.1 The saddle point conditions 

For each of the 2" integrands in (6.30) the condition of a stationary phase 
leads to n saddle point equations 



_d_ 

d7, 



27rz/a^j.(sj; Sj-i) + 27rz/a,,^^, (s^+i; Sj) 



I 



0, je{l,...,n}. (6.31) 



Here, the indices rjj G {S, L} account for the 2" different sequences of short 
and long arc operators under the trace. We shall treat all these equations 
simultaneously by noting for any solution s_ of (6.31) not only the configuration 
of saddle points but also the corresponding sequence of types of arcs, s = 

((Si,r7i), .., {Sn,Vn))- 

In order to obtain a geometric interpretation of the saddle point conditions 
we fix the positions := r{sj) and extend the definition of the angles (6.5) 
and (6.6) to a sequence of n points. 

aj := arcsin (6-32) 
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= - sm(/3° - ttj) 








= +sin(/3° + q;j) 




= +cos{Pj + aj-i) 




= - sin(/3j + ttj-i) 




= — COs(/3j — ttj-l) 




= +sin(/3j - ttj-i) 



Table 6.1 

Components of the incident and reflected velocities. For the geometrical interpre- 
tation see Fig. 6.2. 

and 



■= <(ni;rj+i -r^) 
This definition imphes 



cos(/?0) 



(3j+i :- <(nj+i;rj+i - r^) . 



(6.33) 



kj+i-rjl 



sin(/?j) 
sin(/5°) 



and 



cos(/3j+i - = n^+i sin(/3j+i - = n^- x n^+i 



(6.34) 



(6.35) 



Again, aj determines the angles of the incident and the reflected velocity 
vectors with respect to the direction given by r^+i — r^. It follows that the 
normal and tangential components of the velocity are given by the expressions 
in Table 6.1. They allow to state the derivative of the action with respect to 
the arc length s along the boundary in a particularly convenient form: 




1 - 



Tj — rj_i)t,- t,- X r 



2p 

cos(Q;j_i) sin(/3j) + 



2p J 



2p2 



Similarly, one finds 



— as(sj+i; Sj) 



1 



ds 



-aL(sj;Sj-i) 



( + cos(q;j) sin(/3°) 
+ cos(Q;j_i) sin(/?j) + 



Tip 
1 



np 



— aL(sj+i; Sj) = — - cos(aj) sin(/?°) - 



ds 



np 



'i+i \ 
2p J 

2p , 

2p i 



(6.36) 

(6.37) 
(6.37a) 
(6.37b) 



73 



As a result, an explicit expression for the jth saddle point condition is obtained 
in terms of the vectors rj_i, r^, fj+i, and hj. Naturally, the condition depends 
on the type of the two operators involved. 



2p 



— sm{Pj) cos(q;j_i) + sin(/3°) cos(q;j) if (77^-1, rjj) — (S, S) 

— sin(/3j) cos(q;j_i) — sin(/3°) cos(aj) if (77^-1, rjj) — (S, L) 
+ sm{f3j) cos(Q;j_i) + sin(/3°) cos{aj) if (77^-1, 77^) = (L, S) 
+ sin(/3j) cos(q;j_i) - sin(/3°) cos(q;j) if (77j_i, 77^) = (L, L). 

(6.38) 



The left hand side of this equation can be written in terms of the angles appear- 
ing on the right side after adding and subtracting the expression (r^ hj)/{2p). 

(r,+i - r,) ^ (r,- - r,_i) n, ^ ^^^^ ^.^^^^^ ^ ^^^^^^^ sin(a,_i) (6.39) 



2p 2p 

Combining the last two equations, the saddle point condition assumes a form, 

sm{j3j + aj_i) = sin(/3° - aj) if (77j_i, 77^) = (S, S) 

sin(/?j + aj_i) = - sm{p'^ + aj) if {r]j_i, r]j) = (S, L) 

sm{Pj - aj_i) = - sin(/5j - aj) if {rjj^i, r]j) = (L, S) 

sin(/3j- - aj_i) = sin(/?° + a^) if (77j_i, 77^-) = (L, L) , (6.40) 

which should be compared to the expressions in Table 6.1 for the compo- 
nents of the classical velocities. One observes that the equations (6.40) simply 
amount to the condition 

t.-v,,,- = t, <.,. , (6.41) 

for j = 1, . . . n, and any rjj e {S, L}: The tangential component of the classical 
velocities which correspond to the saddle point configuration s are continuous 
in the point of reflection. Since the modulus of the velocity is a constant of the 
motion, the trajectory is either continuous in this point or the normal compo- 
nent changes its sign. In the first case the trajectory penetrates the boundary 
which we call an unphysical solution. In the second case, the trajectory corre- 
sponding to the saddle point configuration obeys the law of specular reflection 
in r^-. 

From the fact that (6.41) must be satisfied simultaneously at the n points Vj 
it follows that any saddle point configuration corresponds to a closed, peri- 
odic orbit. However, by no means is this orbit necessarily a physically allowed 
classical trajectory. Figure 6.3 sketches the two different types of saddle point 
configurations which appear in magnetic billiards. Here we choose n = 5, ie, 
the saddle points correspond to periodic orbits of period 5. Clearly, both of 
them are unphysical trajectories. The one on the top features a specular refiec- 
tion at r2. Then the boundary is penetrated at giving rise to a full cyclotron 
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Fig. 6.3. Typical saddle point configurations appearing in the semiclassical evalu- 
ation of the trace (6.30). Both correspond to unphysical trajectories. The upper 
configuration has no relation to an orbit of the classical problem, while the lower 
one corresponds to a physical trajectory which is dressed by an additional cyclotron 
loop. 

loop. After one more reflection (this time from the exterior) at = r2 the tra- 
jectory arrives at its initial point. It performs one more cyclotron orbit without 
even displaying a boundary point at r2. This saddle point is a legitimate solu- 
tion of (6.41) belonging to a dense and two-dimensional set of stationary points 
(since the boundary points ri and T2 may be shifted independently without 
changing the picture.) It has clearly no relation to a physical periodic orbit. 
The saddle point shown on the bottom part of Fig. 6.3, on the other hand, 
does exhibit the boundary points of a physical periodic orbit (with period 4). 
Nonetheless, the depicted trajectory is unphysical since it leaves the interior 
domain, performing a cyclotron loop between the third and forth boundary 
points. Obviously, there is an infinite number of these unphysical saddle points 
attached to any proper, physical periodic orbit. They merely dress the orig- 
inal orbit with additional cyclotron loops. It might be expected that these 
unphysical contributions can be re-summed, leaving behind only the contri- 
butions of physical periodic orbits of the interior and exterior problem. This is 
a difficult task, due to its combinatorial nature in conjunction with a number 
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Fig. 6.4. Sketch of an interior ghost orbit, see text. (The normals are pointing 
outwards.) 

of ambiguities. A saddle point configuration may, for example, incorporate an 
interior and and exterior periodic orbit at the same time, leaving the question 
undetermined whether to assign the contribution to the interior or to the ex- 
terior problem. These problems are resolved immediately by the splitting of 
the operator (6.24) into interior and exterior types. Here, it is the Heaviside 
functions introduced by the splitting which guarantee that only those saddle 
points contribute for which the corresponding classical trajectory is directed 
into the correct domain at each point of reflection. As a consequence, the un- 
physical solutions discussed above are erased from the sum. The remaining 
saddle points will be denoted by 7!"^ and 'yi^l, respectively. They correspond 
to the periodic orbits of period n found in the classical interior and exterior 
billiard problem. 

Strictly speaking, the set of saddle points 7*^"^ which are directed into the 
correct domain at each point of reflection still includes the so-called ghost 
orbits. These are periodic orbits which leave (and necessarily re-enter) the 
proper domain without exhibiting a component of the saddle point (ie, a point 
of reflection) when leaving it. The left side of Fig. 6.4 shows the situation. Like 
in the case of non-magnetic billiards [114] these saddle points finally do not 
contribute to the sum over the traces. This is because for any ghost orbit of 
period n one finds another of period n + 1, with an additional boundary point 
at the position of re-entrance (right side of Fig. 6.4). These two contributions 
differ by a factor (—1) due to the additional boundary point and therefore 
cancel. 

6.3.2 The pref actors 

The next step is to transform the prefactors in the trace integral (6.30). Even- 
tually they should combine with the determinant of the matrix of action 
derivatives introduced by the stationary phase approximation (A. 30). The 
resulting expression should then be given in terms of the physical properties 
of the attributed periodic orbit 7*^"-'. 

We start with the evaluation of the mixed derivatives of the actions in terms 
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of the angles (6.32), (6.33) characterizing the jth part of the trajectory. If the 
arc is of the short type one obtains the formula 

as(.,+„ .,) = (^cos(a,) sm(/3, ) - ^ j 

'^tan(Q;j) sin(/3j+i) sin(/3^) 



27rp' 



- cot(Q;j) cos(/?j+i) cos(/?") + sin(/?j-+i - /?") 

^ ^ (vs,-ftj)(vs,+iVi) _ .g_^2) 
27rp2 sin(Q;j) cos(Q;j) 

Here, the expressions given in (6.37) and Tab. 6.1 were employed, as well as 
^ cos(Q;j) = — sin(Q;j) — arcsin ( ■^-^-'-t^ — J (6.43) 



dsj+i dsj+i \ 2p 

■ tan(Q;j) sin(/3j+i) 



-1 sin(aj) {Yj+i-Yj)tj+i _ 1 



and 

d 



■sin(/?° 



2p cos(Q;j) |rj_|_i — Vj\ 2p 



kj+i-r^'l 



tj+i X n^- (r^+i - rj) x (r^+i - r^) t^+i 



(6.44) 



^ -1 cos(;9,+i) cos(/3°) 
2p sin(Q;j) 

Hence, the mixed derivative (6.42) is essentially determined by the normal 
components of the velocities at the initial and the final point of the corre- 
sponding arc. Note that this expression is manifestly positive if the arc is part 
of a physical trajectory meaning that it lies either in the interior or in the 
exterior at both points (see also Fig. 6.2). 

If the jth part of the action corresponds to a long arc one obtains in a similar 
fashion 

1 cos(/?° + aj) cos(/3j+i — ctj) 



■aL(,s,+i,s ' ^ 



dsjdsj+i 27rp2 sin(aj) cos(q;j^ 

1 (vLj%)(vLj+ini+i) 
27rp^ sin(Q;j) cos(Q!j) 



(6.45) 



The form of this formula is analogous to (6.42), except for the difference in 
sign. It follows that equation (6.45) is manifestly negative if the angles aj,f3j, 
and (3j+i describe a segment of a physical trajectory — again due to the change 
in the orientation of the velocity vector. 
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The mixed derivatives (6.42) and (6.45) allow the transformation of the prod- 
uct under the trace integral (6.30) into a symmetrized expression. For a given 
saddle point 7'^"^^ we denote the geometric part of the total action by 

n 

A(7W):=^a,^. (.,+!,. ,) , (6.46) 

cf (6.1), (6.2). In addition, the number of long arcs appearing in 7*-"^ will be 
called l^. The product under the trace integral (6.30) can now be stated in 
terms of A(7*^")) and For the interior operators, it assumes the form 

J=i ^sin(Q;j) cos(Q;j 



2 



nj=i Vi)) g2.i.A(7-0 (6.47a) 

n^=i (sin(a,)cos(Q;j))5 

1 

n"=i (- (vf^^. (v{^,+i ,-+i}n,-+i 



YVj=i (sin(aj) cos(Q;j))2 
^ -Q f sV^^^i^^^^^i^l ' e^i (^--^-) e2-^^(^i:t) (6.47c) 

.d^A(7S;J) 2 o,^i^.AM") 



(27r)in 



^^u^Hlint; g-if^,g2.i.A(7-) . (6.47d) 



Here we used several times the fact that the saddle point configurations de- 
scribed by 7int correspond to physical, interior periodic orbit with n reflections. 
First, we noted the positivity of the factor — (vJJ^.no)j to write it as a product 
of square roots (and shifted one index by one). Second, the reflection condition 



V 







was employed to get from (6.47a) to (6.47b). As a result, the prefactors are 
symmetric in Sj and Sj+i, which allows stating them in terms of the mixed 
derivatives of the classical action. Finally, given the sign of each factor in 
(6.47d) the sign of the product can be taken out. It is (— )^^, due to the 
long arcs in 

Upon evaluating the trace of the exterior operators one is led to the same 
expression (6.47d), with 7^^^ replaced by 7ext . This is because the additional 
sign in the definition (6.24) of the exterior operators cancels the change in the 
orientation of the normals relative to the velocity vectors. 
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p 


cyclotron radius (p > L)j 


(2.33) 


a 


relative distance or the initial and the nnal point (0 < a < ^) 


(6.5) 


P 


relative direction of normal at the initial point (0 < /J*^ < 27r) 


(a a\ 

(6.6) 


n 

p 


relative direction of normal at final point (0 < /3 < 27r) 


la a\ 

(6.6) 


as, aL 


geometric part of the action of a short (long) arc 


IR -\\ 

(6.1) 


7, 7^ ^ 


physical periodic orbit (with n reflections) 


page 76 


A(7) 


geometric part of the action of the periodic orbit 7 


(6.46) 




number of reflections (repetitions) in 7 


page 79 


/^7 


Maslov index (number of conjugate points in 7) 


page 80 


M(7) 


stability matrix of 7 


(6.50) 



Table 6.2 

Important geometric quantities 
6.3.3 Performing the trace 

Now, with the prefactors written as mixed derivatives of the action in (6.47d) 
we can follow the standard procedure to derive the semiclassical trace formula. 
We apply the stationary phase approximation to the n-dimensional trace inte- 
gral. At first, the contributing saddle points are assumed to be isolated. This 
amounts to the assumption that the corresponding classical billiard dynamics 
is hyperbolic [3] . The case of an integrable system is treated afterwards. 

Combining eqs (6.30) and (6.47), together with (A. 30), yields 

;(p'"T]= E 



tr 



n 1 lli=i P dsjdsi+i 



^7 P 

76 < 



det 

\ dskdsi J 



7e 



y e^-'^^^^') e-^^'^- . (6.49) 

^ - 'tr(M(7)-2|^ 



^7 



The factor n/r^ appears because the sum is taken over all n-periodic orbits 
of the interior billiard rather than over all contributing saddle points. Each 
n-periodic orbit (with repetition number ) corresponds to n/r^ distinct 
saddle points s, which are related by a cyclic shift of their components. For 
the last equality in (6.49) we used once more the fact that 7!^^ is a classical 
periodic orbit of a billiard problem. This implies a general relation between 
the derivatives of the generating function of the billiard map A(7("^)) and the 
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stability matrix M(7(")) [120], 



det 



, dskdsi 



k,l 



= (-r[trM(,<"))-2]n^''^(^'""' 



7=1 dsjdsj-i 



(6.50) 



Its modulus was taken to derive (6.49). The integer fi^ := + z/^ denotes the 
total number of conjugate points. Here, is given by the number of negative 
eigenvalues of the determinant in the dominator. It counts those conjugate 
points along the trajectory, which are due to the focusing and defocusing 
effect of the boundary. The remaining, trivial conjugate points, which show 
up at each long arc (after an angle of tt), are taken into account by i^. 

For later reference, let us mention that the dual partner orbit of 7'-"^ denoted 
as ^("■^ , has 



/ij — 2n — ^-y 



(6.51) 



conjugate points: As discussed in Section 3.2.3 the dual orbits consists of the 
arcs complementary to those of 7^"^ and has opposite orientation. Prom (6.2) 
we find A(7(")) = n — A(7(")) and it follows that z/y = n — z/^ since every 
element of the matrix of second derivatives in (6.49) is multiplied by (—1). By 
definition we have £y = n — £^ leading to (6.51). Note also that the stabilities of 
dual periodic orbits are equal, trM(7) = trM(7), which follows from equation 
(6.50). 



The trace formula for the spectral counting function 



Inserting the expression for the trace (6.49) into (6.29) we obtain the fluctu- 
ating number counting function attributed to the interior map operator. 



]\q-skip(int) 



1 °o fzc'l" 

ilmV IZL tr 



1 



TT 



TT 



n=l 

00 



l^pinty 



(6.52) 



n=l 



7e 



trM(7) 



-e 



27rii/A(7)-if/i^ 



E 



7e{7mt} r^ 



trM(7) 



TT 



T sin (27ri/A(7) - -/x^ 



(6.52a) 



It is naturally associated with the interior problem, since the sum includes 
all periodic orbits |7int} of the interior billiard problem (with the number 
of reflections). If we are originally interested in the Dirichlet spectrum of the 
interior billiard we have to choose the upper sign in (6.52). In this case each 
reflection is associated with an additional phase shift of tt. The lower sign is 
to be taken if the spectral problem was originally formulated for the exterior 
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spectrum. In this case (6.52) provides the spurious interior Neumann spectrum 
which is included by the double layer equation. It differs from the Dirichlet 
spectrum merely by the fact that there is no phase shift associated with the 
reflections at the biUiard boundary. 

The second factor of the spectral function (6.27) yields a number counting 
function which includes the trace over powers of the exterior operators. In 
complete analogy to the treatment above one obtains a periodic orbit sum 
like equation (6.49). As the only difference, the sum is over all the periodic 
orbits 7ext of the exterior classical billiard map, 

^skip(ext) _ i j2 ^ tr[(P^^*)"] (6.53) 

Y: r sin (^27rz.A(7) - . (6.53a) 



^ 7e{7ext} r. 



trM(7) 



Like above, an additional phase shift of vr is associated with each reflection 
if the original double layer equation was formulated for the same domain as 
the orbits are taken from (upper sign in (6.53a)). Again there is no shift if 
the periodic orbit sum represents the spurious solutions of the double layer 
equation which belong to the complementary domain (lower sign). The fact 
that the trace formulas for Dirichlet and Neumann boundary conditions differ 
only by a phase is also known from the theory of non-magnetic quantum 
biUiards [76]. 

We conclude that for either the interior or the exterior Dirichlet problem the 
fluctuating number function is given by 

^^n^) = - E — r sin (27TuA{j) - ttu, - ^/.,) , (6.54) 

^ r^|trM(7) - ^ ' ^ ^ ^ 

where the sum is over all periodic orbits in the respective domain. 

This final result in complete agreement with the standard trace formulas. 
One could have used them without the preceeding derivation. However, the 
exclusion of the non-physical trajectories would remain an act of faith. The 
derivation above provides a sound basis for the intuitively sound results. 



6.3.4 Geometric interpretation 

At this point a brief discussion of the geometric meaning of a trajectory's 
scaled action is in order. We start with the observation that the actions of 
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short and long arcs are given by identical expressions once the parameter 



^' pki+i-Til 

is introduced to describe the j-th arc. Unhke the angle aj (6.32), it is not 
just a function of and r^+i but it contains information on the type of the 
arc through its sign: aj is negative for short arcs and positive for long ones. 
The geometric parts of the actions of short and long arcs, (6.1) and (6.2), now 
assume the common form 

a(ri+i; Tj) = i + arcsin(aj) - aj^l-a] - ^^^^^^^ ' (^-56) 

which is a remarkable simplification.]^ It allows to show immediately that a 
periodic orbit 7 (of period n) exhibits a geometric action (6.46) 

A(7(")) = E a(r,^„ r,) = ^^^^ , (6.57) 
which is given by the length of the trajectory, 

n 

^1 P E 2 arcsin((Tj)) 
dz/ L J 

and the enclosed area, 

^/^ := ^/^P°^y T E ( 9 + arcsin(a,) + a^^\-a]\ . (6.59) 
j=i / 

Here, is the area of the polygon defined by the points of reflection {rj} 

and each of the summands in (6.59) is equal to the area enclosed by the j-th 
arc and the cord connecting its initial and final points, cf (2.44). (Overlapping 
parts of the enclosed area are counted according to their multiplicity.) 

Equation (6.58a) follows bearing in mind that p and o", are functions of u, cf 
(2.37). It illustrates the fact that the excursion time of a trajectory is given 
by the derivative of its action with respect to energy. Using the proper scaled 
energy E = 2v (cf the discussion of (2.37)) we obtain the scaled time of flight 



The derivation of the trace formulas would have been considerably more compli- 
cated, had we introduced this parameterization earlier. 



— cosl^ajj 
+ cos(ai) 



if "short" arc 
if "long" arc 



(6.55) 



(6.58) 
(6.58a) 
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T-y of the periodic orbit, 



d(2z/) I 



27ri/A(7(")) 



( ^ + arcsin((7j; 
i=i ^ ^ 



(6.60) 



Density of skipping states 

The formula for the fluctuating part of the density of skipping states (3.15) 
follows by taking the derivative of the number counting function (6.54) with 
respect to u, 



It. 



TT 



trM(7) - 2 



TT 



- cos I 27rz/A(7) — nn^ /x^ 



(6.61) 



It must be emphasized, however, that the applicability of this expression is 
rather restricted so far. It is valid only for the interior billiard and only if the 
entire phase space consists of skipping trajectories (ie, for weak fields only). In 
all other cases any attempt to include the cyclotron contributions "by hand" 
yields unsatisfactory results [71,124]. 



Magnetization density 

Another derivative of the action occurs in the definition of the scaled mag- 
netization density (3.26) which was discussed in Sect. 3.4. We find that it is 
determined by the area enclosed by the trajectory (6.59), 

= ±p=< ■ (6-62) 

The scmiclassical expression for the fluctuating part of the scaled magneti- 
zation density is obtained by applying the derivatives in (6.62) to the trace 
formula for Nqsc, cf (3.26). Assuming that all periodic orbits are isolated and 
of the skipping type we find 

m°«^(i/) = ±2 ^ s^^Kh'^T:) ^ [27ri/A(7) - Trn^ - ^/x^) . (6.63) 

7 r^|trM(7) -2|^ ^ ^ ^ 

Hence, compared to the density of skipping states (6.61) each periodic orbit 
contribution to the scaled magnetization density includes the enclosed area in 
units of 6^7r, ie, the magnetic moment of the classical orbit rather than the 
scaled time of flight. Again, the expression (6.63) is only applicable for the 
interior problem at weak fields. The corresponding, less intuitive scmiclassical 
expression for the conventional magnetization at weak fields may be found 
in [66]. 



d62 ^di/ 



27ri/A(7) 
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6-4 Trace formula for the integrable case 



In the previous section the classical billiard map was assumed to be hyperbolic. 
We now shift to the other extreme, the disk biUiard, which exhibits integrable 
motion. 



6.4-1 The disk biUiard 

The periodic orbit formula for the density of states in the interior of the 
magnetic disk was derived recently by Blaschke et al. [71]. These authors used 
the trace formula by Creagh and Littlejohn [125] to account for the continuous 
circular symmetry of the disk. 

In the following, we derive the trace formula starting from the boundary in- 
tegral equation. This demonstrates how the integrable case is treated in the 
framework of the boundary map operators and yields an explicit formula in a 
straightforward manner. Moreover, the exterior case is easily included in our 
treatment. 

Many results of the last section still apply. In particular, the factorization of 
the spectral function (6.27) does not depend on the type of motion, hence, we 
can start directly with the equations (6.52) and (6.53) for the interior and the 
exterior counting functions. However, the trace of powers of the map operators 
cannot be evaluated like in the hyperbolic case since the periodic orbits are 
not isolated but appear in continuous families [113]. 

The classical motion is governed by one parameter, the ratio 

Td := - (6.64) 
P 

between the radius of the disk R and the cyclotron radius. For weak fields, 
Fd < 1, any two points on the boundary can be connected in the interior 
only by short arcs and in the exterior (only) by long ones. The field is strong, 
Fd > 1, if complete cyclotron orbits fit into the interior. The skipping motion 
then displays both types of arcs in the interior and the exterior, and two points 
on the boundary are no longer necessarily connected by an arc. 

It is advantageous to use the polar angles (p = s/R. To be definite, we shall 
choose the angles always such, that adjacent points differ at most by tt. Simple 



84 



geometry tells that the positive angle a, as defined in (6.5),[^ obeys 

sin(a;) = Fdsin^ ^'^ 2^"^ ) ' (6.65) 
Moreover, we note the relation 

-Fj I sin((y9 — (/9o)| ^ sin(a) cos(a;) for Fd ^ 1 (6.66) 

which is needed in proving almost all the equations below. Finally, geometry 
tells that the normal components of the refiected velocities are given by 

-vjs^ no = ^ QF^sin((^ - (^o) ±sin(a)cos(a)^ , (6.67) 

for the short arc and long arc, respectively. They allow stating the prefactors of 
the map operators (6.20) and (6.21) explicitely in terms of the angle increment 

LP - LfQ. 



6.4-2 Operators for the integrahle map 

Upon choosing the symmetric gauge, X = 0? oiis finds that the actions of short 
and long arcs are merely functions of the difference of the initial and the final 
coordinate, 

as(v5 - '^q) ■= as(^V5; ^V^o) = - (a + sin(a) cos(a) + ^F^ sin(v? - v^o) 

(6.68) 

and likewise aL(v5 — v^o) '■= ^l{R'^] R'^o) = 1 — as(v'o ~ V')- This reflects the 
integrability of the classical motion. 

For the special case of the disk billiard the map operators (6.24) can be related 
directly to the magnetic generalization of the T operator [116]. Following [116] 
we define two operators, Ts and Tl, entirely in terms of the actions of a short 
and long arc, rj G {S,L}, respectively, with kernels 

(27ri)2 V dv^dv^o / 
Evaluating the mixed second derivatives of the actions, one finds that they 



15 



For the sake of clarity we use a (6.5) rather than a (6.55) in this section. 
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may be stated in a form 

d^as , . 1 (sin(a)cos(a) + ir|sin((^-99o))^ 

7^- (6-70) 



dif dipQ 27r sin(Q;) cos(a) 

dV 1 (sin(a) cos(a) - |rd^sin(y^ - (^p))^ 

-((^ - (^o) = (6-71) 



d(/7d(/7o 27r sin(Q;) cos(q;) 

which permits the direct comparison with equation (6.67). It follows that the 
operators P™* and Ps""* (cf eq (6.24)) are given essentially in terms of Ts: 

int/D D X W \^ /©(V'-'^o) ifrd>l 

Ps (i?<^, i?<^o) = ts (<^; <^o) p r . (6.72) 

-fill llid<l 

pf (fl^. %o) = t. i {^®(^» - ^) « > ; (6.73) 

They vanish whenever there is no classically allowed trajectory connecting the 
initial and the final point in the considered domain. Similarly, the operators 
Pl* and Pl^* are given as restrictions of Tl. 

int/D D N w ^ ^ /©(V-V'o) ifrd>l 

-K I U II i d < i 

pTiR'P, fl^c) = -t. (^; ^o) ^ If - ; (6.76) 

-ft II II id < 1 

Here we assume \ip — ipo] < n (as throughout this section). 



6.4-3 The explicit trace formula 

To obtain a semiclassical expression for the number counting function we start 
by calculating the kernel of the iVth power (P"it)^ at coinciding initial and 
final point sq. It is given by a (A^ — l)-dimensional integral, 

(p-^)"(^o,.o)=/nj(pr+PL")(^.,^.-i)]^^^ (6.76) 

with fixed sn = -so- This integral may be evaluated by the stationary phase 
method. For the same reason as above (Sect. 6.3.1) only the saddle points 
contribute which correspond to a physically allowed trajectory. However, they 
are now required to start and end at the point .sq. Each saddle point is char- 
acterized by the constant angular increment Aip, the jth component given 

by 

^j^^o+j^^, i = 0,... ,7V-1 . (6.77) 
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For given there is a finite number of possible increments 

M 



N 
int 



±27rf; M 



1, 2, ... M, 
1,2,...M, 



ifrd>i 

max if Td ^ !• 



(6.78) 



Here, the second index M has the meaning of a winding number .[^ It gives 
the number of times the trajectory encircles the origin. The maximum value 
is given by 



[arcsin(l/rd) N/ti] if > 1 
[N/2] ifrd<l, 



(6.79) 



where [•] indicates the integer part. The stationary phase approximation (A. 30) 
brings about a (A^ — l)-dimensional matrix of second derivatives. Its deter- 
minant is easily calculated since the difference between adjacent angles is 
constant: 



det 



difkdifi 



{s!\/\^)f-^ det 



(2 -1 
-1 ■•. 



A;,/=l...Af-l 







-1 2 



:iV(a"(A^)) 



N-\ 



(6.80) 



The number of negative eigenvalues is z/a = or z/a = — 1, respectively, for 
positive or negative sign of a"(A(y9) (ie, for long or short arcs). 

Taking the square-root of (6.80) cancels all but one of the prefactors in the 
integrand of eq (6.76). Altogether, the kernel (p™* + PlM (sq, Sq) is given by 



(27ri 



E 



d2(27rz/as(Av9)) 



+ 



d2(27rz/aL(A(/?)) 



^Af27rii/as(A</3)-if (Af-l) 



(6.81) 



d(^2 



^N2TTiuai^{Aip)-i^N 



e(rd - 1) 



It is a sum over all families of interior periodic orbits where each family is 
represented by the orbit starting at Sq. 



The nth power of the exterior operators, (Ps^* + Pif*) , assumes the same 
form except for the Heaviside function which appears in the short arc term 



We use capital letters for the indices A^, M in this section to avoid confusion with 
the radial and angular momentum quantum numbers, see (6.98). 
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of the sum. Naturally, the summation is now over the exterior periodic orbit 
families, the respective increments given by the set 



''i.H:!! !•!• ■!!■■■" '.^'.l (6.82) 



yN _ ) ^" at: 

±27rf; M = 1,2,...M^3. if Td < 1 . 



As the last step in forming the trace tr{(P™*)^} we have to integrate Sq. Since 
the expression (6.81) does not depend on sq this simply adds the factor 27iR/b. 

It follows that the fluctuating number function due to the skipping orbits 
(6.52) assumes the form 



nS«°" = (^)'I; E ( (6.83) 



TT 



N=2 Av>er^. 



int 

sin(A(/7) + sin(Q;) cos(q;) . /_ tt tt 



■sin 27rz/A^as(A9?) + N-^ 



(sin(Q;) cos(q;)) 2 V 2 4 

ir? sinfAoj) — sinfa) cosfo;) / . tt 7r\^, ,1 

+ ' ^ /. , , ^ ^ in (2nuNa^{A^) + N- - -) Q{T^ - 1) , 

(sm(Q;) cos(Q;j)2 V ^4/ J 

with a = arcsin(rd sin(| A(^|/2)). Analogously, the periodic orbit sum for the 
exterior problem is given by 

^ ^ ^ N=2 ^^^vg, y 

sin(Q;) cos(a) + sin( Av?) / tt 7r\ 
V , sm 2t:vN as(A(^) + + j e(rd - 1) 

(sm(aj cos(Q;jj 2 \ z 4/ 

sin(Q;) cos(q;) — ^r? sin( Aa?) . / ^ ^^tt 7r\ 1 

+ V , \,i (27r^AraL(Ay.) + N---\ . 

(sm(Q;) cos(q!))2 V z 4/ j 



r/ie conventional density of states 

The scmiclassical expression for the density of states is obtained by taking 
the derivative of the number function with respect to the energy. In order 
to compare with the result of Blaschke and Brack, which is in units of the 
conventional energy E, we have to take the derivative 
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Applying this differential to (6.83) yields the fluctuating part of the density 
for the interior problem 





cos \2niyN as (Ac/?) + N- + -j 

os(a) ( , . ^ TT 

cos \27ruN aL( A(/7) + N- 



n 




(sin(a) cos(q;)) 2 



4 



(6.86) 



This periodic orbit formula is identical to the result in [71]. It approximates 
the quantum spectrum of the interior magnetic disk only for weak fields Fd < 
1, when all trajectories are of the skipping type. For strong fields, Fd > 1, 
complete cyclotron orbits occur in the interior. One might wish to include the 
latter "by hand" into the periodic orbit sum. However, it was shown in [71] 
that energies close to the Landau levels cannot be reproduced this way. Rather 
than trying to refine the semiclassical approximation, we shall define a new 
spectral density of edge states below which will resolve the problem of the 
bulk contributions. 



6.5 The separable case 

We proceed to quantize the disk billiard for a second time - now using the 
separability of the quantum problem in a specific gauge. This way closed 
expressions for the spectral functions may be obtained, which yield explicit 

formulas for important quantities, such as the magnetization. Moreover, by 
formulating the relation to the periodic orbit formula derived in the prccccding 
section we can examine the effect of general boundary conditions on the trace 
formula. 



6.5.1 The disk billiard revisited 

The magnetic disk turns into a separable problem if we choose the symmetric 
gauge, X = 0, and place the center of the disk at the origin. In this case the 
canonical angular momentum L is conserved and the eigenstates are charac- 
terized by the quantum number 



In the second part of (6.87) we state the scaled angular momentum in terms 
of the radial distance c of the center of motion, cf (2.12). Along with the 
cyclotron radius p the latter determines whether the classical motion is of the 
skipping type. This is the case for it! — p < c < R + p. Hence, a quantum 
state (of energy u) corresponds to classically skipping motion if its angular 
momentum quantum number m is bounded from above and below by 

mm<« = + 2^/l}R (6.88) 

and 

mmin = max {r? — 2-\/i/-R, —v^ , (6.88a) 

respectively. Here, the scaled radius R := R/b enters as the only external 
parameter. 

We start with the traditional Bohr-Sommerfeld quantization method and pro- 
ceed to discuss its relation to the periodic orbit formula of Sect. 6.4. The exact 
quantization in terms of special functions is discussed in Appendix A. 3. 



6.5.2 Semiclassical quantization 

Using polar coordinates (r, •&) , the ansatz 

^(^,^)^^^e^-'^ (6.89) 
Jr/h 

transforms equation (4.1) into the form of a one-dimensional Schrodinger equa- 
tion for the radial function ip{f). 

-lV"{r) + (^-^ ^ - uj if if) = (6.90) 

It may be solved to leading order in using the standard WKB technique, 
see eg [126,34]. 



The semiclassical wave function 

It follows that in the energetically allowed region the resulting semiclassical 
wave function has the form 

^^-\r, ^) = A/-,,, '-'i^^^'^^^^-p . (6.91) 



90 



Here, the phases ^^^isk ^d?sk obtained by an integration starting at the 
interior and exterior classical turning points of the radial motion, respectively. 



,^mt I ~\ 1 nr~^^ 7^0 ^ I , . / 2z/ + m-r^ 

<P(jig]^(z/, m, r) = -JAur^ — [r'^ — m)'' — {u -\ ) arctan 

2 2 



( 



yjAur"^ — (r^ — m) 



((2z/ + m)r'^ — m? \ vr / m — Iml 
, ^ — I +- + ^ 

raJAur'^ — (f^ — m)^^ 



m 

— arctan , 

2 \ Tn . Aur"^ — ('■r2 — m\'^ / 2 \ 2 



(6.92) 



and 



$^^\(z/, m, r~) = TT h/ + ^ - ^^^^{u, m, r~) . (6.93) 



As for the normalization factor Adisk) we find [126] 



-2 vr 



/ 2u + m- 



(ATdisk)"' := X ^ 2 I , — I , (6.94) 



where again the upper sign stands for the interior problem. 



A spectral function 

Allowing for general boundary conditions (3.6) at the disk radius r = R, we 
obtain the quantization condition 



±cot($:,f/r(^,m,i?)-^) 



3 



i±A) {AuR' - (R' - m^y 

(6.95) 



2^r{AuR^ - (i?2 - mf) + {±K)R^(2u + m - R'^ 

The boundary condition enters on the right side through the dimensionlessP^ 
mixing parameter A (3.7) which vanishes for Dirichlet boundary conditions. 
In order to transform the dependence on the boundary condition into a phase 



The dimensionless mixing parameter (3.7) is introduced for convenience. Strictly, 
it is not an independent variable but should be replaced by 2-^/7' A/6 everywhere (to 
avoid energy dependent boundary conditions). This distinction does not matter, 
ultimately, since we are only interested in the derivatives at A = 0, see (6.99), (7.3). 
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shift CCA, we define 



/ 



a\{i','m, R) :— arctan 



(±A) (4i/i?2 - _ 



2^r(AvE? - - m)2) + (±A) R^{2p + m - R^ 

(6.96) 



The semiclassical quantization condition (6.95) is then readily brought into a 
form, 



cos ( ^d'isk^''*!^) T OiA{i', m, R) 



TT 



0, 



(6.97) 



which permits a spectral function ^ to be written in terms of two quantum 
numbers, the number of radial nodes n, and the angular momentum m. 



■- $ 



int / ext 
disk 



R 



i^,m, — ] To^A[u,m,— 



R 



n + - TT 



(6.98) 



with n e No,mniin ^ tti < (6.88). Its zero in u yields the semiclassical 

energy of a state with given radial and angular quantum numbers n and m. 
Although the energies are defined implicitly by (6.98), the spectral function 
yields explicit formulas for the infinitesimal change of the energies as an ex- 
ternal parameter is varied. For the derivative of the energy with respect to 
the boundary mixing parameter at Dirichlet boundary conditions (A = 0) we 
obtain 



di/ 
dA 



\ 



1 - 



(i?2 - m)2 
AvR? 



A=0 



TT I 2u + 111 - n- 

— =F arctan _ _ =— 



(6.99) 



This short formula is further compressed below and needed soon. 



The magnetic moment 

It was shown in Section 3.4 that the scaled magnetic moment of a quantum 
state in the magnetic billiard is determined by the derivative of its energy 
with respect to the magnetic length, see (3.25). Prom the semiclassical spectral 
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function (6.98) we find 



(,2 d (sc) 

j^s,disk — ip arctan 



d62 d ^(sc) 2 / 2;/ + /;/ - ^ ' 



2 V^4i/i?2 -(R2-rny^ 



Alternatively, the expectation value may be calculated directly, using the semi- 
classical wave function (6.91). We obtain indeed, after lengthy transforma- 
tions, 

i{^(==>|(f X v).y„|l^"=') = /im - if')^'")] Mfdf 

2 TT / 2z/ + m - i?2 



arctan 



2 y^J^vR^ - {R'^-mf 



in agreement with (3.25). In the above radial integration (which is limited 
by the disk radius and the interior or exterior turning point, respectively), 
the strongly fluctuating cos^-term was replaced by its mean. The fact that the 
exact relation (3.25) is reproduced shows that this approximation is consistent 
with the semiclassical one. 



The hulk states 

States with angular momenta beyond the bounds given by (6.88) are not in- 
cluded in the spectral function (6.98). Classically, they correspond to cyclotron 
motion. The semiclassical energies of these bulk states are determined by the 
condition that the two wave functions (6.91) defined from the interior and 
exterior turning points must match. They are given by the Landau energies 
v — n -\- \, and the wave functions are readily shown to exhibit a magnetic 
moment of — i^. 

However, the exact quantization, which is given in Appendix A. 3, does not 
distinguish between edge and bulk states and the bulk energies exhibit devia- 
tions from the Landau energies. An asymptotic treatment of these exponential 
corrections to the bulk energy is given in Sect. 7.3.1. 



6.5.3 Relation to the periodic orbit formula 

In the preceeding section 6.5.2 the semiclassical quantization was carried out 
according to the traditional Bohr- Sommerf eld rule for separable systems. It 
is based on the quantizing tori, ie, those invariant manifolds in phase space 
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whose scaled actions are integers. This should be contrasted to the periodic 
orbit formula for the magnetic disk derived in Sect. 6.4. The latter is a sum 
over the rational tori, whose classical frequencies are commensurate [113]. In 
order to sketch how the trace formula is connected to the Bohr- Sommerf eld 
quantization we follow the work of Berry and Tabor [112] who derived the 
trace formula for general integrable systems. In particular, this permits us to 
show how the trace formula is modified if one allows for general boundary 
conditions (3.6). 

The semiclassical spectrum is given by the energies 1/(71, m) , which are implic- 
itly defined as the roots of ^^jg^ (6.98). We may write the spectral density as 
a sum over the two quantum numbers, 

= E '^('^0 - i^{n, m)) 

n,m 

00 . 

= d(i/o) + E / e''^'(^"+^"*)5(i/(n, m) - vo) dn dm 

N,M=-oo 

^d{uo)+ E / ^ e2-W-o,m)+Mm)^^ ^ (g_;^Q2) 



N,M=— 00 ' 



where the Poisson summation formula (eg [127]) was employed to transform 
the sum into an integral. (Boundary corrections which are to higher order in 
u are neglected.) The sum excludes the term with N = M = 0, which yields 
the mean density d. Upon integrating n the 5-function selects the real valued 
"number" of radial nodes which is known explicitely from above. 

n{uo, ^) = ^ (^di!k^''*(^o, m, R) T aA{iy, m, R) - (6.103) 

We evaluate the remaining integral in (6.102) by the stationary phase approx- 
imation. The phase shift cka should be neglected in the saddle point condition 



dm 



*dfs^*TH = 27r^^A(^ (6.104) 



M 



since $ is of order v (while cka is of order 1). A detailed calculation shows, that 
the angles A(/? selected by (6.104) are indeed given by the sets 'Pj^t/ext defined 
in (6.78) and (6.82) (modulo 2ti). It is now convenient to characterize the 
corresponding skipping trajectories by the signed sine of the angle of incidence 

„2 _ d2 _ 2 

e := ^ = n X V , (6.105) 

2Rp ' ^ ^ 

such that the former quantum number m is given by the real value 

m = ^^ + 2v^^£. (6.106) 
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One finds, after a lengthy calculation, that 

2*5^''* + ;^27rm ^7ru + 2u arcsin((7) - 2iyaVl - + B? sin( A(^) 

= 27ri/a (6.107) 

with a = ^{sR + p)/c defined in (6.55), and "a" the geometric action (6.56) 
of one arc. Transforming the summation in (6.102) to positive N we obtain, 
observing (A. 29), 

9 1 _d_(F,int/ext 



AfgN.MgZ: 

n 71 



^int/ext 



X COS (N2niya + N- + - sgn {NdX^l"''') T '^Na^ 



We note the derivatives 



d' 



2 



C:(^" = -r^^T^ (6.109) 



dm2 2i?cv^r^2 

and 

^Cr* = I + arcsin(a) = \ ^ (2vrz.a) . (6.110) 

The last equality permits to integrate the spectral density immediately. It 
yields the oscillatory part of the number counting function, 

^o) = (-) E E ^-^ (6-111) 



X sin i^-KuaN + ^A^ + ^ sgn(cr) =F 2iVaA^ , 

which may be compared to the trace formulas (6.83) and (6.84) obtained from 
the boundary integral equations. The agreement of the prefactors follows after 
a tedious discrimination of the various cases (interior/exterior, short/long arcs, 
and R ^ p.) As the only diff'erence compared to the Dirichlet result of Section 
6.4 we observe the non-vanishing phase factor ^2Na\ for finite A. 



The effect of general boundary conditions 

This result suggests that, compared to Dirichlet boundary conditions, the only 
effect of a finite mixing parameter is the appearance of an additional phase 
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shift at every point of reflection, 



^2a\ = —2 arctan 



A. 



\ 



\ 



1±A 



1 p/R + e 



-2 arctan (a Vl - + 0{K^) , 



as A ^ 0. 



(6.112) 
(6.112a) 



Here, we stated (6.96) in terms of the geometry of the periodic orbit (6.105) 
and of u. One might be tempted to "generahze" the result (6.112) to arbitrarily 
shaped billiards, by replacing the disk radius R by the radius of curvature at 
the point of reflection. However, the phase shift at a point of zero curvature 
(which is given in appendix A. 6) is not reproduced correctly this way. Only 
the limiting expression (6.112a) for small A matches with its zero curvature 
analogue. The latter is determined merely by the (unsigned) angle of incidence 
with respect to the normal at the point of reflection. 



n V 



(6.113) 



Its form coincides with the non-magnetic result [76] . This generality suggests 
that at small A any billiard exhibits the additional phase (6.112a) at the points 
of reflection. All what will be needed below is this dependence to flrst order in 
A. It shows up in the derivative (6.99) which we can now rewrite in terms of 
the geometric quantities a and s, see (6.55) and (6.105), describing the length 
of the arc and the angle of incidence, respectively. This way the dependence 
of the energy on A assumes a particularly simple form. 



dA 



A=0 



I + arcsin((T) 



(6.114) 



7 A spectral measure for edge states 



In the following, two different quantitative deflnitions for edge states are in- 
troduced. We discuss their relation, the asymptotically smooth form of the 
corresponding spectral densities and their semiclassical interpretation. 



7.1 Bulk states and edge states 



In Chapter 3 we discussed the existence of two types of states in the spectra 
of magnetic billiards. The wave functions of a few typical representatives were 
given in Chapter 5, in Figs. 5.5 - 5.8. Observing these images one might think 
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that it is an easy task to separate the spectrum into two disjunct sets, edge 
states and bulk states, respectively — similar to the classical trajectories which 
are either skipping or cyclotron orbits. However, in general there is no way to 
perform such a strict partitioning. Rather, a general wave function may share 
some features of both, edge states and bulk states, to a certain degree and 
there is a gradual transition taking place between the characteristics of the 
two types of states. 

On the other hand, there is a clear need for an objective way to separate 
the edge from the bulk contributions in the spectrum. Bulk states are very 
uninteresting. They do not contribute to transport and tend to accumulate 
in the vicinity of Landau levels. Moreover, their number often dominates the 
spectrum. In the exterior a mean density of states cannot be defined (as a 
derivative of a mean counting function) due to the infinite number of bulk 
states showing up in the vicinity of each Landau level. As a consequence, the 
oscillatory part of the spectral density cannot be extracted — which seems to 
impede any statistical or semiclassical analysis of the exterior problem. Also in 
the interior the accumulation of bulk states severely complicates the analysis 
of the spectrum in terms of the classical skipping motion. 

To the best of our knowledge, no general and objective criterion for what 
constitutes an edge state has been proposed so far.[^ Clearly, any reasonable 
definition must take into account the fact that there exist transitional states 
between pure edge and pure bulk states [128,48]. We illustrate this gradual 
transition by referring to Figure 5.11 on page 62. It displays an exterior spec- 
trum as a function of the boundary mixing parameter A (3.7). In the level 
diagram one observes that the infinitely many states which accumulate near 
the Landau levels are hardly affected by changes of the boundary condition. 
These are bulk states. The extreme insensitivity of their energies with respect 
to A is explained by the fact that bulk wave functions are not localized at the 
the boundary. They approach it with an exponentially damped tail giving rise 
to exponentially small energy shifts. (This will be discussed quantitatively in 
Sect. 7.3.1.) Other states depend strongly on A because they are localized at 
the boundary. They are naturally associated with edge states. The fact that 
states may have a transitional nature can now be seen in the right part of 
Figure 5.11. One observes a sequence of bulk states which originate from the 
Landau level and gradually turn into edge states with a strong dependence on 
the boundary, ie, a large slope. Obviously, it would be inappropriate to split 
this sequence at an arbitrary point into two distinct parts. 



In [77] Akkermans et al. propose special "chiral" boundary conditions. In the case 
of a disk billiard (and only there) this leads to a gap in the level diagram which 
may be interpreted as a separation into edge and bulk states. A disadvantage of 
this approach is that the obtained spectrum is unrelated to the standard Dirichlet 
spectrum. 
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In order to quantify the notion of edge states we propose to attribute positive, 
real valued number w„ > to each eigenstate ipn which gives a measure for 
the degree to which the state has the character of an edge state. This way a 
density of edge states can be defined which applies in the interior as well as in 
the exterior and which consistently accounts for the gradual transition from 
edge to bulk. As compared to the standard density (3.9) each 5-contribution 
is weighted individually in our definition: 

oo 

dedgei^) ■= ^ ^ri) ■ (7.1) 

n=l 

Hence, the sum still extends over all states in the spectrum, but for a proper 
choice of the quantum weights Wn the bulk states are effectively suppressed by 
their small values. There are a few requirements which are naturally imposed 
on the definition of the quantum weights Wn'- 

(i) The mean density of edge states dedge must be well defined in the exterior. 

(ii) The interior and the exterior mean densities should be equal to leading 
order. 

(iii) As a sequence of bulk states approaches a Landau level their weights 
must decrease at least exponentially. 

(iv) Last not least: The weights should admit a semiclassical interpretation 
which complies with the intuitive notion of edge states. 

To make the last requirement more specific consider the semiclassical periodic 
orbit formula for the oscillatory part of the density of edge states. It should be 
a sum over the interior or exterior skipping periodic orbits which differs from 
the expression (6.61) for the standard density at most by a classical weight 
Wj attributed to each periodic orbit contribution 7. 

dZU^) = - E ^ r cos (2nuAi^) - nn, - ^^,) . (7.2) 

^ r^|trM(7) -2|' ^ ^ ^ 

Similar to the quantum weights w„ which must consistently fade out the bulk 
contributions, the classical weights should vanish gradually for periodic 
orbits which are increasingly close to being detached from the boundary.]^ 

In the following sections we introduce two different definitions of the quantum 
weights Wn satisfying the abovementioned requirements. The first one, which 
is very convenient from a mathematical point of view, is discussed in Sect. 7.2. 
It has the property that it renders the leading term of the mean edge density 



Equation (7.2) is stated for completely chaotic classical dynamics. In the case of 
integrable dynamics the corresponding periodic orbit sum is modified analogously 
by the same classical weight. 
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dedge proportional to the circumference JSf of the biUiard rather than the area 

The second definition, given in Sect. 7.4, is the most natural choice from a 
physical point of view. It has the property that the interior weights assume 
unit value for large cyclotron radii, ic, it approaches the standard density if 
bulk states cannot exist in the interior. Consequently, the interior mean edge 
state density equals the standard mean density for this second definition. The 
relation between the two different definitions of the quantum weight is also 
discussed below. 



7.2 A spectral density based on the boundary conditions 



Figure 5.11 suggests that the slope in the level diagram provides a quantitative 
criterion for the degree to which a state is of the edge type. We therefore 
propose to weight each Dirichlet eigenstate \ipn) by the derivative of its energy 
Un with respect to the boundary mixing parameter at Dirichlet boundary 
conditions A = 0, 



di/, 



dA 



A=o 2^/u dA 



(7.3) 

A=0 



which is positive valued. It will be shown in the sequel that this definition 
complies with the requirements for a definition of edge states stated above. In 
particular it admits a semiclassically meaningful interpretation, as discussed 
in Sect. 7.2.2. 

For the mean density we obtain the simple expression (cf Sect. 7.3.2) 

=Sf 1 1 

4dgeH = ^'^^ T 2 ' (7-4) 

where the upper sign stands for the interior problem. The leading order term is 
proportional to the circumference of the billiard. As argued below, the sec- 
ond order term may be related to the mean curvature of the billiard boundary 
(which is positive from the interior and negative from the exterior) . 

Before discussing the various asymptotic properties which come along with 
the definition (7.3) wc present a few examples of edge spectra of interior and 
exterior billiards. They provide a first indication that the quantum weights 
succeed to sort out the bulk contributions consistently. A more quantitative 
check of this assertion will then be given in Chapter 8 where we perform a 
statistical analysis of the edge spectra. 

An edge spectrum {(i^njWn)} consists of the energies and the attributed 
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Fig. 7.1. Weighted spectra of the interior (top) and the exterior (bottom) magnetic 
disk (area = vr, 6 = 0.1). Each dot (open and filled) corresponds to an eigenstate 
IV'n) with the energy given by the abscissa. The ordinate indicates the attributed 
quantum weight Wn defined in eq (7.3). It serves to distinguish edge states (with large 
Wn) from bulk states. The latter accumulate at the Landau levels v = N+^, G N, 
and are characterized by vanishingly small weights Wn- A sequence of transitional 
states emanates from each Landau level and connects with the edge states. As an 
alternative criterion, the angular momentum quantum number permits to decide 
whether the state corresponds classically to skipping motion (full dots) or cyclotron 
motion (open dots), see text. 
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weights Wn. Figure 7.1 gives an example of an interior (top) and an exterior 
(bottom) edge spectrum at a strong magnetic field. The spectra belong to a 
disk billiard of unit radius and were obtained from eqs (A. 23) and (A. 25). 
In these plots each point belongs to one eigenstate and indicates the weights 
versus the energy. One observes how the weights segregate edge states with 
large Wn from the bulk states. The latter accumulate at the Landau levels 
t/ = iV + |, G No, with vanishingly small weights. They are highlighted 
in Fig. 7.2 which shows the same data as Fig. 7.1 on a logarithmic scale. A 
sequence of bulk states can be seen emanating from each Landau level and 
gradually turning into edge states. 

Since the disk billiard has a second quantum number we can compare our 
characterization of edge states based on weights with a classical criterion. As 
discussed in Sect. 6.5.1 a state corresponds classically to skipping motion if the 
angular momentum quantum number lies within the bounds given by (6.88). 
In the Figs. 7.1 and 7.2 we indicate those states with constants of the motion 
which belong to a skipping trajectory by a full dot. The others are given by 
a large open dot. One observes that the effective separation produced by the 
weights complies with the classical criterion. At the same time it seems more 
appropriate to formulate the separation in terms of a continuous quantity. 
This is the more so since a second quantum number does not exist for shapes 
other than the disk. 

Figure 7.3 shows the exterior edge spectrum of an ellipse billiard (which is 
not integrable in the magnetic field) . We took the same area and magnetic 
length h as for the disk in Fig. 7.1. Comparing the ellipse spectrum to the 
disk one observes that they resemble in their gross features. In particular, the 
bulk states behave very similarly. However, for the ellipse there are additional 
structures showing up in the distribution of the weights of edge states. These 
can be related to features of the classical (mixed chaotic) phase space, as will 
be shown below. 



7.2.1 Edge state counting functions 

Upon integrating the density one obtains the edge state counting function 

Nedge(z^) / 4dge(i^') di^' = X) w„ 0(1/ - i/„) , (7.5) 
-'^ n=l 

which is a weighted staircase. It jumps by Wn at the corresponding spectral 
point z/„. Again, the sum formally includes the bulk states. We expect their 
contribution to be effectively eliminated by the rapid decay of the weights such 
that the edge state counting function should bear no marks of the Landau 
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Fig. 7.2. Weighted spectra of (a) the interior and (b) the exterior magnetic disk in a 
semilogarithmic plot (same data as in Fig. 7.1). The logarithmic scale highlights the 
states which correspond classically to cyclotron motion (open dots). One observes 
that the the respective weights decrease exponentially. 
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Fig. 7.3. Weighted spectrum of the exterior eUipse billiard (with eccentricity 0.8, 
area = tt, magnetic length 5 = 0.1). It should be compared to the exterior disk, 
Fig. 7.1 (bottom). While the bulk states are very similar, one observes that the edge 
weights no longer lie on smooth curves but tend to cluster. These structures can be 
related to the classical (mixed chaotic) phase space. 
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Fig. 7.4. Spectral counting functions for the magnetic disk (R/b = 10). (a) Total 
number of interior states. (b,c) Number of (b) exterior and (c) interior skipping 
states according to the angular momentum criterion. (d,e) Weighted number of 
edge states for the (d) exterior and (e) interior problem. When averaged the curves 
are well reproduced by the smooth counting functions, (3.10), (3.17), (3.16), and 
(7.6), respectively (not shown). The inset gives the counting functions for the first 
four Landau levels. The small kinks seen in Ns^jp and Nedge are damped away at 
higher energies. 

levels. According to (7.4) its smooth part is given by 

Nedge(^) = ^^^^ T +0(1). (7.6) 

Note that the leading order exhibits the same functional dependence as the 
phase space estimate of the skipping states for the periodic straight line prob- 
lem (3.18). The only difference is an additional prefactor of |. 

In Figure 7.4 we compare various spectral counting functions obtained from 
the magnetic disk spectra given in Fig. 7.1. Curve (a) shows the total num- 
ber of states in the interior. It exhibits distinct steps at the energies of the 
Landau levels. In the exterior a total counting function does not exist but the 
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Fig. 7.5. Solid line: Oscillatory part (7.7) of the edge state counting function for the 
exterior magnetic disk (Kjh = 10). (The unspecified constant part of Ncdge was not 
subtracted.) Dotted line: The oscillatory part convoluted by a Gaussian of width 
(7 = 0.1 minus an offset of 3.0. 

angular momentum criterion (6.88) enables counting the exterior states of the 
skipping type, see curve (b). The corresponding number of interior skipping 
states is indicated by curve (c). As one expects these two counting functions 
hardly exhibit steps at the Landau levels but they show a different functional 
dependence (given by cqs (3.17), (3.16)). In contrast, the weighted exterior 
and interior edge state counting functions, curves (d) and (e), respectively, 
display the same mean values to leading order. Their average is reproduced 
by eq (7.6) and no marks of the Landau levels are visible. 

To examine more closely the suppression of the bulk states we plot the fluc- 
tuating part of the edge state counting function, 

NS^ge(^) = Nedge(^) - Nedge(^) ■ (7-7) 

Figure 7.5 shows this quantity as obtained from the exterior spectrum. One 
observes that the fluctuating function hardly exhibits a signature of the Lan- 
dau levels. The dotted line depicts N°^gg(i/) after convolution with a narrow 
Gaussian which smoothes out the fluctuations. Its oscillations are due to the 
remnant contributions of the bulk states which are slowly damped out for 
higher energies. The fluctuating part of the exterior ellipse spectrum shown 
in Fig. 7.3 looks very similar [106]. 
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1.2.2 The semiclassical density of edge states 



The standard spectral density (3.9) is given by tlie derivative of tiie number 
counting function N(i/) with respect to energy. The edge state density (7.1) 
may also be formally defined as a derivative, now with respect to the boundary 
mixing parameter A, at Dirichlet boundary conditions (A = 0), 



dA 



^ _ b dN(i/) 

A=0 ^ " 



(7.8) 

A=0 



Hence, the periodic orbit expression for the oscillatory part of the edge state 
density can be deduced immediately once we have the semiclassical formula 
for Nose at hand. For the time being, we restrict ourselves to hyperbolic sys- 
tems. Combining the results of the previous chapter (cqs (6.54), (6.112a), and 
(6.113)), the number of states based on the skipping part of phase space is 
given by 

NSi'^(^; A) = r (7.9) 

^ T r^|trM(7) -2 ' 

X sin (2nuA{-f) - -/i^ - 2A^\hjVj\] + 0{A'^) . 
V 2 .^^ ) 

Compared to (6.54) the leading order dependence on A is included, as dis- 
cussed in Sect. 6.5.3. (See Tables 6.1 and 6.2 for the definition of the various 
quantities in (7.9).) Since the semiclassical bulk states do not depend on the 
boundary condition their contribution vanishes when taking the derivative. 
Using (7.8) one obtains the semiclassical trace formula for the edge state den- 
sity at Dirichlet boundary conditions, 

^:^ge(^) = ! E ^^^'^""'""'^ , COS (27r^A(7) - vrn, - ^/.,) . (7.10) 



^ T r^|trM(7) 



This expression should be compared to that of the unweighted density of states 
(6.61) which was obtained by taking the derivative with respect to the energy 
V. It exhibits the scaled time of flight r-y (6.60) in the numerator. Hence, the 
periodic orbit sum (7.10) differs from the standard spectral density only by 
the prefactors 

^^L^ . (7.11) 

They attribute an individual classical weight to each skipping periodic orbit 
7. The classical weights are given by the time averaged value for the normal 
component of the velocity |nv| at the points of reflection. They vanish for 
cyclotron orbits. Similar to the quantum weights, the lead to a gradual 
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transition from edge to bulk contributions. It is easy to see that in the hmit 
of a "grazing" trajectory of increasingly many short arcs variations in the 
curvature of the boundary may be neglected and the classical weights 
approach a constant value. In the opposite case of an orbit which is almost 
detached from the boundary the weights vanish since the normal components 
of the velocities approach zero at a finite time of flight in the denominator of 
(7.11). 

It is instructive to compare the distributions of quantum and classical weights. 
A direct comparison is not possible since the classical and the quantum weights 
are associated with different objects, eigenvalues and periodic orbits, respec- 
tively. In Fig. 7.6 we compare the distribution of classical weights to the cor- 
responding weighted quantum spectrum. The data were obtained for the inte- 
rior elliptic billiard, and are given in both cases as a function of the classical 
cyclotron radius p. The shade in the distribution of classical weights gives 
the probability for obtaining a certain weight if the trajectories are chosen 
randomly with respect to the invariant measure. It was approximated numer- 
ically by the histogram over a finite number of trajectories taken uniformly 
from phase space .[^ Remarkably, one observes that the characteristic features 
of both distributions coincide. This shows that the quantum weights may be 
considered the expectation values of an observable which has a classical limit, 
ie, they measure a classical property. This holds in spite of the fact that the 
Wn are defined in terms of the boundary condition, which has no classical 
analogue. 

The bifurcating structures seen in Fig. 7.6 are due to stable periodic orbits. At 
the bifurcation points periodic orbits 7^"^ with a fixed number of reflections 
exhibit the smallest possible cyclotron radius (n^ = 6 in the case of the 
rightmost structure). As the cyclotron radius increases, the orbits turn into 
pairs with either longer or shorter arcs. (Some of the corresponding islands of 
stability in phase may be identified in Fig. 3.1, left column.) 

7.3 Asymptotic properties of edge and bulk states 

We proceed to briefly discuss the leading order behavior of the bulk energies 
and the smooth part of the edge counting function. Both estimates are ob- 
tained in the semiclassical limit such that it is legitimate to substitute the 
boundary by the straight line with periodic boundary conditions discussed 
in Appendix A. 6. The finite mean curvature and variations of the boundary 
curvature are expected to appear only as higher order corrections, see the 
discussion in Sect. 3.3.3. 

In the numerical calculation we could use general (non-periodic) trajectories to 
approximate the periodic orbits which are dense in phase space. 
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Fig. 7.6. Weighted edge spectrum (7.1) (top) and phase space distribution of the 
classical weights (7.11) (bottom) for the interior ellipse. To ease comparison, also 
the quantum spectrum (calculated at constant 6 = 0.1, as in Fig. 7.3) is given in 
terms of the classical cyclotron radius { p = bx y^.) One observes that the quantum 
weights tend to mimic the structures in the distribution of classical weights (which 
are due to stable islands in phase space, cf Fig. 3.1). [figure quality reduced] 
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Fig. 7.7. Bulk state energy shifts Ai/„, and weights Wn, for the magnetic disk 
{R/b = 15.0111) at the 50th Landau level (on a double- logarithmic scale). The 
interior (• ) and exterior (O ) bulks states are approximated by the periodic line 
(x), cf (7.12) and (7.13). 



7.3.1 Bulk state energies and weights 



The energy shift of a bulk state for general boundary condition is derived in 
Appendix A. 6 by a uniform approximation, see (A. 51). An asymptotic expan- 
sion, which amounts to the WKB approximation in the energetically forbidden 
region, yields the expression 



Aum ~ ^ exp ( - 2u[q^^ql-l - log (g„ + .Jql-l))) (7.12 



27r 

\ 1 1 exp { N -\ 2 ( " ) ) as m — > cxo 



1 / 27ibm \ ^^^^ / 1 / nbm 



(7.12a) 



with G N the Landau level and g™. the quantized distance of the cyclotron 
center to the boundary, as defined in (A. 49). We observe that the bulk state 
energies approach the Landau levels exponentially fast (indeed, like a Gaus- 
sian) as the integer m increases, ie, as the distance of the cyclotron motion 
from the boundary grows. 



The weights of bulk states follow likewise, by taking the derivative of equation 
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(A. 51). Essentially, they decay as fast as the shifts of the bulk energies: 



1 

\ 2 



In Figure 7.7 we show exact bulk state energies and weights in a double- 
logarithmic representation. They belong to the interior and the exterior disk 
at R/b = 15.0111 and to the 50th Landau level. The crosses indicate the zero 
curvature estimates according to eqs (7.12) and (7.13). One observes that the 
asymptotic weights and the spacing between the asymptotic energies match 
approximately the exact values, and lie between those of the finite curvature 
case. The approximation is improved when both m and N increase. 



7.3.2 The mean edge counting function 

As a second application, the periodic straight line problem allows the straight- 
forward derivation of the leading term of the mean edge state counting function 
(7.6). We simply identify the transverse quantum number n as a partial count- 
ing function for states with fixed quantum number m. An explicit formula for 
n, which includes the dependence on A follows from ^[inl — (A. 40). The sum 

Nskip(i^,A) = ^n(m,A) (7.14) 

m 

yields the total number of states corresponding to skipping motion. Taking 
the derivative with respect to A and replacing the summation by an integral 
we obtain the leading order of the smooth edge state density, 

fi_S)':^+0(l) 

J-^ V u J nb 

(7.15) 

The same result can be derived from the magnetic disk problem discussed in 
Sect. 6.5.2. It follows from (6.97) that the semiclassical spectrum is obtained 
by requiring that the phase 

V^'/^'^iu, R, A) = rn, R) T «a(^, m, ^) - ^ (7.16) 

is an integer multiple of tt. Hence, the smooth counting function of the skipping 
states is obtained by 

" m 



dedgei^) 



dNskip(i^;A) 



dA 



A=0 



2^ 



+0(1) . 
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where the sum can be replaced by an integral over the interval |m — < 
2v^R. Using the (6.96) one gets 



i^int/cxt 



dA 



^ 1 

-y1 



A=0 



27r6 



(7.18) 



which reproduces (7.15). This is further evidence that the periodic line problem 
yields the leading order terms consistently. The second order term in (7.4) is 
not obtained this way. It will be deduced in the next section by relating the 
quantum weights of the disk to the magnetic moments of the states. 



l.J^ Edge magnetization as a spectral measure 



The preceeding section showed that our first definition (7.3) of the quantum 
weights yields an efficient and mathematically natural way to separate edge 
from bulk. However, the physical interpretation of the mixed boimdary con- 
ditions is not immediate [76]. We therefore propose an alternative definition 
of the weights which is physically more accessible. It is obtained from the 
expectation value of the magnetic moment of the state. 

In Section 3.4 we introduced the edge magnetization (3.28) of an interior bil- 
liard. It gives the scaled excess magnetization that is induced by the presence 
of the billiard boundary and was defined as 

Medge(^) = E - ^n) ■ (7.19) 

n=l 

Like the edge state counting function (7.5) this is a weighted staircase. The size 

of the steps are now given by a derivative with respect to the magnetic length 
rather than A. Since the Landau levels do not depend on b the bulk states 
contribute merely to a negligible degree to (7.19), as discussed in Sect. 3.4. 
Hence, it is reasonable to extend the definition (7.19) of the edge magnetization 
to the exterior problem. 

The exterior edge current shows an orientation which is opposite to the interior 
(see Sect. 5.2). Therefore, one expects A^edge to turn negative in the exterior 
- which is indeed found. Moreover, the mean is finite in both cases and given 

by 

-r— . . 1 Q 1 j5f 3 ,_ 

The interior case (upper sign) follows from eq (3.32) while the exterior one 
(lower sign) is suggested by symmetry and confirmed empirically. Like in the 
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case of the edge counting function (7.5) the moduU of the mean interior and 
exterior edge magnetizations are equal to leading order. This suggests to use 
the edge magnetization density 

niedge('^) = ^ Medgeil^) = ^ ^ (^ ~ ^ri) , (7.21) 

which was introduced in Sect. 3.4, to define a physically motivated spectral 
measure for the edge states: 

1 oo 

Ci(^) ±- medge(^) = E S{u - U^) (7.22) 

^ n=l 

The index (M) is used to distinguish this magnetization based density of edge 
states from the former definition (7.3). The weights are now given by 

:= ±^^^ = ^mSM^ . (7.23) 

Again, these positive quantities may be obtained as derivatives of the ener- 
gies with respect to an external parameter (which is the magnetic length h 
in the present case). At the same time they are expressed in terms of the 
(symmetrized) expectation values of the scaled magnetic moment (3.25). The 
weights vanish as the Landau levels are approached since for interior and exte- 
rior bulk states the scaled magnetic moment approaches the diamagnetic value 
('V'nl^^l^n) ~^ ~^ from above and below, respectively. As for the smooth edge 
state density, the expression is obtained immediately from (7.20): 

It coincides with the standard mean density for the interior case. Moreover, 
the interior weights approach unity as the magnetic field B is decreased and 
bulk states no longer exist in the interior. This is seen immediately if we write 
the weight (7.23) in terms of conventional units 

^.r = ±(||(^n|rxv|V'„) + l) . (7.25) 

In the limit the interior weights assume unit value and the interior 

edge state density turns into the standard density. This limit does not make 
sense in the exterior since in this case the spectrum ceases to be discrete; at 
finite field the exterior weights remain always non-negative. 

The expression (7.25) shows that the weighted spectrum may be obtained 
immediately by measuring the magnetic moments of the states. The scaled 
magnetization density 3.26 follows from the conventional magnetization den- 
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sity (3.21) by a simple multiplication with the magnetic field: 



m{u,b')^Bm{E^%,B^^) (7.26) 

This relation to an experimentally measurable quantity is a clear advantage 
of the present definition of the edge state density. It is achieved at a price 
— the leading order of the mean edge density is now determined by the area 
of the billiard rather than by its circumference, see (7.24). It indicates that 
with this measure the quasi one-dimensional character of the edge states is not 
accounted for to the same degree as by the former definition of (iedge- However, 
it does an equally good job in consistently suppressing the bulk contributions. 
Moreover, in the case of a disk both spectral densities, dedge and rf^e; 
identical up to a factor. This is seen from (A. 27) which leads to the equation 

. (7.27) 

This relation is as surprising as fortuitous and does not hold for general billiard 
shapes. Nonetheless, it allows to deduce the second, constant term of the mean 
density (7.4) by comparison with the smooth edge state density (7.24). Being 
the next order after the circumference term it is determined by the mean 
curvature which is equal for all simply connected boundaries, /p n{s)ds = 
±27r (according to the Gauss-Bonnet theorem). 

It follows from (7.27) that for the disk the magnetization based edge spectra 
differ from the spectra shown in Fig. 7.1 merely by a geometric transformation 
and there is no need to reproduce them here. As for the ellipse. Figure 7.8 
shows the weights as obtained from the edge magnetization. Like in Fig. 7.6 
the structures in the distribution are reproduced by the probability density of 
the corresponding classical weights. 



Semiclassical edge magnetization 

The semiclassical periodic orbit formula for the complete magnetization den- 
sity is given in (6.63). Likewise, one finds that ±iri°|gg/z/ has the form of the 
trace formula for the standard density with each periodic orbit now weighted 
individually by 



w, 



7 



27rz/A(7) 



a. 



fj+i X 
2p2 



27ri/A(7)] YJ'jli + arcsin((7j)) 



(7.28) 
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Fig. 7.8. (a) Edge state spectrum using the edge magnetization weights (7.23) com- 
pared to (b) the phase space distribution of the classical weights (7.28) for the 
interior ellipse, [figure quality reduced] 
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see (6.58a) and (6.62). Like in the case of w^^ the classical weights w'^^ vanish 
as trajectories get almost detached from the boundary (since the numerator 
approaches 2^, while — > oo). 

Equipped with well-defined spectral densities for edge states and the corre- 
sponding trace formulas we can now proceed with a statistical and semiclassi- 
cal study of edge state spectra. We shall not only consider the statistics within 
a given edge state spectrum (in Chapter 8) but look also for cross correlations 
between different, classically related spectra (in Chapter 9). 



8 Properties of edge state spectra 

We apply the weighted spectral densities discussed in the previous chapter 
to analyze two aspects of the interior and the exterior edge spectra. First 
we perform a statistical analysis comparing the results to the predictions of 
random matrix theory. Second, we check the validity of the semiclassical trace 
formula taking the disk billiard as an example. 

8.1 Universal auto- correlations 

One of the central goals in the field of quantum chaos is to understand how 
the statistical properties of the quantum spectrum reflect the nature of the 
underlying classical dynamics [5]. We extend these studies to magnetic billiards 
by making use of the spectral measure of edge states introduced in the previous 
chapter. It was constructed to focus on the non-trivial part of phase space 
which is determined by the billiard boundary map (3.4). 

As a first point we check whether the edge spectra of both interior and exterior 
magnetic billiards display the universal characteristics of random matrix the- 
ory (RMT) if the corresponding skipping motion is hyperbolic. Our quantity 
of choice to characterize the spectrum statistically is the spectral form factor 
K{t). It is sensitive to correlations of the eigenenergies beyond the mean level 
spacing [107]. The standard form factor was already used in Chapter 5 to study 
the two-point correlations in the unweighted spectra of interior billiards. For 
edge spectra K{t) is readily defined in terms of the 2-point autocorrelation 
function of the edge density, 

^.oM = /Cge(^'+0^:^ge(^'-0^l(^'-^o)di.'. (8.1) 

Here, we included a normalized Gaussian window function gi to pick up a 
spectral interval centered at uq. 



115 



Before comparing to RMT it is advantageous to remove the trivially system- 
dependent properties of the spectrum by "unfolding" it [107]. This is a trans- 
formation of the spectrum which renders it dimensionless and of unit mean 
density. Dealing with a weighted spectrum the unfolding procedure must trans- 
form both the energies and the weights. The natural choice involves the smooth 
edge state counting function Nedge and the average weight (w^) / (w) in the 
spectral interval considered: 

i^n 7^Nedge(i^n) and Wn--^j-^Wn. (8.2) 

Here, the first and second moments of the weights, 

oo 

{w) = ^Wn 5'(Nedge('^n) " t'o) (8.3) 
n=l 

and 

oo 

{w'') = E 5(Nedge(i^n) - i>o) , (8.4) 

n=l 

are taken locally in the spectrum in terms of the window function gr (a nor- 
malized Gaussian of width ag.) As a. result of this unfolding, both the weights 
and the weighted density have unit mean. 

Since we are dealing with a discrete spectrum, the form factor must be aver- 
aged to be well-defined. The standard procedure is to take the spectral average 
over non-overlapping parts of the spectrum, 

K{t) = e'^''^ R,,{i>) g,{i>) di>j (8.5) 

as indicated by the triangular brackets. According to the spectral ergodicity 
hypothesis [129] this should be equivalent to an ensemble average for hyper- 
bolic systems. 

If we choose the widths of the Gaussians gi and g2 as Ogjypi and (7p\/2, re- 
spectively, the Fourier transform of the autocorrelation function leads directly 
to the power spectrum. The form factor is then given by the weighted sum 



where the Fourier transform of g is denoted by g. 



) , (8.6) 



The previous discussion holds for both definitions, (7.3) and (7.23), of the 
spectral density of edge states. However, it is necessary to keep the type of 
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Fig. 8.1. Form factors (8.6) of the interior (a,b) and exterior (c,d) edge state spectra 
for the asymmetric stadium (a,c) and skittle (b,d) bihiard at p = 1.2. The shapes 
are defined in Fig. 5.1. The functions fohow the RMT predictions of the GOE and 
GUE ensembles [107], respectively (dashed lines). The heavy lines correspond to 
stronger spectral averaging than the thin lines (cxp = 10 and ag = 3, respectively.) 
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Fig. 8.2. Distribution of the quantum weights Wn > 0.1 of the interior (shaded) and 
exterior (transparent) skittle spectrum at p = 1.2. The histograms show peaks whose 
positions are well reproduced by the phase space estimates (8.7) and (8.8) (indicated 
by the arrows). Unlike the interior case, the exterior distribution shows a tail due to 
the transitional states which ranges to the small weights. (For normalization (bulk) 
states with weights smaller than 0.1 had to be disregarded.) 

the underlying classical motion unchanged during the spectral averaging. This 
is conveniently done by taking the spectrum in the semiclassical rather than 
the conventional direction, see Sects. 2.3 and 3.3.2. In this case the quantum 
weights Wn are simply obtained by taking the derivatives with respect to A 
at fixed p, see Appendix A. 7. In the present section we use only this first 
definition of a spectral density of edge states since it is not possible to define 
magnetization-based weights for the semiclassical direction. 

Figure 8.1 shows the form factors for the interior (top) and exterior (bottom) 
edge state spectra for the asymmetric stadium (left) and skittle (right) bil- 
liard, respectively. The spectra were obtained in the semiclassical direction, 
at fixed p = 1.2, ie, for the same situation as in Fig. 5.4. The weights were 
obtained by numerical differentiation with respect to A. We observe that the 
interior form factors follow the RMT prediction of the Gaussian Orthogonal 
and Gaussian Unitary Ensembles, respectively, as expected from the specific 
symmetry properties of the Hamiltonians. In the interior case this is not sur- 
prising. To ensure essential hyperbolicity of the classical motion the value of 
p had to be chosen large such that the interior phase space consists only of 
skipping trajectories which cover it ergodicly. As a consequence, one expects 
that all interior states are edge states to an equal degree. Indeed, the interior 
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weights arc distributed narrowly around a mean value w, given by the ratio 
of weighted and unweighted mean densities, 



as can be observed from the shaded histogram in Fig. 8.2. The weights do 
not provide additional information in this case, which explains why K(r) 
reproduces the RMT prediction, like in the unweighted case. 

In contrast, the standard form factor - like any other standard statistical 
function - does not even exist for the exterior spectrum, which is dominated 
by infinitely many bulk states. Nonetheless, we find that the exterior spectrum 
closely obeys the predictions of random matrix theory (bottom row of Fig. 5.4) 
if viewed in an appropriate way, ie, by means of the edge state density. This 
way, a crucial test for the consistency of the spectral measure of edge states 
is passed. The quantum weights succeed to filter out selectively the relevant 
edge states, which in turn exhibit the universal characteristics expected for 
chaotic motion. 

The distribution of the exterior weights is given by the transparent histogram 
in Figure 8.2. One observes that the distribution of large weights is peaked 
like in the interior case. Again, the peak position is well described by the ratio 
of weighted and unweighted densities. 



with the mean unweighted density now given by the phase space estimate 
(3.15) of skipping states. Unlike the interior case, the distribution has a tail of 
transitional states which ranges to the infinitely many bulk states with small 
weights. 



8.2 The action spectrum 

We turn from the statistical analysis of edge spectra to their semiclassical 
description. Here, the main purpose is to show that the trace formula for the 
edge state density - which rates each periodic orbit with a classical weight - 
succeeds in approximating the exact edge spectrum. 

We choose the disk billiard for which an explicit periodic orbit formula is 
readily obtained from equation (6.111). For the exterior case and — R/p < 1 




(8.7) 




(8.8) 
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Fig. 8.3. Action spectrum of the exterior disk at p = 2R. The positive values give 
the Fourier transform (8.10) of the exterior edge density (absolute values). The 
positions of the peaks are well reproduced by the trace formula (8.9) (negative 
values) - except for the small peaks at integer t which are remnants of the bulk 
states. The peak heights match well in most cases; they are expected to fit better 
if a spectral interval larger than G [0; 48] is used. 

we find, see (7.8), 



^osc I 
"edge I 



2u 

7T 



oo 

E 



Ari/2 



E 



N=2 

sm(a) cosfa 



sm a 



X 



(sin(Q;) cos(a)) 2 



ir? sin( Ao9) / , , , TT 

^ cos 27ii^N aL(A(^) + A^- 



with a defined by equation (6.65). (The term | sin(Q! — ^)| corresponds to 
the normal component of the velocity, |nv|, in (7.10).) Moreover, the exact 
quantum spectrum of the disk is calculated relatively easily in terms of the 
roots of special functions, see App. A. 3. We calculated spectral intervals large 
enough so that the Fourier transformation of the spectral densities. 



^edge' 
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resolves the classical actions t of the underlying periodic orbits. Here, the 
function /i is a suitable window centered on the midpoint z/q of the spectral 
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Fig. 8.4. Action spectrum of the interior disk at a cyclotron radius p = 0.4 x R 
small enough to enable bulk states. The Fourier transform (8.10) of the interior 
edge density (positive values, v £ [0;60]) is well reproduced by the trace formula 
(negative values). Note that in the top part, which shows the remnant peaks of bulk 
contributions, the y-axis has a different scale. 
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Fig. 8.5. Definition of the Bunimovich stadium used in Sect. 8.3 and Chapter 9. 

interval. This action spectrum may be readily compared to the semiclassical 
prediction based on (8.9). 

Like in the previous section, it is convenient to take the spectrum in the semi- 
classical direction, at constant p. In Figure 8.3 we show the action spectrum 
for the exterior disk at a cyclotron radius p — 2R (positive values). The 
corresponding prediction of the trace formula (8.9) is given by the negative 
values. One observes that the peak positions match very well with the pre- 
dictions of semiclassical theory. The only exception are the small peaks at 
integer actions which are not reproduced semiclassically. They are remnants 
of the infinite number of bulk states. The peak heights are well reproduced 
most of the time, except if two peaks overlap too strongly. These deviations 
are expected to fade as a larger spectral interval is used and the the widths of 
the peaks decrease. This is also seen in Figure 8.4 where we present the exact 
and semiclassical action spectra of the interior magnetic disk - based on a 
large spectral interval {u e [0; 60] at p = 0.4 x R). Here, the cyclotron radius 
was chosen small enough for bulk states to exist in the interior. One observes 
again that the latter are very efficiently suppressed in the action spectrum 
giving rise only to the small peaks at integer values (shown in the top part 
of Fig. 8.4). In the Fourier transform of the unweighted density, in contrast, 
the bulk states obliterate the edge contributions such that not a single ac- 
tion is resolved (not shown). Blaschke and Brack [71] analyzed semiclassically 
the spectrum for the unweighted interior problem. The contribution of the 
bulk states was estimated and added by hand resulting in an unsatisfactory 
agreement between the semiclassical and the quantum spectra. 

In conclusion, we find that the the semiclassical trace formula succeeds in 

reproducing the quantum edge state density. It does so by weighting each 
periodic orbit contribution with a classical weight which vanishes for cyclotron 
orbits. This removes the bulk contributions analogous to - and consistent with 
- the quantum weights of the edge state density. 



122 




Fig. 8.6. Weights as obtained from the edge magnetization for the stadium bilhard 
at 6 = 0.2 and high energy. Note that the values of the weights differ in the (top) and 
exterior (bottom) by a factor of about ten. Notwithstanding, the mean edge state 
densities are equal to leading order. The classical cyclotron radius which corresponds 
to this part of the spectrum is large, p E [2; 2.32], giving rise to essentially hyperbolic 
classical motion. 

8.3 Using the edge magnetization 



Finally, let us demonstrate that the edge state density may as well be defined 
in terms of the magnetization as discussed in Section 7.4. 

Choosing the Bunimovich stadium billiard (defined in Fig. 8.5) we calculated 
the interior and exterior magnetization spectrum in the high-energy direction, 
aX h = 0.2. The selected spectral interval u G [100; 135] corresponds to large 
cyclotron radii p G [2; 2.32] giving rise to essentially hyperbolicP^ classical 
motion. Quantum mechanically, the problem exhibits one unitary and one 
anti- unitary symmetry (rotation by vr and reflection at one axis, respectively). 
Hence, the spectrum decomposes into two symmetry classes (a feature which 
will be used in the next chapter) while each class should obey the character- 
istics of the Gaussian Orthogonal Ensemble [107]. 

The weighted spectra are shown in Figure 8.6. Note that the weights are 
very different in magnitude, although they lead to the same average edge 

The term "essentially hyperbolic" means that although there might be small inte- 
grable parts in phase space their combined area is much smaller than the uncertainty 
product (6^7r)^. 
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Fig. 8.7. Form factor of the exterior Bunimovich stadium (Fig. 8.5) computed from 
the edge magnetization spectrum shown in Fig. 8.6, bottom part. 

magnetization (7.20). This is explained by the different areas ^kip of the 
interior and exterior skipping motion since the mean weight is asymptotically 
determined by the ratio 

= ^^^^ ^ ^ , (8.11) 

see (3.15) and (3.32); (in the interior case ^/skip = ^\ Similar to Fig. 8.2 the 
distributions of the magnetization weights are localized at u?*-^' (not shown). 

Figure 8.7 presents the form factor (8.5) of the exterior magnetization spec- 
trum restricting the energies to a single symmetry class. As one expects the 
form factors follows the GOE prediction (dashed line). This indicates that the 
weights (7.23) based on the magnetization succeed to filter the bulk states 
consistently. They perform as well as the weights (7.3) based on the boundary 
condition. 



9 Spectral cross correlations: The interior-exterior duality 



The previous chapter focused on the correlations within a given interior or 
exterior edge spectrum and their relation to the corresponding classical dy- 
namics. We now turn to a different type of question, namely, whether one can 
relate the interior and the exterior spectra belonging to the same billiard. 

From a spectral theory point of view there is no apparent reason why the 
spectra of the magnetic Laplacian defined on complementary domains should 
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have anything in common. However, it was shown in Section 3.2 that the 
classical periodic orbits of the interior and the exterior problem are in general 
intimately related. They come in dual pairs with equal stability and their 
actions adding up to an integer multiple of the action of a cyclotron orbit. 
Since the semiclassical spectra are determined by the sets of periodic orbits 
one expects that the correlation in the classical dynamics carries over to the 
quantum case inducing a relation between the interior and exterior spectra. 

This observation being made it comes as a surprise that one does not find 
any signature of a cross-correlation if the spectra are analyzed with standard 

statistical methods, ie, the interior and exterior spectra of a billiard seem to 
be statistically independent (not shown). Nonetheless, it will be found in the 
present chapter that there exists indeed a strong cross-correlation between the 
two spectra. It can be observed only if a quantitative definition of edge states, 
as developed above, is at hand. 

9.1 A semiclassical theory of spectral cross correlations 

In order to unravel the connection between interior and exterior edge state 
energies a special cross-correlation function is needed. It not only involves the 
Dirichlet energies of the edge states but also relies crucially on the information 
provided by their weights. 

The cross correlation function 

As the first step to obtain the appropriate correlator, we formally extend the 
definition of the edge state density to finite boundary mixing parameters A. 



The dependence of the spectral density on A will be needed only in the vicinity 
of the Dirichlet boundary condition, A = 0, (3.6), where an expansion to first 
order in A is allowed. The spectral density (9.1) can then be written only in 
terms of the Dirichlet energies and Dirichlet weights. 



4dge(i^;A) -— N(i/;A) 



(9.1) 



dedge(^;A) = ETr(A)<^(^-^n(A)) 



n=l 




(9.2) 
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which follows from (7.3) and the properties of the 5-function. The cross- 
correlation function is now defined as an integral over energy and boundary 

parameter 

Cii^o) =11 Cgf ^(^;A)<S"H^;-A)MA)5(^-^o)dAd^ , (9.3) 

with normalized Gaussian window functions h and g. Here, h serves to restrict 
the integration over A to the range where the linear approximation in (9.2) is 
valid and may have a width of order one. The function g is needed to regularize 
the pair distribution d°^^^^^^\i') dl^g^'^^^u) . It selects a narrow energy interval 
centered around the energy i/q and should have the width of a few effective 
nearest neighbor spacings. 

Inserting expression (9.2) the cross-correlation function turns into a double 
sum over the interior and exterior edge spectrum, 




C(^o) = E 77^9 h h - , (9.4) 



where the primes label the exterior energies and weights for the sake of brevity. 

This pair correlation function is far from being standard since it relies heavily 
on the weights attributed to the individual levels. However, this function is 
the most natural choice to accentuate the spectral cross-correlations which 
originate from the underlying classical interior-exterior duality. This will be 
shown below. The important point to note is that due to the small width 
of g only a few pairs of interior and exterior spectral points will contribute 
appreciably at a given i/q. It is the pairs with equal weighted distances from 
the left and right, respectively, to the reference energy i/q. Here, the energy 
differences are scaled individually by the reciprocal weight attached to each 
spectral point. The function h, in contrast, limits the absolute energy distance. 
Note also that the prefactor in (9.4) ensures that pairs which include at least 
one bulk state do not contribute appreciably to the sum. 

The term Cbg in (9.4) subtracts the background. It is approximated by 

<:^bg = 4dge(^^0) lE/^f^^^) -^edge(^o)) , (9.5) 



if we neglect the width of g and disregard the fact that the interior and ex- 
terior mean edge densities differ in the higher order terms. We shall discuss 
the correlation function further after we derive its main properties using the 
semiclassical approximation. 
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The semiclassical correlator 



We turn now to the semiclassical evaluation of the correlation function using 
the periodic orbit formula (7.2) discussed in Section 7.2.2. It apphes to com- 
pletely chaotic systems. One obtains a double sum over the skipping interior 
and exterior periodic orbits: 

C(i/o) = fdugii' — Uq) -^y^ — '■ r ' — ' r 

J ^ 7r2:^,r^|trM(7)-2|5 ry |trM(7')-2|5 

X I cos ^27rz/(A(7) + A(7')) - 7r(n^ + ny) - ^(//^ + ^y)^ 

+ COS ^27ri/(^A(7) — A(7')^ — 7i{n^ — ny) — -^(/^i ~ A*7')^ 

x/.(^El%v,l + iE|n;v;ij| (9.6) 

Here, h is the Fourier transform of the window function h, and the exterior 
quantities are again marked with a prime. The width of h is small compared 
to the sum over \nj-Vj\ (which is of order rij). As a result, the second term in 
the curly brackets of (9.6) is suppressed. In the first term of equation (9.6), 
h reduces the sum effectively to those pairs with approximately equal sums 
of angles of incidence J2j I'^j'^jl = J2j I'^'j^jl- The dual pairs of periodic orbits 
discussed in Section 3.2.3 have precisely this property. Hence, the only system- 
atic contribution to the correlator will come from these pairs. In Sect. 6.3.3 
we discussed the relations between 7 and its dual partner orbit 7, which may 
be summarized as 

Hi) + H~) =n^ = n- , 
trM(7) = trM(7) , 



7 ' 



(9.7) 



If wc retain only the contributions of the dual pairs the cross-correlation func- 
tion simplifies to a single sum over interior (or exterior) periodic orbits. As- 
suming global classical duality we obtain 

Ci'^o) = ^ E cos(27rn,(^o - |)) 9^ . (9.8) 

The restriction of the double sum (9.6) to the dual pairs is tantamount to the 
diagonal approximation used in the semiclassical evaluation of the autocorre- 
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lation function [108]. In the present case, the actions of the chosen pairs of 
periodic orbits complement each other to an integer n^, while in the usual 
diagonal approximation it is the resonant terms, A(7) — A(7') = 0, which give 
the dominant "diagonal" contribution. 

In deriving (9.8) the energy dependence of the amplitudes of the trace formula 
could be neglected since the variation of the energy was assumed to be small on 
the classical scale in (9.3). If uq is taken large (i.e. we are in the semiclassical 
regime of the spectrum) the classical quantities in (9.8) will hardly change 
as 1^0 is varied. By grouping together the contributions from all the periodic 
orbits with the same number of reflections we obtain 

oo 

<^M= _E /(n)^(n) cos(27rn(i/o-|)) , (9-9) 

with 

Assuming ergodicity, the weighted sum over classical n-orbits (9.10) can be 
calculated as a phase space average. For large n it takes on the universal value 
/(n) = n/8. At the same time, the number of reflections is geometrically 
bounded from bellow, n > rimin, for a given cyclotron radius. Hence, 

f{n) ^ \ ^ (9.11) 
as n > nmin ■ 

Equation (9.9) makes a clear prediction on the form of the cross-correlation 
function. Even if the classical dynamics changes slowly as uq is varied the 
infinite sum (9.9) will be appreciable only at energies uq = N + N E 
No, where the cosine terms are stationary. We therefore expect the cross- 
correlation function to display pronounced, equidistant peaks at large energies. 
Their positions are expected to coincide with the Landau levels (although they 
have nothing to do with bulk states) and their appearance provides a direct 
quantum manifestation of the existence of classically dual orbits. 

If the billiard exhibits a discrete symmetry the semiclassical theory suggests 
also a natural way to test that the predicted structures in C(fo) are not 
artefacts (eg due to the bulks states). In this case the cross correlation between 
exterior and interior spectra with different symmetries is derived in a similar 
fashion as (9.9). However, now we have 
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where counts the number of times the periodic orbit 7 crosses the symmetry 
hne [130]. Since in the sum will be even or odd with equal frequency the 
terms cancel on average and no correlation signal is expected. 



Action cross correlations 

Next, we consider the Fourier transform of the cross-correlation function (9.3) 
which highlights its fluctuating part. The semiclassical theory predicts a se- 
quence of equidistant 5-spikes at integer values, 

r 1 °° 

D{t) = / C{uo) e-2-'^°* duo^- J2 i-rf{n) g{n) 6{n - t) . (9.13) 

It ii'm.in 

They correspond to the sums of the actions of dual pairs, which complement 
each other to integer values, starting from the minimal number of reflections 

'^min- 



Using the edge magnetization 

Let us turn to the question whether the correlation is also seen if one uses the 
spectral density d'^^^ based on the edge magnetization (7.21) as the spectral 
measure. It has the advantage of being easier to measure both numerically and 
in experiments since one does not have to change the boundary conditions. The 
theoretical treatment is completely analogous to the above with the correlation 
function now involving an integration over the the magnetic length parameter 
b rather than the boundary condition: 

X g{u - uq) du . (9.14) 

The linear expansion of the dependence of the energies on P yields a double 
sum like eq (9.4) with the Wn replaced by the magnetic weights 'it;^'^' (7.23). 
Semiclassically, the integration over 5b'^ selects those pairs of interior and ex- 
terior orbits which satisfy 

«;;-^'T^ = w;-rVy , (9.15) 

ie, 

(^^ + ^^)p = 2(^^ + ^) , (9.15a) 

which is again the dual pairs, 7' = 7. (This may be seen from eqs (7.28) 
and (6.60) since the signs of the cxj (6.55) and the order of the points of 
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reflection are reversed as one goes from an orbit to its dual. Geometrically, 
it is evident that dual orbits with reflections satisfy ^ + .5^^ — '^'^pN and 

These predictions for the cross-correlation function are not restricted to purely 
chaotic dynamics, although the bouncing map was assumed to be hyperbolic, 
so far. For the (integrable) disk one obtains a completely analogous result. 
The function /(n) is not universal in this case, but the prediction remains 
that C(t'o) is peaked at the energies of the Landau levels. Below, in Sect. 9.3 
it will be shown that semiclassical correlations can be predicted even for mixed 
chaotic systems (without resorting to periodic orbit theory). 



9.2 Numerical evidence 

In this section we provide numerical evidence supporting the above semiclas- 
sical predictions. We start by analyzing the edge spectra of the ellipse billiard. 
The underlying classical motion is not completely chaotic but still we expect 
most of the predictions of the semiclassical analysis to hold. We choose this 
example because we accumulated the most extensive numerical data for this 
system. 

In Figure 9.1 we show the cross-correlation function (9.4) for the ellipse bil- 
liard at magnetic length b — 0.1. It was calculated from the edge spectra 
shown in Figures 7.3 and 7.6. The corresponding classical dynamics exhibits a 
strict one-to-one correspondence between the interior and the exterior classi- 
cal dynamics up to z/ = 21.6. Beyond this energy, when the cyclotron radius is 
greater than the minimum radius of curvature, the classical duality still holds 
in a substantial part of phase space. One observes that C(i/o) is strongly fluc- 
tuating and displays pronounced, equidistant peaks at energies uq — N + ^, 
as predicted by (9.8). In Figure 9.2 we focus on these dominant structures by 
plotting the cross-correlation function in terms of UshUt = z^o(iiiod 1) around 
one half. To check that the resulting correlation signal is not an artefact or due 
to the accumulation of bulk states we make use of the fact that the spectra of 
the eUipse decompose into two symmetry classes. As shown above (9.12) one 
expects C(fo) to be structureless if one correlates edge spectra belonging to 
different symmetries. This is clearly supported by the numerical results shown 
as a dashed line in Figure 9.2. 

The Fourier transform (9.13) of the cross-correlation function exposes the sums 
of actions of the contributing pairs of periodic orbits. The absolute value 
of D{t), calculated for the same spectrum as Fig. 9.1, is shown in Figure 
9.3. In this case, the periodic orbits of the de-symmetrized ellipse have at 
least nmin = 4 reflections. One observes that \D{t)\ displays distinct spikes at 
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3000 



2000 - 



1000 




Fig. 9.1. Cross-correlation function (9.4) for the elliptic billiard (eccentricity e = 0.8, 
h = 0.1, (Tg = 0.001, (Jh = 1, positive part). The pronounced peaks at = N + 

G No, indicate the existence of non-trivial correlations between interior and 
exterior edge states. (The figure remains unchanged if one removes all bulk states 
from the sum (9.4) by imposing a threshold on Wn] not shown.) 




0.4 0.425 0.45 0.475 0.5 0.525 0.55 0.575 0.6 

^shift 



Fig. 9.2. Cross-correlation function from Fig. 9.1 summed over integer shifts of the 
argument, c(t'shift) = Sn^("' + ^^shift); in order to focus on the positions of the 
peaks. In the double sum (9.4) the energies were taken within the same symmetry 
class (solid line) and between different symmetry classes (dotted line.) 
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20 



t 



Fig. 9.3. Fourier transformation D{t) (9.13) of the cross-correlation function C{uo) 
given in Fig. 9.1 (absolute value). The peaks at integer t correspond to the combined 
actions of dual periodic orbits (cr^ = 5 x 10~^, ah = 4.) 

integer values. The real parts of the peaks have signs (— )", as expected from 
eq (9.10) (not shown). As predicted by the scmiclassical theory the dominant 
peaks start at nmin = 4 which is a clear proof for the classical origin of the 
edge state correlations. The tiny peaks at t = 1, 2, 3 vanish if one decreases 
the width of the window function g (which in turn deteriorates the statistical 
significance of the result). They are due to the remnant contributions of the 
bulk states, and disappear if one removes the bulk states from the correlator 
sum by setting a threshold on the weights (not shown; the remaining peaks 
would not change by this procedure). 

We repeated the calculation of the cross-correlation function of the ellipse 



spectrum now using d[ 



(M) 
edge 



as the spectral measure. The resulting function 



exhibits peaks at the Landau energies similar to Fig. 9.1 (not shown). Its 
Fourier transform is given in Figure 9.4. Again, the peaks are located at integer 
values starting at t — A. This shows that the edge magnetization density niedge 
succeeds to unravel the cross correlations similar to the edge density dedge — 
a reassuring but not a surprising result. 

The ellipse spectrum considered so far exhibits generic, mixed chaotic dynam- 
ics with relatively large integrable parts in phase space. As the last point, we 
demonstrate that the correlations do exist also in a system which is completely 
chaotic. We choose the spectrum of the stadium billiard defined in Fig. 8.5 
and use the edge magnetization to define the spectral density. As discussed in 
Chapter 8, the spectral interval shown in Fig. 8.6 corresponds to cyclotron radii 
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Fig. 9.4. Fourier transform of the ellipse cross correlation function (9.14) for the 
ellipse defined in terms of the edge magnetization density. The graph should be 
compared to Fig. 9.3. (cTp = 5 x 10~^, ah = 0.5) 

large enough to ensure that the corresponding classical dynamics is essentially 
hyperbolic. Figure 9.5 gives the corresponding cross-correlation function. Like 
in Fig. 9.2 the variable is plotted modulo one in order to focus attention on the 
peaks. Again, we observe a clear cross-correlation signal for pairs within the 
same symmetry class (solid line) while the reference calculation from different 
symmetry classes shows no peak (dashed line). 



9.3 The pair relation 



The peaks in C(z/o) were attributed to the complementarity of the actions of 
dual orbits. Quantum mechanically their occurrence implies that there exists 
a pairwise relation between individual interior and exterior edge states. This 
follows from the discussion of the quantum correlator (9.4) above. We have 
noted that pairs of edge energies contribute only if they have the same weighted 
distance to the reference energy from the left and right, respectively. Since the 
peaks appear at z/q = + | the interior energies Ui and exterior energies i/- 
must appear in pairs which satisfy 



(9.16) 
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0.4 0.425 0.45 0.475 



0.5 0.525 0.55 0.575 

I^shift 



0.6 



Fig. 9.5. Cross-correlation summed over integer shifts of the argument hke in 
Fig. 9.2. The data belongs to the stadium billiard in Fig. 8.6 (i^ = 100 — 135, 
ag = 5 X 10~^, ah = 0.2, using the edge magnetization density (7.21).) A clear 
cross-correlation exists between energies of the same symmetry class (full line), 
while there is no signal if the energies are taken from different symmetry classes 
(dashed line.) 



w 



N+1/2 



N+3/2 



w' 



Fig. 9.6. For every correlated pair of interior and exterior edge energies, i' and ly' , 
there exists a Landau level + ^ such that the distances - scaled individually by 
the reciprocal quantum weights, w and w' - coincide. 



with integer N, see Fig. 9.6. Although this is not an exact relation, it will 
be the more precise the larger and the closer the two energies are, since the 
semiclassical approximation (7.2) and the linearization (9.2) then hold the 
better. 



It is clear from (9.16) that the information provided by the ratio of the indi- 
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top 


32.5367 


0.826 


0.0445 


0.0444 


0.506 


32.4775 


middle 


33.5248 


0.489 


0.0507 


0.0501 


0.533 


33.4733 


bottom 


32.5082 


0.286 


0.0288 


0.0248 


0.508 


32.4874 



Table 9.1 

Energies and weights of the correlated pairs in Fig. 9.7, with the primes indicating 
exterior states. 

vidual quantum weights plays a crucial role in unraveling this pair correlation. 
It explains why standard correlation functions, which involve unweighted den- 
sities, do not show any signal in general. Moreover, the fact that the quantum 
weights enter reciprocally in (9.16) explains how a pairwise relation between 
interior and exterior states can exist in spite of different local unweighted den- 
sities. It is consistent with the mean edge densities (7.4) being equal in the 
interior and exterior. Using the weights (7.23) which are based on the magne- 
tization we obtain the same formula. It implies that for correlated pairs the 
individual ratios of the weights defined by (7.3) and by (7.23) are approxi- 
mately equal. This is indeed observed numerically 

For a given interior edge state it is of course not known, a priori, which is 
the associated Landau level + | and the exterior weight. Therefore, even 
in an asymptotic sense it is not possible to infer an edge spectrum given 
the complementary one by just using the relation (9.16). However, having an 
interior and exterior edge spectrum available, one can decide whether they 
belong to the same billiard. In the spectra considered so far we could easily 
spot single pairs of edge states by just using the relation (9.16). Examples are 
given in Table 9.1. 

An alternative derivation of the pair relation 

An independent semiclassical derivation of the pair relation (9.16) can be 
obtained without invoking periodic orbit theory by inspecting the semiclassical 
map operators P''^* and P'^'^* of the interior and the exterior, see Sect. 6.2. 
We present here the derivation for the magnetization based weights (7.23), 
the calculation for the weights (7.3) is quite analogous. If follows from the 
discussion of (6.27) that the interior Dirichlet number counting function can 
be written in terms of the eigenphases d(,{v] fo^) of the unitary map operator 

pint 

NSlp(^; b'') = E 02. {e,{u: b') + tt) (9.17) 
e 
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where 027r(^) is the unit staircase function at integer multiples of 27r. The 
interior edge state density (7.22) is then given by 



Ci(-;&^) = -^^Nl-,(.;6^) 

= -77 E ^(^; b') 52. {9,{u; b') + tt) , (9.18) 

with 52.(0) = ^027r(6') the 27r-periodic (5-function. The unweighted spectral 
density (3.9) is obtained in the same way by taking the derivative of (9.17) 
with respect to u. Comparing the two densities we conclude that the edge 
state weights are given by 



(9.19) 

(y)=. mod 27r 



The same equations hold for the exterior quantities (labeled by a prime) with 
a plus sign in (9.18) and (9.19). Now we make use of the duahty relation (6.28) 
between the interior and the exterior map operators. It implies 

Qi{v; b^) + e'^{u; b^) = 27r (^i/ - ^) + 2nM . (9.20) 

with integer M. Take a pair of interior and exterior energies u and u' which are 
determined by the ith eigenphase, ie, 9i{i') +7c = 2nM and 9'^{i'') + 71 = 2-kM' 
with integer M, M' . If we expand^the eigenphases to first order around the 
Landau level i/ + = M + M'-M-|we obtain from (9.20) 



(9.21) 



with Ai/ = V — uq and At/' = z/' — z/q. The first equality holds if the distances 
to the Landau level Ai/ and Ai/' are sufficiently small. To the same degree of 
approximation we can replace the Landau energy z/q by the eigenenergies v and 
v' and using (9.19) we get the pair relation (9.16) in terms of the magnetic 
weights (last equality). The fact that it can be obtained without resorting 
to periodic orbit theory shows that the pairwise cross-correlation is a generic 
semiclassical feature of dual magnetic billiards and is not related to the type of 
the classical motion. We note that the duality relation (6.28) also imphes that 
in the vicinity of a Landau level the operators P™* and P''^* are approximately 
inverse to each other. Hence, for correlated pairs of eigenstates the normal 
derivatives at the boundary are expected to be approximately equal (see the 
numerical test below). 
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-2-10 1 2 



Fig. 9.7. Pairs of correlated interior and exterior wave functions. The energies and 
weights are given Table 9.1. (Ellipse billiard at 6 = 0.1; the shading is proportional 
to the modulus of the wave function, and the boundary is indicated by a dotted 
line.) [figure quality reduced] 
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Fig. 9.8. Dual pairs of classical periodic orbits in the ellipse billiard, at p = 0.57. The 
top and bottom orbits are stable, while the middle one in unstable. Their classical 
weights (7.11) correspond to the quantum weights (7.3) of the states in Fig. 9.7. 
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Fig. 9.9. Boundary functions |n| = — i^n^l of the correlated wave func- 

tions depicted in Fig. 9.10 (along one half of the boundary). Solid line: interior 
(v = 110.6567), dotted line: exterior (i/' = 110.4841). The difference is given as a 
dashed line. 

Correlated wave functions 

We proceed to present three pairs of correlated wave functions of the ellipse 
billiard. The interior states were chosen to have different locations in Figure 
7.6 which displays the distribution of quantum weights on page 108 (top part). 
At energies corresponding to p ~ 0.57 we took states with weights lying in 
the top branch of the rightmost bifurcation structure, in the middle, and in 
the bottom branch, respectively. The respective correlated exterior states were 
identified using the pair relation (9.16). Table 9.1 lists the data for the three 
pairs. 

Figure 9.7 displays superimpositions of the interior and the exterior wave 
functions. One clearly observes that the top and bottom wave functions are 
localized on dual periodic orbits. In agreement with this observation one finds 
that the structures of increased density of classical weights in Fig. 7.6 (page 
108, bottom part) may be attributed to periodic orbits which bifurcate as 
the cyclotron radius p is increased. The top and bottom wave functions were 
taken from a fork which belongs to orbits with period 6. Two of these orbits 
are shown in Fig. 9.8 along with their dual partners. 

The middle wave functions in Figure 9.7, in contrast, are localized on a chaotic 
region in phase space confined by un-destroyed invariant tori. For comparison, 
an unstable pair of dual classical orbits from this region is given in the middle 
part of Fig. 9.8. Note that it exhibits the same spatial extension as the wave 
functions. Here, the correlation of interior and exterior wave functions is not 
evident from the visual inspection. 

The most convincing evidence for the correlation between interior and exte- 



139 



-2 -1 1 2 

Fig. 9.10. A typical pair of correlated wave functions taken from the spectrum 
displayed in Fig. 8.6 (superimposed; v = 110.6567, z/' = 110.4841, h = 0.2). The 
stadium-shaped boundary is not drawn but visible as a regular nodal line. Figure 
9.11 shows the pair on a larger scale, [figure quality reduced] 




Fig. 9.11. The pair of correlated interior and exterior wave functions from Fig. 9.10 
on a larger scale. The black line indicates the stadium-shaped boundary. The circular 
scars in the exterior wave function match the classical cyclotron radius p 2± 2.10. 
[figure quality reduced] 
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rior cigcnfunctions comes from a striking similarity of their respective normal 
derivatives at the boundary. This is a feature shared by all pairs of correlated 
wave functions, including those which are based on a chaotic part of the phase 
space. To show this clearly we compare in Fig. 9.9 the normal derivatives of a 
typical pair of correlated states taken from the stadium billiard. One observes 
that the moduli resemble each other, even though they arc rather irregular. 
The difference of the interior and the exterior values is indicated by the dashed 
line. As a consequence of (9.16) one expects in general that the difference is 
the smaller the closer two edge states are in energy. If the edge energies hap- 
pened to coincide this would have to take place on a Landau level according 
to (9.16) and the continuation of one wave function would simply yield the 
other, ie the normal derivatives would coincide. 

In Figure 9.10 we show the wave functions of the pair from the stadium billiard 
for which the boundary functions are shown in Fig. 9.9. Although they exhibit 

the typical irregular pattern of wave function based on a chaotic part of phase 
space, one can notice that the interior and exterior structures match. In the 
exterior wave function scars of periodic orbits may be discerned if viewed with 
some measure of imagination, see Fig. 9.11. Here, the circular structures match 
the classical cyclotron radius p = 2.10. 



10 Conclusions 

The main goal of this review was to present a practical and coherent theoretical 
framework for the study of magnetic edge states allowing in particular the 
investigation of their semiclassical properties. An important step was to set 
up a quantitative and physically meaningful definition for the notion of edge 
states. It became the crucial ingredient for identifying the quantum analogue 
of the classical interior-exterior duality. In the remainder we discuss a few of 
the unsettled problems and comment on possible new directions of research 
motivated by the progress made so far. 

The boundary integral method presented in Chapter 4 was formulated only for 
smooth boundaries. This restriction was crucial for the regularization of the 
hypcrsingular integral operators, which is necessary at non-Dirichlet boundary 
conditions, and a general treatment of corners or cusps is still missing. Also 
the effect of the general boundary conditions on the smooth number counting 
function was computed indirectly only. Balian and Bloch and others developed 
a systematic method for the computation of these quantities [75,76] which 
should be extended to the magnetic case. The implementation of this program 
is technically far from being simple and awaits a proper treatment. 

As a natural direction of further research one should ask whether the proposed 
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spectral measure is applicable and useful in other areas, specifically for the 
physics of the (fractional) Quantum Hall effect, where the concept of edge 
states is frequently employed without a clear definition. The introduction of 
confining walls which extend to infinity adds an essentially new dimension to 
the study of edge states. It requires a scattering theory appropriate to the 
continuous spectrum encountered in such systems. The models investigated 
so far dealt with the simplest geometries [131,132]. Once an obstacle (ie, a 
billiard) is placed near a wall or between parallel walls, scattering resonances 
will appear. Their analysis is important from an experimental point of view 
and theoretically challenging since they introduce resonant transitions between 
counter-propagating modes along adjacent leads. In the field-free case the scat- 
tering matrix is often discussed in terms of boundary integrals, and we expect 
that the techniques developed in Chapter 4 may provide convenient tools for 
the corresponding magnetic scattering theory — the same holds for the devel- 
opment a semiclassical scattering theory. These problems are bound to find 
applications for magneto-transport experiments involving ballistic mesoscopic 
devices. 

The quantum interior-exterior duality belongs to a more general class of prob- 
lems emerging in the field of quantum chaos. Here, one considers systems 
whose classical dynamics are related in some way and asks for the quantum 
correlations implied by the classical correlations and how they emerge in the 
semiclassical limit. Recently B. Gutkin analyzed the boundary operators for 
a line partitioning the 2(i-sphere [133] and found that the corresponding semi- 
classical map operators are related by an analogue of Equation (6.28). More 
complicated surfaces or the consequences of other classical relations were not 
yet analyzed but are a natural direction of further research. 
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A.l Green function in angular momentum representation 

In Section 2.5 the magnetic Green function was obtained by a direct evalu- 
ation of the Fourier integral. In the present appendix we derive its angular 
momentum decomposition. This lets us correct some erroneous results in the 
literature and discuss the irregular Green function. Moreover, the solutions of 
the radial Schrodinger equation will be needed below, in App. A. 3. 

The symmetric gauge must be employed since only this choice renders the 
angular momentum a constant of the motion. In polar coordinates the inho- 
mogeneous Schrodinger equation (2.50) then assumes the form 

G. = -^<5(f-fo) . (A.l) 
An ansatz in terms of the difference of polar angles, 

1 oo 

G.(r,ro) = — E e^"^(^-''°)G„(f,fo) , (A.2) 

m=— oo 

which cannot be justified a priori, separates the radial and the angular coor- 
dinates. For f 7^ fo the functions Gm solve the radial Schrodinger equation in 
the plane 

Gm(r,ro)^0- (A.3) 
The definition 

G^(r,fo)=fHe-^"/2^^(^2) (^.4) 
leads to an equation for g^n, 

zaLi^) + (1 + M - z)g^{z) - - + ^ ^9m{z) = , (A.5) 

which is known as Rummer's differential equation and is solved by the reg- 
ular and irregular hypergeometric function, iFi and U, respectively [36]. For 
energies different from the Landau levels it follows that a pair of independent 
solutions ui,U2 of the radial Schrodinger equation (A.3) is given by 

u,{r) = fl"^l e-^"/2 iFi 1 + |m|, (A.6) 

and 

u,{r) ^ fH e-^'/^ uQ - 1/ + ~ 1 + |m|, f^) . (A.7) 



-^idl + ^d,) + ^{r + i^r^u 
4 r 4 r 



^2 J- o [r —m) 

dj + -df - - — z^r^ + 4i/ 
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Both are real valued solutions, ui is bounded at f = and diverges as f ^ cxo. 
The function U2-, on the other hand, decays like a Gaussian in this limit but 
displays a (logarithmic) singularity as f ^ 0. 

Another fundamental system of equation (A. 3) is obtained if one replaces ^2 

by 



~He+^V2uQ + i.+ H_t^,l + |^|,_f2^ . (A.8) 



This is a complex valued solution [134] which we call "irregular" . Apart from 
its logarithmic singularity at f — > 0, it diverges exponentially as f — > oo. 

Both Ui and U2 are needed to form a solution Gm of the inhomogeneous equa- 
tion (A.l) since the 5-function implies a discontinuity of the derivative, 

diGmih + 0, fo) - d,Gm{h - 0, fo) = . (A.9) 

^0 

The requirement that the Green function must vanish as f — > oo, together 
with its continuity at f = fo, leads necessarily to the form 

r< 1 \ui{f)u2{fo) ifr<ro 

Gm{r,rQ)^ ———I A.IO) 
roW^(ro) \u2{r)ui{ro) if r > ro , 

with Wronskian W = Uiu'2 — u'iU2- In total, the Green function in angular 
momentum decomposition and symmetric gauge is given by 

p / 1 I |m|— m 

'-'^^^'^o^- 47r^f^J \m\\ \h^) """Pi 262 

^ /I \m\ — m ^ 11 \ 
X iFi \2~^^~2 '^ + l"^l'^<j 

,l + |m|,^>j (A.ll) 

with Mk^n{z) and the regular and irregular Whittaker functions [36], 

and 

-,^J and 2;> := max ( -, . (A.13) 
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Note that this expression differs shghtly from the (incorrect) expressions in [32] 
and [26; eq (6.2.26)]. 

An independent solution to the inhomogeneous problem (A.l) may obtained 
if one drops the requirement that the Green function should vanish as f — > oo. 
It involves the irregular solution (A. 8) and leads to the Green function 

-^"'•^^ = ^ ^^^^^ &r - m 



m=— oo 



^ /I jml + m , , 
X iFi [2+^+^^ ,l + \m\,-z< 

^Vl \m\+m ^ . . \ / A -, ,x 

xU(^- + i/+^ ,l + \m\,-z>j (A.14) 

which we call "irregular" . This expression was derived by Tiago et al. [33] . Un- 
like the regular Green function (A. 11), this one diverges exponentially once the 
distance between initial and final point exceeds one cyclotron diameter. This 
property renders the irregular Green function impractical for most purposes. 



A. 2 The null field method 



The null field method is an alternative scheme to quantize magnetic billiards 
in the interior [33]. We include it for completeness although its practical use 
is limited. 

Let us start with equation (4.4). In terms of the irregular Green function it 
reads 

G(''^'^)(r;ro)9„/,^*— = 0, (A.15) 

where we chose ro G ]R2\P, Dirichlet boundary conditions, and the symmetric 
gauge. Rather than transforming this into an integral equation we place Tq on 
a (large) circle with radius Rp which is centered at the origin and surrounds 
the billiard domain. 

Now assume that the billiard boundary is given as a function r{9) of the polar 
angle and expand the unknown boundary function in a Fourier series, 

d^/,r{r{0))=Y.^'''ci. (A.16) 
e 

Using the angular momentum decomposition (A.14) of the irregular Green 
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function equation (A. 15) assumes the form 



with 



and 



^ e-^°-a„S^,Q = , (A.17) 

t,m=—oo 



I I |m|+m 1 

^-^ — — '"H^ 



2t /rm\IH /r2(^) 



7o 



= / e ^ ■ ^ —TTT exp 



x,F.(l + .+ H+^,l + H,-i;g))d9 (A.19) 



/ ^ 
^0 



2 2 ' ' " 62 



^ ^"pv 262 



^ ,'1 Iml — m , I , r'^{9)\ ,^ ,. , 

XiFi ,l + |m|,— ^ d^. (A.19a) 



2 ' ' 2 ' 62 

Equation (A.17) holds for all polar angles ^o- For negative arguments the 
function U is known to be complex and non-zero [134]. Therefore we can divide 
by (irn for all Rp which leaves the condition for the existence of a nontrivial 
solution Q to 

det(S^^) = . (A.20) 

This is a spectral equation which was derived by Tiago et al. [33] (except for 
a misprint in their paper). 



A. 3 Exact quantization of the magnetic disk 



We briefly describe how to quantize the interior and the exterior of the mag- 
netic disk. As discussed in Sect. 6.5, the problem is separable in the symmetric 
gauge. For this choice the exact solutions of the free Schrodinger equation are 
given above, see (A. 6) and (A. 7). It follows that the interior and the exterior 
eigenfunctions of the disk are specified uniquely by their behavior at the origin 
and at infinity, respectively, and by the angular momentum quantum number 
m. Since the interior wave function (at energy u) must be regular at the origin 
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it has the form 



(A.21) 



For the exterior wave function, which must vanish at infinity, we have 



imi? 



\m\ 



e 2b2 u( _ _ + 



L^— ,1+H;^j . (A.22) 



Upon applying the general boundary conditions (3.6) at the disk radius r — R 
one obtains the spectral functions 



R 

Cdisk( i^;m,A, — 



with 



R 
'I 





\m\b\ 




2R ) 




m 


il- 







, 1 |m| — m , , R^ 



-, 1 + |m|; 



(A.23) 



^(a, 6; z) 



iFi(a, fe; 2;) for the interior problem 
U(a, &; z) for the exterior problem. 



(A.24) 



{djK indicates partial derivation with respect to the jth argument.) Unlike 
the semiclassical case (6.98), one cannot predetermine the radial quantum 
number here but has to search for all zeros at given angular quantum number 
m. The derivatives of the energies with respect to external parameters are 
given explicitely by derivatives of the spectral function like in the semiclassical 
case. For variations in the boundary condition we find 



di/ 
dA 



R dsK 



A=0 



diK 



(A.25) 



with the arguments of K like above. Similarly, the derivatives with respect to 
the magnetic length are given by the quotient 



i?2 d.K 



62 diK 



(A.26) 



dz/ 



We note the relation 

d62 b dA 

which holds in the semiclassical case as well. 



, ^R dp 



(A.27) 



A=0 
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A. 4 The stationary phase approximation 



The method of the stationary phase yields asymptotic expansions of integrals 
with rapidly oscillating integrands hke / g{x)e^'^^'^^^^Mx. One can show that for 
large u the leading order contribution stems from the stationary points of the 
phase /. After an expansion of the phase to second order around these points 
and the use of the Gaussian integral 



1 

f e'^'^'dx = (^^j ' eisgn(a)f ^^.28) 

one finds that for functions f,gE C°°(M) where / has a finite number of 
non-degenerate stationary points xj, ie, f'{xj) = 0, the asymptotic expansion 
reads [135] 

(A.29) 

as z/ — oo. For functions of an iV-dimensional argument, f,g & C°°{M.^), an 
analogous form can be found [136]: 

^(x)e2--/Wd^x ~ f-) ' E ^ e2W(x,)-i.,f _ (A.30) 

xj |det/"(xj)|2 

Here, uj gives the number of negative eigenvalues of the matrix f"{xj). 



A peculiar d -function 

As an immediate application the stationary phase approximation allows to 
show that the complex function 




(A.31) 



has the property of a one- dimensional Dirac 5-function, 

/4(0de = l (A.32) 
/ giO SeiO dC = g{0) (1 + 0(£)) as £ ^ 0. (A.32a) 

This follows from (A. 28) and (A.29), with z/ = 1/e and / = x'^/2, and is 
not easily proven otherwise. The product of (A.31) for the two Cartesian 
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components of the vector r yields the two dimensional (5-function 

= S{r - ro) , (A.33) 



hm — — -r - exp 



2eb^ 



which shows up in (2.48). 

A. 5 The product relation of the map operators 

We show that the relations (6.25) for the products of the interior and exterior 
map operators (6.24) hold semiclassically. They were needed to prove the 
factorization of the spectral function. Since possible saddle point contributions 
are excluded by the vanishing prefactors and the duality condition the only 
relevant contribution to the product stems from regions where the initial and 
the final point are close. Hence, we are allowed to replace the boundary locally 
by a circular arc. It follows that the expressions (6.72) - (6.75) derived for the 
disk billiards may be employed to show that the kernel of the product (6.25) 
acts like a 5-function. Assuming the angle between the initial and the final 
point 6ip = {s — so)/R to be small we find 

(p^pr ) is, so) = 

I b , /d2(27ri/as), , , d2(27ri/aL) , , 

^27ri!/ (as (-/p' - ipo ) +aL (</3o +<5v' - )) 



X e 



^i^e^'^^^/d^'^e^^-^^-l®^^'-^^^ ^^^'^^^ (A.34) 

TT i? J ^ dip' |l if Td < 1 . ^ ^ 

with Fd = R/{hv^), see (6.64). The dependence on 5ip was expanded linearly 

in the phase and neglected in the prefactor. The latter is cancelled by the 
change of the integration variable to Ml •= '^d^p^^^a,i^{ipQ — if'). Likewise, one 
finds for the second combination of interior and exterior operators: 

(pl ) (^, ^o) = - ^ e-- / d^' p- e-- / - ^ ' 

VPL PS y V , o; ^ ^ j |q if Td < 1 , 

(A.35) 

with := ■nd^QSi.s{^Q — (f'). The sum of the kernels assumes the form of 
a semiclassical 5-f unction once the integration is carried out. The ranges of 



e((^' - ^oM'p' -fo + S) if Td > 1 
1 if Td < 1 



d^(7raL). 



dip' dipo 



-1 
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If Td > 1: 

: 9^0 - ^ ^ ^^o - + ^ <Pq+'<P 

us: ^r|cos(^) . Td + ^rl -Td + ^rl . ^r|cos(^) 

If Td < 1: 

(fi' : (po — T^ ^ y'o - <^o + >■ (^0 + 71" 



^s: -^.r^^ > rd + -r2 -rd + -r2 . --r^ 



^2 



2 



Table A.l 

Ranges of the integrations in Equation (A.36). 



integration differ for weak and strong fields. They can be found in Table A.l. 
Setting Tp = 1 arcsin(l/rd) one gets 



(^int^ext I _int„cxt 
Ps Pl +Pl Ps 



TT R 



So) = 



-ir2cos(^) 



r? 



- ^ ' duse'""^^'^ ifrd>l 



1 p2 
I 

1 T^2 



if Td < 1 



TT 



1 sin ^) 



b 



(A.36) 



This proves the identity (6.25). In a similar fashion one finds that the product 
(6.26) does not contribute semiclassically. 



To show the semiclassical unitarity of the interior map operator, pi"t(pi"t)t — 
id, we note the kernel of its adjoint explicity. Using (6.2), (6.14), (6.20), (6.21), 
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and (6.24) we find 



{^^')\s;so)^{^^'y{so;s) (A.37) 

= -e-^""'""^-^ I -(vsn)(g,go)^ e(-vsn)e^"''^'^s(5;3o)eix-ixo 
(27ri) 2 I (sin(Q;) cos(q;)) 2 

+ n)(5,go)^ e{-VL n) e-'i e'^^'^^^i^^^'^") e^*-'^° 

(sin(Q;) cos(q;)) 2 

If the integral corresponding to the operator multiplication pint^pmt-^f jg eval- 
uated semiclassically one obtains a finite contribution only if the initial and 
the final points are close, like in (A. 34) above. We may again replace the 
boundary locally by an arc of constant curvature and we have in this case 
— vg/Ln(s;so) = VgyLno(s; sq). Comparison with (6.24) shows that the oper- 
ator corresponding to (A.37) assumes the form (P™*)t — _g-27rii/pext rpj^g 
unitarity of P"^* follows now immediately with (6.28), and the same holds for 

pext 



A. 6 The straight line with periodic boundary conditions 



In this appendix we discuss a model system which allows studying the transi- 
tion from edge states to bulk states asymptotically. In order to remove the ef- 
fects of a a finite curvature we deform the boundary F of a billiard to a straight 
line of length =Sf . In addition to the (mixed) boundary conditions along the 
straight line we prescribe periodic boundary conditions at the end points of 
the line and perpendicular to F. This is clearly no longer a billiard problem in 
its proper sense and there is no distinction between an interior and an exte- 
rior. Nonetheless, the classical and the quantum problem is well-defined, with 
a discrete quantum spectrum. This simple system permits discussing of the 
asymptotics of bulk and edge states in a straightforward fashion, see Sect. 7.3. 

The problem is separable in the Landau gauge (2.11) and may be solved anal- 
ogous to the disk problem in Section 6.5. Now it is the longitudinal canonical 
momentum (ie, the transverse component of the scaled center of motion) which 
is the second constant of the motion. Due to the periodic boundary conditions 
it is quantized, 

Cy nb , . 

= -m , (A.38) 

with integer m (here we put the boundary on the x-axis) . The transverse part 
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if of the wave function obeys 

v'\z) + (v-\z')^{z) = Q, (A.39) 

with z :— 2{y — Cy)/h. The semiclassical and exact solutions of this equa- 
tion yield spectral functions hke in the case of the disk, see Section 6.5 and 
Appendix A. 3, respectively. We report only the results. 



A. 6.1 Semiclassical quantization 



For given longitudinal and transverse quantum numbers, m and n, the semi- 
classical energies of skipping states are determined by the roots of the spectral 
function 



Cline ( ^> A> y 



TT f Tim h 

— -\- arcsin 



Trm b 



a 



line 



^\ ( 3\ 
i/, m, — - TT n - . 



h ) 



irm b 

(A.40) 



The phase shift 

Q;"'"'(i/,m,^) = arctan I A 



1- 



(A.41) 



is determined by the boundary condition A. 



A. 6. 2 Exact quantization 



Equation (A.39) is solved by the parabolic cyhnder functions. It follows that 
the exact spectral function has the form 



5 line 



i/;m,A, — \ 



D. 



- 27rm— 1 + A 



Trm— D„ 1 ( — 27rm— 



+ 



1 



D„ , 1 1 — 27rm 



(A.42) 



where is Whittaker's form of the regular parabohc cylinder function [36]. 



A. 6. 3 The uniform approximation 

As the most important point, we are interested in a semiclassical description of 
the situation when the corresponding classical trajectory is just detached from 
the boundary. Since the WKB approximation of the wave function fails close 
to the classical turning points we have to resort to a uniform approximation. 
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see eg [34]. It yields the asymptotic wave function in the whole region around 
one classical turning point, z^, — —I^Jv^ in terms of the (action) integral 



w\z) :- 



f 


1 ,2 


'dz' 









- sinh i 2 arccos 
2 



arccos 



(- 



V 



+ arcsin 



+ 



2V^; V207 



1 - 



(A.43) 

if 2; < -207 

if -2yjv < z < 2y/i), 



which we define to be positive for any z. In the uniform approximation the 
two independent solutions of (A. 39) are given (for —00 < z < 1\fv^ in terms 
of the Airy functions [36] 



w{z) 



\V — 4 



Ai - sgn(z - z^ 



-wiz, 



(A.44) 



and 



w{z) 



(P2{z) = N — Bi - sgn(2; - z^) ( -w{z 



The general solution may be parametrized by an angle e 
(f{z) = cos(q;„) (fii{z) - sin(Q;„) (fi2{z) 



'2' 2J 



(A.45) 



(A.46) 



This form is particularly convenient. By virtue of the asymptotic expansions 
of the Airy functions [36] we regain the WKB wave functions in both the 
energetically forbidden region, 



4-2 — 2 \z 



and in the energetically allowed one. 



i^iz) 



1 



— -T COsiw{z) -J-Oiu 



{z -C Zi,, ) 

(A.47) 

{z > z^.) . 

(A.48) 



Note the factor one-half in (A.47) which arises in a non-trivial fashion when 
connecting the WKB solutions of the two regions [34] . 

The eigenfunctions turn into bulk states once the longitudinal quantum num- 
ber m is large enough to leave the boundary in the energetically forbidden 
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0.1 



0.05 



-0.05 




ctan(A) 



Fig. A.l. Energy shifts of the bulk states (A. 51) for the 4th Landau level and pa- 
rameter values = 5.39724, b = 0.25, m = 14 ... 20) which allow the comparison 
with the right part of Fig. 5.11. As m is increased the energy shifts Ai/ get exponen- 
tially small (7.12) and the boundary mixing parameter for which there is no energy 
shift approaches the Neumann condition. 



region, 



Tvbm 



^>1. 



(A.49) 



Prom the uniform approximation (A. 44), (A. 46) we find that in this case the 
angle is determined by the ratio 



tan(Q;u) 



Ai 


((Hi) 


-A^g^-lAi' 


((H^) 


Bi 


((Hi) 





(A.50) 



By comparing the asymptotic expression (A. 48) of the wave function in the 
allowed region with that of a Landau state (which has no phase shift) one 
obtains the energy shift Ai/ compared to the Landau energy. 



Ai/(m, A) = — arctan 

TT 



Ai 


((H^) 


-A^g^^-lAi' 


((H^) 


Bi 


((H^) 


-AV^^-lBi'((H^) 



(A.51) 



Figure A.l shows the energy shifts for the fourth Landau level as a function 
of the boundary mixing parameter. Here, the value of =Sf/6 was chosen to 
correspond to the situation of Fig. 5.11 in Chap. 5. We observe that the 
bulk state behavior is reproduced qualitatively even at this low Landau level. 
A quantitative comparison of the bulk energy shifts (A.51) with a quantum 
spectrum is given in Sect. 7.3.1. 
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For quantum numbers m which put the boundary into the energetically al- 
lowed region (|g„, | < 1) the angle a„ is semiclassically given by the phase shift 
(A. 41) obtained above. For these states the energy derivative with respect to 
A reads 



dX 



A=0 



I + arcsin(g^) 



(A.52) 



It coincides with the limiting expression of a large disk if we set |c| 
in (6.114), (6.105) and let ^ oo. 



R + Cy 



A. 7 Scaled spectra 



We collect a number of formulas for spectra defined in the semiclassical di- 
rection. As discussed in Sections 2.3 and 3.3.2 those spectra are obtained by 
decreasing the magnetic length b at fixed cyclotron radius p (unlike conven- 
tional spectra where p is increased at fixed b). Since in both cases the spectra 
are noted in terms of the scaled energy u = p^ /b^ the superscript- (p) is used 
to indicate spectra taken at fixed p. 

Scaled spectroscopy [137,138] has the advantage that the classical dynamics 
remains fixed as the spectral variable is increased. This allows to ensure that 
the underlying classical motion is chaotic throughout the spectral interval 
(Sect. 5.1) and to extract classical actions easily by Fourier transformation 
(Sect. 8.2). 

However, one should be aware of the fact that the spectrum obtained this way 
does not belong to a single self-adjoint operator. Rather, a sweep through a 
family of operators (parametrized by an effective Planck's constant) is per- 
formed as the spectral variable is increased. Clearly, the energies are real and 
the eigenvectors are still proper solutions of the Schrodinger equation but the 
latter are not orthogonal. Moreover, it may happen that two energies coalesce 
and vanish as an external parameter is varied. 

Most of the formulas in the main part of this article hold as well for spectra 
at fixed p after the substitution b p/y^- In particular, this is the case for 
the spectral functions and the trace formulas which are to leading order in v 
The smooth number counting function (3.10), for example, reads 



N^") (^,) =^ ^,2 _ ^ ^ 1 _ (A.53) 



However, care is needed in the case of the spectral densities. The density of 



155 



edge states is now given as 

oo 

dlU^) = j:^i'^S{u-ui^^), (A.54) 



with the weights now defined at constant p, 

" ■ dA A=o' 
Here, we obtain the mean edge counting function 



List of important Symbols 

Most important: 



2 



cyclotron radius: p p 

scaled energy: = To 
magnetic length: b o"' 



(A.55) 



Latin symbols: 



^ area of the billiard domain — 

area enclosed by the trajectory of 7 (6.59) 

^e^kip area determining the phase space of skipping 

orbits, Fig. 3.6 

Ai(2;) Airy function [36] 

A (7) geometric part of the action of the periodic orbit 

7 (6.46) 

A(r) vector potential at arbitrary gauge (2.10) 

A(r) scaled vector potential at arbitrary gauge 

(A(f) = 2A{bi)/{Bb)) 

ALan(i'), Asym(r) vector potential in Landau gauge (2.11) 

(symmetric gauge (2.12)) 
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as(r; Tq), aL(r; Tq) geometric part of the action for the short (long) 

arc (2.56), (6.1), (6.2) 

Or, ol annihilation operator of right (left) circular 

quanta (2.18) 

B magnetic induction (B — 'V x A) 

B biUiard bounce map (3.4) 

Bi{z) Airy function [36] 

b magnetic length (2.14), (2.33) 

C{uo) cross correlation function (9.3), (9.4) 

c, c (scaled) center of cyclotron motion (c e R^) 

(2.22) 

V domain of the interior billiard (P C R^) 

D{t) Fourier transform of C(t'o) (9.13) 

Df;{z) parabohc cylinder function (Whittaker's 

form) [36] 

(i(z/) standard spectral density (3.9) 

4dge(i^), o^ed^gU^) spectral density of edge states (7.1), (7.8), (7.22) 

dedgei^): ^^Idgll^) smooth spcctral density of edge states (7.4), 

(7.24) 

^edge(^) fluctuating part of the spectral density of edge 

states (7.2), (7.10) 

^edge(^) spectral density of edge states in the semiclassical 

direction (A. 54) 

dskip(i^) smooth spectral density of skipping states (3.15) 

^skip(^) fluctuating part of the spectral density of 

skipping states (6.61) 

E (kinetic) energy 

E proper scaled energy {E = E/{t\oj) = 2z/), page 15 

f{n) weighted classical sum over n- orbits 7^") (9.10) 

iFi{a,b;z) regular confluent hypergeometric function [36] 

W generating function of the billiard bounce map 

(3.5) 
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G, ^(r; To) free Green function at energy v, with Tq the 

initial point (2.51) 

q(sc) (^r; ro) semiclassical free Green function at energy u, 

(2.58), (2.59) 

Gl{z), G^{z) gauge independent part of the (regularized) free 

Green function at energy u (2.61), (2.63), (2.64) 

G°*^'^'^) (z) gauge independent part of the semiclassical free 

Green function at energy u (2.60) 

g{z) normalized Gaussian window function, g{z) = 

(27r(72)-5 exp ( - -zV(2(7|)), with "small" ag 

g{t) Fourier transform of g(z) 

H magnetic Hamiltonian (2.8) 

H scaled Hamiltonian (2.36) 

h{z), h{t) normalized Gaussian window function, cf g{z), 

with "large" width ah, (and its Fourier transform) 

j(r) probability current density (2.32) 

K(r) form factor (8.5), (8.6) 

K{a, b;z) confluent hypergeometric function (A. 24) 

jSf circumference of the billiard domain (=Sf = |r|) 

length of the trajectory of 7 (6.58) 

L canonical angular momentum (2.12) 

jC magnetic Lagrangian (2.1) 

£ scaled Lagrangian (2.34) 

A^(z/) scaled magnetization (3.27) 

M.edge{^) scaled edge magnetization (3.28), (7.19) 

Aledge(i^) smooth edge magnetization (7.20) 

Mk,fi{z) Whittaker function [36] 

Mmax maximum winding number in magnetic disk 

(6.79) 

M(7) stability matrix of 7 (6.50) 

m angular (or longitudinal) momentum quantum 

number (6.87), (A. 38) (or else integer) 
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mmax, "^min maximum (minimum) angular momentum 

quantum number corresponding to skipping 
motion in the disk (6.88) 

me particle mass 

m(z/) scaled magnetization density (3.26) 

mcdgc(^) edge magnetization density (3.30), (7.21) 

medge(i^) smooth edge magnetization density (3.32) 

m°^'^(i/) fluctuating part of the scaled magnetization 

density (6.63) 

H normalization constant of the wave function 

N(i/) spectral number counting function (spectral 

staircase) (3.8) 

N(i/) smooth number counting function (3.10) 

Ncdgc('^) edge state counting function (7.5) 

Nedge(i^) smooth part of the edge state counting function 

(7.6) 

N°|gg(i/) fluctuating part of the edge state counting 

function (7.7) 

Nskip(i^) smooth counting function for skipping states 

(3.14) 

N^gJP(i/) fluctuating part of the counting function for 

skipping states (6.52a), (6.54), (6.83), (6.84), (7.9) 

n-y number of reflections in 7, page 79 

n normal vector of billiard boundary, pointing 

outwards (3.2) 

^int' ^ext angular increment of A^-orbit in interior (exterior) 

disk (6.78), (6.82) 

P semiclassical map operator (6.17) 

p(s. So) kernel of semiclassical map operator (6.20), 

(6.21), (6.24) 

Ps Birkhoff coordinate conjugate to s, page 27 

p, p (scaled) canonical momentum vector (2.2), (2.15) 

Q boundary integral operator (4.10) - (4.13), (4.20) 
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q particle charge {q B > 0) 

qm quantized relative distance of the center of motion 

from the boundary (A. 49) 

q(r; fq) boundary integral kernel (4.25) - (4.28) 

q(sc) j-j,. scmiclassical boundary integral kernel (6.10), 

(6.13), (6.15) 

R disk radius 

R scaled disk radius {R = R/h) 

r-y number of repetitions in 7, page 79 

r, f (scaled) particle position vector (r e M^), (2.15) 

s curvilinear coordinate on boundary {sj = s{rj)), 

(3.1) 

T Larmor period (T = 2tt/uj) 

Tcyc cyclotron period {T^yc = ^T) 

i scaled time {i — cut) 

t tangent vector of billiard boundary (3.2) 

\J{a,b; z) irregular confluent hypergeometric function [36] 

U(r;ro) free quantum propagator (2.47) 

V scaled velocity {v — v/{ujh) — 2p/h) 

V velocity vector (2.6) 

Vs, Vl normalized velocity vector at point of incidence 

for short (long) arc(6.7). Fig. 6.1 

Vg, Vl normalized velocity vector after reflection for 

short (long) arc (6.7), Fig. 6.1 

W classical action (time domain) (2.40), (2.42) 

Wk,^i{z) Whittaker function [36] 

Wm w!^^ quantum weight of state \ipn) (7.3), (7.23) 

Wj, wi^^ classical weight of orbit 7 (7.11), (7.28) 



Greek symbols: 
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a{r; Tq) relative distance of the initial and the final point 

(0 < a < f ) (6.5), Fig. 6.1 

«(r,;r,+i) (6.32), Fig. 6.2 

a\ phase shift (depending on boundary condition A) 

(6.96), (A.41), (6.112) 

/3(r;ro), /5°(r;ro) relative direction of the normal vector at 

incidence (reflection) (6.6), Fig. 6.1 

pj, /5° angles (6.33), Fig. 6.2 

r bilhard boundary (F = dV) (3.1) 

Fd relative radius of the magnetic disk (6.64) 

7, 7^"^ physical periodic orbit (with n reflections), page 

76 

e parameterization of the angular momentum 

(1< £ < 1) (6.105) 

C relative distance (2.54) 

T) index for type of arc, r] e {S, L} 

Q{x) Heaviside step function 

9 polar angle 

k{s) curvature of the billiard boundary at the point s 

(3.3) 

A dimensionless boundary mixing parameter (3.7) 

A boundary mixing parameter (3.6) 

/x-y Maslov index (number of conjugate points in 7), 

page 80 

u scaled energy (2.37) 

^{u) spectral function (4.21), (A.23), (A.42) 

^'^^^^i') semiclassical spectral function (6.98), (A. 40) 

$ radial WKB phase in the disk (6.92), (6.93) 

(p polar angle in the disk {(p = s/R) 

p cyclotron radius (2.33) 

p radius vector (p = r — c G R^) (2.7) 

(7j arc parameterization (—1 < Uj < 1) (6.55) 
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ag, ah width of the normahzed Gaussians g{i'), h{u) 

scaled time of flight of 7 (6.60) 

X(r) gauge field, page 9 

X(f) scaled gauge field (X = Xo = X(fo)) (2.35) 

\E'(2;) digamma function [36] 

i\) stationary wave function 

cUc cyclotron frequency {uj^ = 2^7 > 0) 

LO Larmor frequency (2.17) 
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